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Abstract \ suggestion that it might be possible to deseribe equilibria between ion exc hangers 


and electrolyte solutions in analogy with membrane equilibria (Part ID) has been followed up to 


determine the influence of some variables on the kinetics of ion ex« hange. It is proved that at low 


concentrations of the solution the mass transfer through a liquid film around the grain is rate- 


determining for the process, whereas at high concentrations the rate determining step is the 


diffusion within the grain (Boyn ef al.). It is predicted that during a given process the relative 


influences of these two diffusion mechanisms need not alwavs be constant 


The quantitative treatment of the process ina column in which no instantaneous establishment 


of equilibrium occurs is generally very difficult from a mathematical point of view. With iso- 


therms having a negative curvature the break-through curve will. howeve r, in the end assume a 


stationary shape, when all concentrations travel through the column at the same rate. 


The treatment of this asymptotic approach is comparatively easy; in order to compute the 


break-through curve it is only necessary to measure the mass transfer coeflicients of liquid film 


and grain; this may be done graphically. 


Résumé— Afin de déterminer Vintluence d'un nombre de variables sur la cinétique de léchange 


dions, ila été donné suite A la suggestion (Part ID) qu'il serait possible sous certaines conditions 


de décrire les ¢ quilibres entre Péchangeur dions et les solutions d’électrolvtes en anlogie avec les 


équilibres de membrane. Ila été prouvé qu’éa une basse concentration de la solution le transport 


de matiére a travers un film de liquide enveloppant le grain est déterminant pour la vitesse de 


marche du processus tandis qu’éA des concentrations élevées, la vitesse est détermines par la 
diffusion a Tintérieur du grain (Boyp et al.). Il est prédit que, durant un processus déterminé 


il n'est pas nécessaire que les influences relatives de ces deux mécanismes de diffusion soient 


toujours constantes. 


Le traitement quantitatif du processus dans une colonne of il n'y a pas d’établissement 


instantane d’équilibre est en général trés difficile au point de vue mathématique. Toutefois dans 


le cas disothermes a courbure négative la courbe de traversée aura finalement une forme station- 


naire, quand toutes les concentrations passent A travers la colonne A la méme vitesse. 


Le traitement de ces approximations asymptotiques est relativement facile: pour caleuler la 


courbe de traversée il est seulement nécessaire de mesurer les coeflicients du transport de matiére 


du film de liquide et du grain; ceci peut étre effectué par moyen graphique. 


Zusammenfassung — Die Vorstellung, dass man Gleichgewichte zwischen TLonenaustauschern 
und elektrolytischen Loésungen in Analogie zu Membran-Gleichgewichten (Teil ID) beschreiben 


kann, wurde zur Bestimmung des Einflusses einiger Veriinderlichen auf die Kinetik des Ton- 


enaustausches weiter entwickelt. Es wurde bestiitigt, das bei kleinen Konzentrationen der 


Ldsung der Stoffiibergang durch eine Fliissigkeitsschicht um das Korn ges hwindigkeitsbe- 
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stimmend fiir den Prozess ist, wihrend bei hohen Kozentrationen der geschwindigkeitsbestim- 


mende Schritt die Diffusion innerhalb des Korns ist (Boyd u.a.). 


Man kann voraussagen, dass 


wihrend eines vegebenen Prozesses der relative Eintluss dieser beiden Diffusionsmechanismen 


nicht immer konstant zu sein braucht 


Die quantitative Behandlung des Prozesses in einer Saéule, in der sich das Gleichgewicht nicht 


unmittelbar einstellt, erweist sich im allgemeinen als sehr schwierig vom mathematischen Stand- 


punkt aus. Wenn die Isothermen cine 
sich ciner stationiiren Lage annihern 


keit durch die Siule strOmen. 


Die Behandlung dieser asymptotischen Anniherung ist verhiltnismiissig cinfach 


Durchbruchskurve zu berechnen 


keitaschicht und des Korns zu messen 


In the preceding part [1] we 


under 


compare d the 


capacity of ion exchangers practical 
conditions with the capacity that can be cal- 
culated from equilibrium considerations. As 
kinetic factors slow down the process the actual 
capacities are lower than those cal- 
the 


article will deal with the kinetic theory of ion 


found by 


culation. Therefore, present part of our 


exchange. 


QUALITATIVE CONSIDERATIONS ON THI 


KINI 


Over 30 years ago DoNNAN [2], in a survey of 


rics oF lon EXxcuanat 


the theory on membrane equilibria, pointed in 
passing to the resemblance between the processes 
occurring on a membrane and those met with in 
ion exchange. The experimental results des- 
cribed in Part IL [3] tends to support Donnan’s 
remarks. It seems useful further to examine this 
resemblance, as it may prov ide an insight into 
the kinetics of ion exchange. 

Formally it may be said that three mechanisms 
may determine the rate of ion exchange, nam« ly 
(a) diffusion of ions from the solution through a 
liquid film to the exchanger surface and 
(b) diffusion of 


particles and (c) chemical reaction between the 


vice 
versa ; ions in the exchanger 
groups on the exchanger and ions in solution. 

authors thought 


from their experiments that the rate was deter- 


Some they must conclude 
mined by a chemical reaction inside the exchanger. 
Thus Nacnop and Woop [4, 5), like Jupa and 
Carron [6], correlated their experiments with a 
However, it is 


unlikely that chemical kinetics should play an 


bimolecular reaction equation. 


important part in ion exchange ; the rate of ion 
reactions in solutions is determined by diffusion 


negative Kriimmung haben, so wird die Durchbruchkurve 
wenn alle Konzentrationen mit der gleichen Geschwindig- 


Um die 


braucht man nur den Stoffiibergangskoeflizienten de Fliissig- 
Dies kann auf graphischem Wege geschehen. 


phenomena, and it is difficult to imagine why a 
similar rule should not apply to ion exchangers, 
in which the diffusion of ions is slower by one 
order than it is in solutions (Hom [7], Boyp and 
SOLDANO [8]}). 

The cause of the improbable results obtained by 
these authors may be that in the mathematical 
treatment of experimental results sometimes a 
chemical reaction cannot be distinguished from 
film diffusion (Boyp, Apamson and Myers [9], 
Apamson and Grossman [10]). 

At present it is assumed that only film and 
particle diffusion can be rate-determining factors 
and this assumption is borne out by practical 
results. There is an important difference between 
these two types of diffusion, in that diffusion 
through a liquid film depends, partly, on the 
velocity of the solution with respect to the particle, 
whereas particle diffusion does not. 

We shall now examine whether it is possible 
a priori to make predictions about the relative 
influences of film and particle diffusion in different 
cases. For the present we shall assume that both 
film and particle diffusion play a part, and that 
equilibrium — viz. a Donnan equilibrium — exists 
at the interface. The assumption of equilibrium at 
the interface implies a transition of ions from the 
the 
which is very rapid as compared with the diffusion 


solution into exchanger, and vice versa, 


through solution or exchanger. 


Then the number of ions transferred per unit 
of particle area is (cf. Fig. 1) 


Ny = ky a(a — = kp (q; — 9) (1) 


In this formula k, and kp are the transfer co- 
efficients in the liquid film and exchanger res- 
pectively; they both denote rates. The general 
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Fic. 1, Changes in concentration at the particle surface. 


applicability of the last term in the above equation 
may be doubted, the more since the particle is 
not flat but spherical, so that much more compli- 
cated diffusion equations are to be expected (ef., 
eg. Boyp et al. [9}). Recently, however, 
GiurncKAUF [11] demonstrated by calculation 
that the simple equation given above is a good 
approximation also for particle diffusion, and that 
it may be used for computing columns, provided 
the position of the equilibria is not too extreme 
(q is the mean concentration in the particle), For 
this reason use will be made chiefly of this approxi- 


mate equation. 


If the number of dissociated groups in the 
exchanger be g, the Donnan equilibrium at the 


interface may be expressed as: 


qj 

or 
(3) 


From equations (1) and (3) it follows that : 


(4) 
kp 


Prevalence of particle diffusion gives : 
N, ~ kp U (2) — (5) 


whereas, if diffusion through a liquid film is 
slower than particle diffusion, we get : 
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= kya (a — aq) (6) 


For exchange of bivalent for monovalent ions we 
may write :— 
Yi 
(l—ajPa 
If x, is small we may then use the approxima- 


aq 


tions «, ‘and 


ky ky 


Ifqg; g. the following equations hold : 


(1 z)vag q 
| g 
kp N@ ky, 
From (7) and (8) it is seen that the relative 


influence of a and g on the reaction mechanism of 


(8) 


bivalent-monovalent exchange is the same as 
with monovalent-monovalent exchange (ef. 
equation (4)). Also equations (5) and (6) hold for 
bivalent- monovalent exchange. 

From equations (4), (7) and (8) the following 


conclusions may be drawn : 


(a) The higher the concentration in solution, the 
greater the chance that particle diffusion is an 
important factor. This has been shown experi- 
mentally by Boyp et al. [9]. A theoretical example 
may further illustrate this point. For k,; and kz 
the following formulae may be given : 


ky, = (9) 
ky = Dp Olr (10) 
Equation (9) follows directly from the first law of Fick 
and (10) from GuurcKaur’s [11] equation : 
kK’ = 15D/r (11) 

in which, for our case k’ is equal to 


particle surface area 


Kp: (12) 
r 


bed volume 


If now Bp 0-5, equations (11) and (12) lead directly to 


equation (10). 


| | 
| ca | 
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q 
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If we further 
d—5-10* em"); r 
equiv it 
= 100 kp. 


assume 
0-05 em 
from (9) 


. 10; 
2000 m 
(10) that 


and 


‘ollows “ane 
foll 1 


get k, ~ OL 


so that, according to equation (4), film diffusion 
then 


If now a 2m equiv, |, we 


must dominate; if a 200 m equiv L., 


a 
ky LO kp, in which case partich diffusion is 


a 
~ 


the rate-determining factor. 


(b) Lf a strongly acid exchanger (U* form) is 


contacted with a buffer solution, a is not constant 
. at the 


throughout the process but varies from a 


start, to a low value. In buffer solutions any 


and 
liquid 


influence of particle diffusion will dis 
the influence of the 


diffusion through the 


film remains. 


(c) The greater a, the more chance that diffusion 
through a liquid film comes to be rate-determining. 
This cannot easily be verified through experiments, 
as it is, in general, nece ssary to work with different 
types of exchanger which, owing to differences in 
porosity, for exampl may render the inter 
pretation of the comparison extreme ly dubious: 
a difference in porosity causes a difference in the 
effective value of r (ef. equation (10)). Moreover, 
D, itself ce pe nds on the ce gree of cross-linking, 
on the capacity of the exchanger and. probably, 


also on the concentration in solution {12}. 


(d) Only for monofunctional acid 
exchangers, such as Dowex 50, can g be equalized 
to the total 


weakly acid exchangers g varies with the degree 


strongly 


capacity, For polyfunctional or 


of loading (“degree of occupation ") of the 
exchanger. If a weakly acid exchanger loaded 


with Na form 
rather the HX form) by means of hydrochloric 


ions is brought into the (or 
acid, g is at first very great and decreases in the 
the With polyfunet onal or 
weakly acid exchangers any of film 


diffusion favour of 


course of process. 
influ nee 
lost in 


will eventually be 
particle diffusion, 


(e) If a buffer solution is led over a column 


containing a polyfunctional exchanger (H* form) 


it is difficult to predict anything (a and g both 
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Pre- 


sumably there are many possibilities in this case. 


increase when break-through occurs). 
As far as is known in weakly acid exchangers the 
influence of particle diffusion is predominant, 
Experiments carried out by McGarvey and 
Amberlite IRC 50 


calcium bicarbonate solutions at least point in 


THOMPSON with and 


this direction, as will be explained in another 
paper. Conway, Green and Reicuenperc [14] 
draw the same conclusion from their experiments 
in which this exchanger was neutralized with 
We 


may therefore state that experiments with buffer 


sodium hydroxide (cf. also Bauman [15)). 


solutions have shown particle diffusion to play an 


mmportant part. 


(f) From what was said in (b) and (d) it is seen 
that in a given experiment the relative influences of 
film and particle diffusion may be 
The 


been found before. 


subje ct to 


variations. literature does not show this to 


have However in our kinetic 
experiments which are described in a forthcoming 


paper such cases have been found. 
Base or ADSORPTIVE 


2. Tue GENERAI 


PERCOLATION IN A CoLUMN, SLowED Down 
BY TRANSFER PHENOMENA 
If there is no equilibrium in the column the 


material balance applies 


Instead of the equilibrium equation we now 


(13) 


Ss | =4 0 


‘ 


get a kinetic equation of the following type : 


(4) 


\ combination of 


(14) 


equations (13) and (14), 
together with the limit conditions. must enable us 
to solve any given case. Such a solution, in which 
ec and g are given as functions of x and ft respec- 
tively, will, even for simple cases, be highly com- 
plicated mathematically and difficult to work 
We 


shall therefore try to find out if, perhaps, there 


with (Thomas [16], [17, 18}). 


are asymptotic solutions which, although limited 
in their applicability, render possible the kinetic 
interpretation of the experiments by virtue of 
their simplicity. 


Ws 
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Below we shall consider— unless otherwise 
a column for which at ¢ 0 and z > 0, 
0: this is to be loaded with a 


stated 
the quantity q 
solution of concentration ¢,. After a short time t 
the relation between qg and ¢ will be as shown in 
Fig. 2, in which the initial state coincides with 


the horizontal axis. 


fo 


Fic, 2. Kinetic relation q, (¢). 


This curve will rise with time, or, in other 

words, | ‘) is positive and | is negative. We 
cf 

shall express each of these differential quotients 

in a formula : 


dq aa a (15) 


If v, is the velocity at which a given liquid con- 


centration moves through the column, then 


dz = v, dt ; and hence : 
(dq), vy dt 4 dt (16) 
and : 


where v, is the velocity at which a given con- 


centration g moves through the column. 


During the loading of a column | | 0, and, 


wl, 


as 0, we may conclude that v, > 
Ali 
This conclusion may also be drawn from the 
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Mw 

. . . 

expression for | which can be derived in an 
tla 

analogous way : 


(v, — v,) (“) (18) 


q t 
ve 
where, for the reason mentioned above, . <@ 
t 
From v1, > U, it consequently follows that the con- 


centration front in the adsorbent will tend to catch 
up with the proceeding concentration front in the 
liquid and will eventually coincide with it. 

2.1. Linear isotherms. This case is completely 
analogous to the problem of heat transfer from a 
flowing medium to a fixed bed. For the sake of 
simplicity, calculations are often made on the 
assumption that the isotherms are straight, but 
in experiments this simplification can only seldom 
be tolerated (Boyn, et al. [19]). Therefore, it 
may suflice here to refer to Section 2.3, as with 
the isotherms showing a negative curvature, the 
front here tends to adopt a stationary shape 
although less strongly. 


2.2. Isotherms with positive curvature. 
(f"(e) > 0; ef. Fig. 3). 


whether the eventual value of g is equal to f (ce). 


The question is now 


Kinetic relation q, (c) if the isotherm has a positive 
curvature, 


Fic. 3. 


Although this seems to be so, it is difficult to 
prove. It may be remarked, however :— 
That also if eventually equilibrium should be 


reached, v, is equal to v,, as :— 


q 
) 
958/59 
|_| 
— 
fic) 
resp | 
| 
| re) 
| | 
q, | 4 
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™ 
vy (19) 


dq df(e) 
Cae), 


That data in the literature point to the pos- 


sibility of circumstances under which, practically 
speaking, equilibrium will be established. For 
instance, GLUECKAUI [20] compares the statically 
measured equilibrium between Cu*+ and H* ions 
for Zeo-Karb, with the curve that can be cal- 
culated from the continuous elution curve (Part 
IT [3]), assuming equilibrium to be established in 
the column. The equilibrium curve and the 
calculated curve coincided. 

Sugata, Bancnero and Wuarre [21] led sodium 
chloride solutions over Dowex 50 in the K* form 
(the K* ion is more strongly adsorbed on this 
exchanger than the Nat ion), and found an elution 
eurve which indicated the establishment of 
equilibrium in the exchanger. The experiments 
of Wueaton and Bauman [22], already mentioned 
in Part IIL [1], point in the same direction. 

It cannot be said that v, vy implies the 
existence of equilibrium. However, experimental 
data do point to a chromatogram of the same 
appearance as that obtained at instantaneous 


establishment of equilibrium, 


2.3. Isotherms with curvature. 
(f ‘(c) 0; ef. Fig. 4). For isotherms having a 


negative curvature it can be prove d that the maximum 


negative 


values of q are given by the chord passing through 
(0, 0) and (c,. f (e,)) (Pig. 4). 


fic) 
resp 


Fie. 4. Kinetic relation q, (c) if the isotherm has a negative 


curvature. 
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Assume that eventually (q," [c] < 0). The value of ve 
derived from the material balance is : 


ny is defined by 


oy 
If, finally, 0 = ,, then ( | . After sub- 
AL 
stitution of this equation in (20) and subsequent 
differentiation it appears that do, de would be posi- 
tive From a physical point of view this situation is 
unreal, in analogy to that discussed in Part 1 [23] 


for the instantancous establishment of equilibrium, 


Consequently, the eventual situation will be 


f 


(22 
ve fo 
from which it follows that : 
(ey) 
q (23) 


Here qg is a function of ¢ alone. 


In the end the concentration front in the liquid will 
coincide with the concentration front in the ad- 
sorbent ; this event is characterized by the fact that all 
concentrations have the same forward velocity from 
that moment. We are then concerned with a so- 
called stationary sharp front. 

If information is desired about the kinetic be- 
haviour of an adsorbent it is advisable to study the 
desorption at f{”' (c) 0 and the adsorption at 
I'(e) 0; this is due to the fact that in both 
cases a sharp front will eventually result, which : 
(a) can be much more easily treated mathematic- 
ally than a widening front and (b) in principle, 
cannot ever lead to a situation indistinguishable 
from equilibrium. 

This had already been found in the experiments 
by Susata et al. [21], mentioned above : when 
sodium chloride solutions were passed over 
K*- Dowex 50 (widening front). situation 
resembling the instantaneous establishment of 


» 
(*) 
VOL. 
9 
] O58 / 
| 4 
Jae) 
ax 
9, 
Co 
| 6 = 


equilibrium was obtained, whereas when 


potassium chloride solutions were led over 
Nat— Dowex 50 (sharp front) a kinetic effect was 
clearly observed. 

For this reason many investigators have 
studied the sharp front (GLurcKAUF and Coarers 
[24], Hiesrer and Vermeuien [25], Micuarns 
[26], [27], Vermeu.en [28], VermeuLen 
and Hurrman [29], Larrpus and Rosen [30] 
and Bappour, Go and Epstrei [31)). 

Unlike other authors we shall not attempt 
to express the equilibria in formulae. The 
quantity which must necessarily be known, the 
driving force, will be determined graphically. 

The discussions given in the previous paragraphs 
have been illustrated in the Table. This Table 
also shows the result of analogous observations 


on isotherms in which inflexion points occur. 


3% Tue Break-Turoucun Curve ror a 
STATIONARY Suane Front 

The object of this section is to give a description 
of the way in which mass transfer coeflicients can 
be determined from a stationary break-through 
curve. In this connextion it is shown how the 
driving force for film and particle diffusion can be 
determined. Furthermore it is proved that an 
inflexion point in the curve corresponds with a 
maximum in the driving foree. Finally, attention 


is drawn to characteristic differences between 


SHAPE OF 


VALUE g (= q; 
8 AS A FUNCTION OF @ CONSTANT 


(ACCORDING TO EQUATION 
3 


RELATIVE NFLUENCE &, 
C AS COMPARED WITH Kk, AS CONSTANT 
A FUNCTION F g 


TATIONARY AT 


D 


Table 1. Characteristics of ton exchangers 


INCRE ASING 


INCRE ASING 


FRONT AT ADSORPTION |” WIDE 
EQUILIBRIUM 
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particle and film diffusion from which, by using 
experimental data, it can be derived which of the 
two mechanisms is rate-determining. 

3.1. The stationary break-through curve in 
chromatographic processes in general. Substitution 


of equation (23) in equation (14) gives : 


oc c 
F 24 
With x = 1, this equation gives the slope of the 
break-through curve of the column. F (c) 


depends on the rate-determining mechanism. 
In case particle diffusion is rate-determining, 
equation (24) (see equation 5) becomes : 

kp a’ (f(e) q) (25) 

where a’ is the particle area per unit bed volume, 
oc 
| can be measured, and the driving force 


‘ 
can be graphically determined from Fig. 5; 


therefore ky can be calculated. 


A significant point in the break-through curve is the 
inflexion point. This corresponds with the maximum 


value of the driving force f (c) q: 


r Jt) 


for (see Fig. 5): 


(26) 


FIRST DECREASING, FIRST INCRE ASING, 


DECREASING THEN INCREASING [THEN DECREASING 


+ — 


+ 


FIRST DECREASING)FIRST INCREASING, 


DE CRE ASING THEN INCREASING |THEN DECREASING 


|STATIONARY IN LiKE O LIKE DI: 
LIKE OT LIKE D 


U<a@<@, @« 


(THEORY ON HEAT ADSORPTION DE SORPTION 

BEST SUITED FOR KINETIC |i, ; 
INVE STIGA oN ED) ) r 


_[|AOSORPTION 


DE SORPTION 


IL. 
/ SQ 
2, c fle ) 
Oo 0 ‘ 
= 
‘ ‘ ‘ i 
| 
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Fic. 5. Determination of the driving force 


diffusion. 


If film diffusion is rate-determining. equation 
(24) (see equation (6) ) becomes : 


The driv ing force 


ky a’ (¢ (27) 


, can be read from Fig. 6 
and k, can be calculated from the measured slope 


of the break-through curve by means of equation 


Here again an inflexion point in the break-through 
curve is determined by the maximum of the driving 


force 


(28) 


Concentration ¢, which satisfies equation (28), corres 


ponds with a concentration c, in the fluid (see Fig. 6.) 


Fic, 6, Determination of the driving force in film diffusion. 
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in particle 


VAN KREVELEN 


Determination of the break-through curves at 
S) and different 
a’ values must vield different results for particle 


different linear velocities V (v 


and film diffusion because k, is also a function of 
the Reynolds number. 

3.2. The stationary break-through curve in ion 
exchange processes in particular. lon exchange 


differs in some respects from chromatographic 
processes in general; for particle diffusion the 


equation becomes : 


F(a 


(29) 


whereas for film diffusion one gets instead of 


equation (27) the following relationship : 


k,a‘a (x 
(%) 


showing that the shape of the break-through 
curve varies in proportion with the ion concen 
tration in solution a (provided that the equili- 
brium curve and hence ,, 

Column design, in which the calculation of the 
break-through 


is not dependent on a). 


curves constitutes an essential 
step, requires that the values of kp and k, be 
known. Our experimental work therefore will be 
aimed at determining these quantities using the 
(30). This 


will be dealt with in Part V of this series of papers. 


basic equations (25), (27), (29) and 


NOTATION 


total ion concentration in solution 

particle area per unit of bed volume 
concentration in the liquid 

127! 
thickness of film L, 


number of dissociated groups on exchanger ML $ 


diffusion coeflicient 


concentration in the adsorbent at the moment 

when equilibrium is established ML ® 

as f (c) for the ion exchanger 

mass transfer coeflicient in liquid film 

mass transfer coeflicient in exchanger 

length of column 

number of ions transferred per unit time and 

MT! 


per unit of particle area 


concentration in adsorbent (mean value) 


as 4 
f(c) 
resp 
f(c)-@ 
A, 
Co 
| of x 
| 
VOL. 
9 
| 
J ) | AD 
Ley) (30) a 
| 
aise 
(27). 
for 
(Cg) 
‘o 
is 
a : 
fic) 
resp. i, a 
gtc) 
Na 
=z 
8 


radius of the grain 


column 


x equivalent fraction of a given ion in solution 


"Na 


(mean value) 


a 


W 


ind VAN Krevecen D 


Kiamer K 
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Bauman W Jon Evchange (Edited by 
J 


Tuomas H. A 
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N 
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s cross-sectional area of column 
time 
1 volumetric rate of flow of liquid L*T 
Uy, t rate at which concentrations gq and ¢ move 
through the column LT 
? distance to the point where liquid enters the 
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Kinetic measurements 


K. Kuamer, J. C. H. Linssen and D. W. van KreEVELEN 


Central Laboratory, Staatsmijnen in Limburg, Geleen, Netherlands 
(Recened 12 December 1957) 


Abstract— Kinetic experiments are described which were aimed at measuring the mass transfer 
coefficients of liquid film and particle diffusion. A graphic determination of the driving force of 
the kinetic process was used. Experiments made with a number of exchangers have shown that : 


(1) Sooner or later the front assumes a stationary shape. 


(2) If film diffusion is rate-determining the coefficient of mass transfer is in good agreement 
with the value that can be calculated from known dimensionless correlations (vAN KReVELEN 
et al.). 


(3) The rate determining mechanism may depend on the degree to which the charging of 
the exchanger has proceeded. 

(4) If particle diffusion is the rate determining step, two type of mechanisms may be dis- 
tinguished : diffusion of ions through the liquid in the gel and exchange in the position of ions 


fixed on the skeleton of the exchanger 


Resumeé— Des essais cinétiques sont communiqués, qui avaient pour but de mesurer les co 
efficients du transfert de masse par diffusion a travers le film liquide et a Vintérieur du grain. On 
a utilisé une détermination graphique de la force motrice du processus cinétique. Des essais 


effectués avec un nombre d'échangeurs ont montreé les faits suivants : 
(1) Tét ou tard le front prend une fornne stationnaire. 


(2) Si la diffusion a travers le film détermine la vitesse, le coefficient du transfert de masse cor- 
respondra bien a la valeur qui peut étre calculée a partir des corrélations connues sans dimensions 


(VAN KREVELEN et collaborateurs) 


(3) Le meécanisme déterminant la vitesse peut dépendre du degré de charge atteint dans 
léchangeur 

(4) Si la diffusion a lintérieur du grain détermine la vitesse, on pourra distinguer deux 
types de mécanismes & savoir: diffusion d'ions a travers le liquide dans le gel ou changement de 
position des tons fixés sur le squelette de léchangeur. 


Zusammenfassung Es werden kinetische Versuche beschrieben, die auf die Messung des 
Stoffibergangskoeflizienten in der Flissigkeitsschicht und der Teilchendiffusion abzielten. Die 
treibenden WKrifte des kinetischen Prozesses wurden graphisch bestimmt. Versuche mit einer 
Anzahl Austauscher hatten folgendes Ergebnis 
Friher oder spiter erreicht die Front eine stationire Form. 
Wenn die Filmdiffusion geschwindigkeitsbestimmend ist, liejt eine gute Chereinstim 
munjoor zwischen dem gemessenen Staffiibertragung-koeflizierten und dem Wert, der 
aus bekannten dimensionslosen Beziehungen berechnet werden kann (VAN Kreve Len 
u.a.) 
Der geschwindigkeitsbestimmende Mechanismus kann von dem Grad abhangen, wie dic 
Beladung des Austauschers fortgeschritten ist. 
Wenn die Teilchendiffusion geschwindigkeitsbestimmend ist, miissen zwei Mechanismen 
unterschieden werden: Diffusion der lonen durch die Fliissigkeit in dem Gel und 
Platzwechsel von Ionen, die in dem Skelett des Austauschers fixiert sind. 
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In Part IV it was shown how the coefficients of 
mass transfer through the liquid film and within 
the grain can be derived from the slope of the 
break-through curve. This paper deals with the 
determination of the above coeflicients on a 
number of ion exchangers. Such data are needed 


in column design (Part VI). 


PROCEDURI 


Kinetic experiments were carried out on ton 
exchangers in a column of 4-5 em?®* cross-section, 
variables: quantity of 


using the following 


exchanger (50 -— 200 ml), total ion concentration 
(0-01 and 0-02 N) and rate of flow 
2-5 and 51 hr). 


Both the method of recording the effluent con- 


(normally 


centration and the interpretation of the experi- 
ments require that the rate of flow be kept 
constant; consequently, special attention was 
paid to this point. From the literature [1, 2] as 
well as from our own experience we know that a 
constant flow rate is often difficult to achieve. 


was 


The rate controller described by Haus. [3 


used ; the principles of this are outlined in Fig. 1. 


Fic. 1. Principles of the rate controller 


The pressure differential across the capillary 
A is Ap. The shaded part of the apparatus is 
filled with mereury. The diameter of the wide 
vessel B is much larger than that of C minus E. 
As a result a variation in preliminary pressure 
hardly effects the height of the mercury column 
in B (as a consequence Ap is constant); on the 
other hand it produces an equally large variation 
in back pressure as the level D of the mercury in C 
corresponds with this variation. To ensure 
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proper operation the dimensions must be so 
adapted to the anticipated maximum pressuré 
variations that the tube E will always be im- 
mersed in the mercury. The accuracy of the 
apparatus ts determined not only by the pressure 
variations but also by the ratio between the 
diameters of B and of C minus E. Fouling of the 
capillary is prevented by mounting a glass filter 
before the inlet. 

The concentration in the effluent is determined 
by measuring its electric conductivity. To this 
end the resistance of a conductivity cell was com- 
pared with a fixed resistance in a Wheatstone 
bridge fed from a 220 V_ line via a_ voltage 
stabilizer. The current in the bridge was rectified 
and indicated on a millivolt recorder via a potentio- 
meter. Calibration was effected by the use of 
solutions of known concentration. The tempera- 
tures of the liquid, the liquid controller and the 
conductivity cell were kept constant within about 
0-1 C; this was done to keep the flow rate constant 
and to enable the conductivity to be properly 
measured, 

The exchangers examined were the normal 
commerical products. As a rule they were not 
separated into fractions, partly in order to ensure 
that the 
possible with those obtained in practice, partly 


results would conform as closely as 


because the homogeneity of the particle size 
distribution is merely of secondary importance as 
compared with the effect of the other variables. 

The specific surface area needed for the cal- 
culations was derived from the particle size, 
determined by screen analysis, and the fractional 
free space in column §. The latter quantity was 
always assigned a value of 0-5. An exact deter- 
mination of the fractional free space is not 
feasible : for example displacement with water- 
immiscible liquids yields values which depend 
on the character of the liquid used and moreover 
deviate considerably from those obtained by 
centrifuging the wet exchanger. 

In all the experiments the concentration profile 
was assumed to be flat. Disturbances may be 
caused by the following factors. 

(a) Axial molecular diffusion; considering the 

low diffusion coefficient in aqueous solu- 


tions this effect may be neglected. 


a. 
OL. 
952/59 
c 
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(b) Awnial “eddy” 
COATES |4 


diffusion; GLUBCKAUF and 


have shown that this does not 


ive rise to difficulties even when dealing 


with small particles. 


Axial mixing due to differences in density 
when the heavier liquid is above the lighter 
no effect of 


5): at the concentration used 


this phenome non could he cle monstrated., 


Shrinkage and = swelling: Byrxe and 


Larpipus [6] noted a deviation from. the 


flat concentration profile during contra 
tion. This may have played a part in our 
experiments with Permutit ES 

that the above workers found 
Ambertite 


whereas we 


It is remarkable 


the strongly basi exchanger 


IRA to 


observed shrinkage 


swell during loading 


Column experiments were carnmed out with 


thre types of exe hangers which are often used in 


water purification practice the polvfunctional 
cation exe hange r Dusarit \ KG. the we 


Band the pol 


ikly basi 
functional 


Wit! thu 


stuched the « 


anion exchanger Dowex 
anion exchanger Permutit ES cation 


exchanger in the H* form we onver 


sion of sodium chloride to hvadro hlor wid and 
with the exchangers the 
this 
pra tical use of these « xchangers in demineralisa- 


Part I\ 7 


investigation has been restricted to the Stationary 


anion adsorption of 


hvdrochlori acid: with the 


corre sponds 


tion processes. As stated in our 
sharp break-through front. 
front 


Dusarit 


It mav be asked whether such a may 


VAG 
equilibrium on these ion 
tA; 


must be 


occur in the above with 


because the Na I 
ox hange rs shows an 
Part II [8}). That 


answered in the 


pre 


inflexion 
this 


big. 
question 
allirmative follows from 


part IV 


shown 


siderations analogous to those given iu 


17 or from Part I [9] where it has been 


that in the case of instantaneous equ librium an 
isotherm with ff" (ce) 0 gives a discontinuous 
with 


rhe 


occur at 


form at 0 ‘ ‘ as is also observed 


rl 
having a 
front 


is defined by 


isotherms negative curvature. 


stationary sharp may then 


0 x x, whe re 


r 


f(x.) 
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The adsorption equilibria for hydrochloric acid 
of Dowex 3 and Permutit ES are represented in 
Fig. 10, Part IT this 


graph both isotherms show a negative curvature 


As can be seen from 
so that the primary condition for the formation of 


a stationary front has been satisfied. 


2. RESULTS 


The results are shown in Figs. 2-6. 


Fic. 2. Kinetic measurements on Dusarit VKG 


In Fig. 2 (Dusarit VKG) the curves obtained at 
100 hed 
51 hr) run parallel; this points 
The total ion 
affects 
the slope of the break-through curve, the influence 


and 200 ml volumes sodium 


chloride and 2 
to a stationary shape (Part IV [7 
consice rably 


concentration solution 


of the linear velocity being less pronounced. This 


a marked effect 


indicates that film diffusion has 


(Part IV [7}). 


Fic. 3. Kinetic measurements on Dowex 3. 


In Fig. 3 (Dowex 3) the curves for 44-3 and 


88:5ml bed volumes do not run parallel in 


> ) 
N 
j 
VOL 
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: 
f 
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“ 
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contrast to those for 88-5 and 177 ml. The latter solution concentration on the slope of the straight 


two curves therefore have a stationary shape. section is greater than that of the flow rate. 


whereas in Fig. 6 this effect is not noticeable. It 


will be shown that film diffusion predominates in 


| / Fie. 5 and particle diffusion in Fig. 6. 


3. DISCUSSION 


3.1. Investigation of the stationary character of the 


CUTTES, 


Kirst the question of whether the experimental 


J y conditions permitted the formation of a stationary 


break-through curve will be examined. 


The front is stationary, for example, when two 


Fic. 4. Kinetic measurements on Permutit ES 


break-through curves obtained with different bed 


depths, but under otherwise identical conditions, 


In Fig. 4 (Permutit ES) all the curves ob run parallel. If x, x, (see equation (12), Part 
tained with different bed depths run parallel I (9|) the difference in time between the two 


again. As will be demonstrated later the log ¢ curves then amounts to : 


or log (« c) vs. time curve for this exchanger is ‘ Te ) 
0 (ry Po), S [ f(x») } 


ai 


partly linear; in Fig. 5 the influence of the t,—t 


since, according to equation (22), Part IV [7]. 


/ 
WH / all concentrations have the sane velocity, which 
/ / is equal to that at instantaneous equilibrium. 

/ With x, x, (x z,) the difference becomes : 


Sf (x,) (3) 


The above formulae hold for the Na H* 
J / exchange on Dusarit VKG. In the case of 


hydrochloric acid adsorption oh anion exchangers 


Fic. 5. Kinetic measurements on Permutit ES, plotted on the following expression is valid : 


semi-log seale for determining ky 
t, (vy S| f (e) | (4) 


‘ 


When a and J have no effect on the shape of the 


break-through curve ie in the case of particle 


diffusion (Part IV [7}), either a or }’ may be varied 


instead of the bed depth. 


The criterion for the stationary character of the 


measured break-through curve is that the 


difference in time (which can be determined 


‘ graphically) at all concentrations be equal to the 


value found from equation (2), (3) or (4). 


Figs. 2, 3 and 4 show that the majority of the 


Fic, 6. Kinetic measurement on Permutit ES, plotted on experiments were carried out in the stationary 


semi-log seale for determining kp. region. 
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With Dusarit VKG (Fig. 2) the time difference 
the 100 200 ml bed 
(0-01 N sodium and 
2-5 | hr flow rate) calculated by means of equation 


between curves for and 


volumes chloride 1) 


(3) 1s 


(7, te) Sf 


ai 


(1) TOO (meguis 
2-8 hr 


10 (mequiv 1) 


taken from 


Experiments show (Fig. 


value for f(x.) has been 


The 
Fig. 4a, Part II [8} 
that the time difference is 2-8 hr with a spre ad of 


2) 


1 per cent, the exact value depending on the 
concentration at which the difference is measured. 

With Dowex 3* (Fig. 3) much longer retention 
times are required to obtain sharp stationary 


break-through curves. Moreover thev are con- 
sick rably wider than thos« required with Dusarit 


VKG and KS. 177 


mil hydrochloric 


Using columns of 
(OLN 


the difference in 


Permutit 
bed 
and 0-575 1 hr 


and volume 


acid flow rate) 


time calculated by means of equation (4) is 


(a, to) Sf (eq) 


O-O885 (1) 2130 (mequiv 1) 


32-8 hr 


10 (mequiv 1) O-575 (1 tu 


The difference read from the graph is 82-6 hr 
with a spread of 0-6 per cent showing that the 
curves are indeed stationary. 

The difference between the curves for 88-5 and 
44-3 ml bed volume derived by means of equation 


(4) Is 


16-4 hr 


0-0443 (1) 


; 2130 (mequiy | 

10 (mequiy 1) 0-575 (1. he 

distance 


the 


per cent, 


As can be seen from the graph the 
the 
average being 14-6 hr with a spread of 7 
It will be that the 

ml stationary 


unsuited for further examination. 


between curves varies considerably, 


clear therefore curve for 


bed volume is not and is 


*For the obtained at O575 1 hr flow rate 


2-1 


curves 


Linssen and D. W. van Kerevenen 


Permutit ES gives a stationary sharp front 
difference between the three 
100 200 ml 
volumes (0-01 N hydrochloric acid and 51 ‘hr 


rapidly. The time 


curves, obtained at 50. and bed 


flow rate), are constant. Using equation (4) we 


tind : 


O-1 (1) 1170 (mequiv 1) 


2.34 hr 


10 (mequiv 1) 5 (ihr) 


(1) 1170 (mequiy 1) 


1-17 hr 


10 (mequiy 1) 5 (ihr) 


The differences read from the yraph amount to 


2-30 hr 1-5 per cent and 1°23 hr 2 per cent 


respec tively 


3.2. Investigation of the rate of mass transport 
Dusarit VKG (see Fig. 2). The 
centration proves to have an appreciable effect 


total ion con 


on the shape of the break-through curve, the 
influence of the flow rate being somewhat less. 
This de pendence of the shape on the variable is 
characteristic of film diffusion (Part IV [7}). By 
measuring the slope of the break-through curve 
[7]) Ay was 


calculated and the value thus found compared 


and using equation (30) (Part IV 


with the one derived from known mass transfer 
The 


first 


correlations were 
considered, the KREVELEN 


KkEKELS [10], the second by Tuoenes [11] : 


correlations. following 


by VAN and 


a D, dD, (9) 
0-33 


Vd,p 


bed 


Bk,d, 


(1 B) D, 
Here a’ area volume 


particle per mil 


(110 em‘), and 


= diffusion coefficient of the Nat ion 
(taken at 1-0 10 * em? /sec {12}). 


The k, values calculated by means of equation 
(30), Part IV 


| 
J 


and equations (5) and (6) are listed in Table 1. 


, (for various equivalent fractions) 


The measured values of k, are in good agree- 
ment with those found from both equation (5) 
and (6). Errors of measurement may be due to the 


following causes : 


and, 
vol 
1952/ 
—t 
a 
: 
a. 4 
14 


Table 1. 


Values of k, (Dusarit VKG) 
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Exchanger volume (mil) 
Flow rate (l/hr) 
Initial concentration (N) 


10° em /see 


Ky 
x 
3-3 
x oo 3-7 
Average 


According to equation (5) 


According to equation (6) 


Part IV, equation (30) 


3-7 5-6 
3-4 39 
3 3-3 3-9 


A small driving force which consequently 
cannot be exactly determined. 

Non-homogeneous particle size distribution. 
The impossibility of determining the exact 


value of 8. 


In equations (5) and (6) the D, - value for 
sodium chloride was used. According to ADAMSON 
and GROSSMAN [13] the value of the diffusion 
coetticient should depend on the relative con- 
centrations and be intermediate between the 
diffusion coeflicients of the participating ions. 
However, and Hecrrenicu [14] believe 
that the effect of the differing ion velocities is 
counteracted by the formation of a potential 
gradient in the film layer, which partly represses 
the effect of the faster ion. The correction 
suggested by ScHLOGL and co-workers requires an 
extra term in Fick’s first diffusion law which 
would render the calculations very complicated. 
Considering the limited accuracy of our experi- 
ments and the relatively small difference between 
the diffusion coefficients, we believe there is little 
use in making these corrections; this view is 
supported by the results published by Se.ke 
et al. (15). 

Later in Part VI it will be shown that an inte- 


gration constant t, can be defined by 


(%,) (7) 


whose corresponding a, value is independent of 


the variables. The value of t, calculated by means 
of equation (7) is 2-8 hr for the first three curves 
of Fig. 2 and 5-6 hr for the fourth. 

According to the graph the corresponding value 
of «, is 0-6; however, the curve obtained at a flow 
rate of 51 hr gives a different value. This curve 
also ought to pass through the intersection of the 
other two. 

Dowex 3 (see Fig. 3)*. Here again the question 
arises which mechanism is rate-determining. The 
fact that the concentration in the inflexion pomt 
is relatively low leaves no doubt that particle 
Other 


workers, hay Ing inspected this process (loading of 


diffusion will be am important factor. 


weakly basic anion exchangers with hydrochloric 
acid solutions) arrived at the same conclusion 
[16]. A plausible argument can be made however, 
that particle diffusion need not be the only rate- 
determining mechanism. 

Supposing that both film and particle diffusion 
play a part, we may write (see equation (1), Part 
IV [7]): 

N, =k, — ¢;) 


q) (38) 


At first ¢,; is either very small or zero (Fig. 10, 
Part Il [8]) and the process may be entirely 
described as film diffusion. This agrees with the 
results of batch experiments by Kunin and 
Myers [17] who found that in the early stage of 
the process the take-up of hydrochloric acid per 

*In all the experiments the particle size of the ex- 


changer was 0-71 — 0-85 mm 


200 200 200 100 
2°5 
0-01 0-01 0-02 0-01 
5 
| 
5 
O-6 50 + O04 35 +02 +08 
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Not until 


en 


unit time ts directly proportional to ¢, 


a large part of the exchanger surface has 


loaded 


than zero 


will « value appreciably larger 


adopt a“ 


and will party le diffusion become 


noticeable. 
The values of k, and ky, at 177 ml bed volume 


rate O-S75 1 he initial concentration 


0-01 N HCI) caleulated by mx ans ol equat 
Part I\ \7 and (5)(D 


a 38-5 em 


(tlow 


') and equation (2. 
respec tively have been compiled in Tab y 
As already stated in Part IV a part le diffusion 


equation as given above is not completely correct 


According to (sit rin AL 1S it better il 


approximation for sharply curved isotherms 


would be 
the best 


be ing : 


In both equations 


kp a (11) 


r2 
Equation (11) therefore can now be substituted 
for equation (11), Part IV [7 

The 


9) and 


values calculated by means of equa 

tions 10) are also given in Table 2 

Table 2 shows that with low effluent concentra 
diffusion 


the particl 


the equation for film ipplies 


high 


diffusion equation Is valid, The value s obtained 


trons 


whereas at concentrations 


Table 2. 


Linssen and D. W. van 


with equations (9) and (10) appear to be more 
constant than those caleulated with « quation (25), 
Part I\ 


this exe hanger are 


considering that the experiments 


with well reproducible this 
cork lusion SCCTIIS rmussable 

The value of Dy, derived by means of equations 
(10) and (11) is 41 Ww” 


same value as the on reported by Conway ef al. 


sec, this be ine the 
of neutralization of the 
weakly basic one r Amberlit 


hydroxide, viz. Dy, 


19] for the analogous case 
with 
sodium cm*, sec. 
whether stationarity of the 


effects the 


To demonstrat 
caleulated value 
stated that 


break throuch curve 
ol the transter com Hlicient it should bn 
curve for 200 ml 


Part IV [7 


w he reas 


with « Pmequiv | the 
bed 
hey 


vields the much hitch r value of 2-8 10 em sec, 


volun equation vives 


1-4 em sec equation (5) 


Strictly speaking, if particle and film diffusion occur 
neither of the particle diffusion equations 


Ie) The 


sitle boy sich 
be 


wav of caleulating the Ky Value corresponding with 


used applied sine« only 


the Ay value found from equation (5) would then be 
by trial-and-error formulating Ky by 


Part IV [7] the 


determined vraphically (see Fig. 7) 


means of 


driving foree may 


llow 


method of caleulation used is 


equation (25 
again te 
ever us Ky ky the 


Permutit ES (see Figs. 5 and6). The expe riments 
with Permutit ES prove d to be less re produc ible 
The fairly consick rable 


than those with Dowex 3. 


shrinkage and swelling of the exe hanger during 


loading and regenerating probably plays a part 
therefore, to use 


here. It did not seem advisable. 


a more involved particle diffusion equation than 


Values of ky and k R (Dowex 3) 


10” (om) sec) 


Com 
N- 10° 


Kquation (27) 


Part IV [7 


equation (25) 


Part IV [7 Equation (9) Equation (10) 


: mis (27) 
[12]: 
VOL 
9 
‘ } 
= (9) 1958/ 
Ali = 
(f (ec) q) (3 4f(e) 1 
(10 
= 
» 1-4 
3 I-4 Is 74 


Driving force in the case of combined tilm and 


Fig 
particle diffusion, 


the one sugvested in equation (8). From equation 
(8S) and the course of the isotherm it appears that 
the first part of the process may be described as 
film diffusion (c, = 0); equation (27), Part IV [7] 


then becomes : 


oe 
12) 
With «, 0, and the front must 
satisfy equation (25), Part IV 17]. 
ve Co 
| kp a (f (eo) q) (€g—e) (18) 


ot), (€9) 


In the former case log ¢ is proportional to ft 
while in the latter case log (Cy cis proportional 
to f¢, \ few 
logarithmic scale, are shown in Figs. 5 and 6, 
In the inflexion point of the break through curve 


both equations (12) and (13) are applicable and 


experiments, represented on a 


reach their maximum value : 


ky a’ kp a’ (Cy (14) 


The lines in the figures have be« n drawn as far 
as ¢,; hk, andk, have been derived from the slopes 
of the lines in Fig. 5 and 6; the k, value so found 
has moreover been compared with ky, calculated 


from equation (5) (a’ t7em ' and D 2-2 


10° em? sec [12]). These values are given in 
Table 3. 

The Dy, calculated from equation (11), Part IV 
[7] amounts to 8 10 
approximately equal to that of the self-diffusion 


sec this value being 
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Table 3. Values 0; ky and kp (Permutit ES) 


Exchanger volume (ml) 
Flow rate (Il /hr) 
Initial concentration (N) 


ky (Fig >) sec) 


ky (equation 5) 10° (om /sec) 


10° (em seu) 


Dowex 2. viz. 


coellicient of the Cl ion in 


354 107 om?® sec [20], 


3.3. Some general conclusions. 


It may bn obse rved that t he calculated ky 


value avuTree asonably with those derive dl from 
dimensionless correlations. This supports the 
view that. in processes where Dusarit VKG in the 


H* form is being loaded with a dilute sodium 


chloride solution, film diffusion is the rate- 
determining mechanism. As previously stated the 
values found for D, are of about the same order 
of magnitude as those determined from similar 
experiments by other workers. The large differ- 


ence (factor 100) between the D,- values for 


Dowex 3 and Permutit ES remains to be ex 
plained ; it probably depends on the groups on 
the exchanger be ing dissociated or not. 


\ similar case has been described by Vermet LEN and 
HvUrrmMan [21] who examined the neutralization of 
Dowex 50 with various amines dissolved in water 
and in a number of organic solvents. In experiments 
using the latter the diffusion ev« ficient was found to 
decrease with dielectric constant of the solvent the 
lowest value being 0-001 of that obtained with the 
aqueous solution It may be expected that the 
devree of disso lation of the ex« hanger will decrease 
accordingly 

The comparative measurements by KuNI~ and 


Barry [22] of the time needed for the Ca* Na 


and K i 


weakly acid cation exchangers also furnish an indica 


equilibria to becom established on 


tion that the dissociation of the ex« hanger plays an 
important part the former equilibrium is reached 
much sooner than the latter 

The difference in D,, between the two types of 


anion exchangers may now be qualitatively 
explained as being due to a difference in the 
mechanism of particle diffusion. We shall assume 
to this end that the particle diffusion observed 


with Permutit ES is a diffusion of hydrochloric 


: k 
ney) 
g(c)| 
> | 2-5 2-5 
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acid through the liquid in the exchanger. 


Equation (13) then takes the following form : 
kpa’' —— Ac, (15) 
) 


gel fluid. 
From considerations on the Donnan equilibrium 


where Cr is the concentration in the 


it can be derived that Ae p must be of the order of 


q- 
With Permutit ES ¢ ¢ © 10% (ey f (ey); see 
Fig. 10, Part Il [8]}). in other words (see equations 


(13) and (15)) : 


kp a 10 kpa (Co c) 


and, as ¢ and (c, — ¢) are of the same order of 
magnitude 


(16) 


in analogy to kpa 15 Dp 


Part IV [7]). 


We may now write for the diffusion coeflicient in 


(equation (11), 


the gel liquid : 


D 
kp a’ liquid ( 17) 


and thus : 


wel liquid 


The D,-value thus calculated is 100 times the 
anticipated D-value (10 5 om* sec) showing that 
the assumed mechanism (diffusion of hydrochloric 
acid through the fluid in the exchanger) cannot be 
of position of Cl 


well-dissociated exchanger 


correct. This leaves change 
and OH ions in the 
as the only ‘ rplanation for this rapid transport, 
With B (Dy 210") Dau guia 
em® sec: hence it must be concluded that with 


contains 


Dowe X 


this exchanger, which practically un 
dissociated amine groups, the normal diffusion of 
acid through the fluid is 


rate-determining. 


hydro hlori er hange 


Haus R. A. J. Sei. Instrem., 1955 32 116 


GLUBCKAU! and Coares J. J. J. Chem. So« 


Byane E. B. and Lapipus L. J. Amer. Chem 
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It will also be clear now why the modified 


method of determining the equilibrium on 
Dowex 3 rapidly led to equilibrium conditions 
(Part II [8]); first, the exchanger was loaded at a 
high concentration, in consequence the con- 
centration in the gel tluid was also high, causing 


loading to proceed rapidly. 


NOTATION 


a total ion concentration in solution ML“ 
a particle area per unit bed volume L 
‘ concentration in the fluid ML * 
D diffusion coeflicient 
d, particle diameter 
f(c) = concentration in adsorbent at equilibrium 
k, = coefficient of mass transfer in liquid film LT 
kp, = coefficient of mass transfer in exchanger i 
Kp for definition see equation (15) LT! 
Ny number of ions transported per unit time 
and per unit particle area MTL? 
q concentration in adsorbent ML 
r particle radius L 
S — cross-sectional area of the column L? 
t time T 
t t ly /i reduced time T 
| linear velocity LT ! 
t rate of flow 
i distance to the point w here the fluid enters 
the column L, 
x equivalent fraction of a given ion in solution 
p fractional free space 
y ~ free space per unit length in column L* 
Viscosity ML 
p density ML 
Indices 


b = inflexion point 

integration constant 

i exchanger interface 

L = liquid 

R exchanger 

r = tangent point (see Part I [9]) 

condition before the liquid enters the column 
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Design and scaling up of ion exchange columns 


K. Kiamer and D. W. van KREVELEN 
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Abstract—This final article in our series on the characterization of ion exchangers deals with 


the design of columns. A comparison between calculated and measured break-through curves is 


piven, 
For those processes in which film diffusion is the rate-determining step the correspondence 


between these two types of curves proves to be very close. In the case of particle diffusion the 


agreement is less good, although still reasonable if allowance is made for the simplifying assump 


tion used in the calculations 
In the treatment of the scaling-up problem a formula for the kinetic break-through « apacity 


is derived. It is shown in what way laboratory tests have to be carried out if they are to vield 


results applicable to the design of industrial filters. 


Résumé Ce dernier article de cette série sur la caractérisation des échangeurs d’ions se rapporte 


iu caleul des colonnes. Des courbes ce 
na travers le tilm exteérieur est déeterminante pour la vitesse, 
as de diffusion a lintérieur 


traversee calculees sont comparces a celles mesurées. 


Pour les processus of la diffusi 
la concordance entre ces deux types le courbe s'avere tres nette ken 
des grains, la concordance est moins bonne, tout en étant assez raisonnable, vu la supposition 


simplitiante pour les calculs. 


wrandissement on dérive une formule pour la capacité de 


kn traitant le probleme de | 
les essais de laboratoire devront étre effectués 


traversee cinétique On montre de quelle maniere 
si Ton veut obtenir des résultats applicables au caleul des colonnes (filtres) tels quutilisés dans 


lindustric 


Zusammentassung — Dieser abschliessende Aufsatz in unserer leihe iiber die Charakterisierung 
Berechnung von Kolonnen. Berechnete und 


von lonen-Austauschern beschiftigt sich mit der 


gemessene Durchbruchkurven werden verglichen 
die geschwindigkeitsbestimmend ist, giht es eine 


Prozesse bet clenen 
der Teilchendiffusion 


liesen beiden Kurventypen. Im Falk 


sehr gute Ubereinstinumuny zwischen 
iber noch befriedigend im Hinblick auf die vereinfachende 


ist die Ubereinstimmuny weniger gut 


Annahme in der Berechnung 


Zur Frage der Massstabsvergrosserung wurde eine Gleichung fiir die kinetische Durch 


in welcher Weise Laboratoriumsversuche auszufiihren 


bruchkapazitat abgeleitet. wird 
Berechnung industrieller Kolonnen abzugeben. 


sind, um brauchbare Unterlagen fiir dic 


IN the previous articles it has been shown how — break-through curves on the basis of the data 


in the case of a stationary sharp available. 


it Is possible 
(2) To compare the calculated with the measured 


front to determine both the driving force graphic- 
break-through curves. 


ally) and the coeftlicient of mass transfer (Parts 
(3) ‘To establish some generally applicable 


IV and V 


The obj ct of this final article Is thre ( fold scaling up rules. 


(1) To give a discussion of the design of ion It is assumed that piston flow occurs. The 


possibility of a spread in residence times of the 


exchange columns and of the calculation of 
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liquid in the column owing to channelling has 


edt (¢o c) dt (3) 


phenomena of this type do not lend themselves d ‘. 


not been included in these considerations because 


to quantitative calculations. 
Substituting dt from equation (1) in equation 


Design or lon EXCHANGE COLUMNS (3). we get 


‘ 


the case of isotherms with a negative 
cd c) de 


In 


curvature a stationary front will eventually be : 
I (c) F (c) 


established: also if the isotherm shows an in 0 


flexion point, part of the front will eventually 
Both integrals can be caleulated as functions 


become stationary. These two cases will be dealt , ; 
of ¢ by graphic integration of ¢, F (c) and (eg — ce) 


with in the order in which they are mentioned ; 
I re spective ly. 


here. 
If the values of thes« integrals are plotted as 


functions of e¢, the pom|l of intersection of the 


(a) The isotherm has a_ negative curvature 


(ft 0). 


resulting curves gives the correct value of c,. 


Integration of equation (1) consequently leads to : 


In this case the slope of a stati mary break- 


through curve (ef. equation (24), Part IV [1]), is 


described by 


For the integration of this equation we choose 


From equation (5) the break-through curve 


c, and t, as limits; ¢, is defined by 


may be ealeulated. Starting from the point 


} 
t= + /, (2 
to every ¢ can be found by graphic integration 


of 1 F (e). 


One of the specifications to which a column ts 


(c,. 1). which ts now known, the ¢ corresponding 


This is the time at which —in the case of in- 


stantaneous establishment of equilibrium — the 


designed is a certain concentration in the effluent. 


= Cp, Which may not be reached before a time ty 


in which v, tp, volumes have passed through, has 


elapsed. Substitution of ¢,, and tp in equation (5) 


gives the value of t,.. The volume of the exchanger, 


LIS. can simply be calculated from equation (2), 


provided l, vos negligible, which is generally 


the case. At any rate equation (2) gives a relation 


between / and S. 
The forms assumed by equations (1), (4) and (5) 


in dependence of F (film or particle diffusion 


Fic. 1. Stationary sharp front are shown in Table 2. 


discontinuous concentration jump takes place (4) The isotherm shows an inflexion point 
(cf. equation (4), Part I 0), 
can now be calculated in the 


rhe value of 


The slope of the break-through curve is now 


‘following way : 
pomowing . described not by equation (1) but by 
Assume the break-through curve to have the 1 y 
shape shown in Fig. 1. The two hatched areas (“) l F (c) (6) 
must be equal dt}, (e,) 


. 

58/59 

9 
= 


while the integration constant ¢, is defined by 
t, = {ly + ISf' (e,)}/v (7) 

and ¢, by 

ede 1 c) de 

Fie), F (c) 


(8) 


Equation (8) can be derived by means of Fig. 2 


Fic. 2. Equilibrium curve and kinetic break-through 
eurve for an isotherm showing an inflexion point. 


Curve (1) represents the break-through curve in the 
case of an instantaneous establishment of equilibrium, and 
eurve (2) is the stationary break-through curve for the 
case when 0 < ¢ c,. It may be assumed that if « Ces 
these concentrations travel at the same rate as they do 
if there is instantaneous establishment of equilibrium. 
It then follows that if ¢« c, the curves (1) and (2) are 
parallel and that the area between these curves is equal to 
the area of the rectangle (c,. (c,. t,) (t,. cp) em). The 
hatched areas in the figure must then be equal, ie 


ty 


edt (ey dt (Cy Cg) dt (a) 


. 
0 te le 


from which follows after substitution of dt from 
equation (6) in equation (9) the equation (8) 


By the same procedure as described under 
(a) the determination of ¢, is now carried out. 


Integration of equation (6) produces : 


de 
(c) 


| 
| (t — t,) (10) 
ca 


f'(e,) 


— 


Using equation (10) we can design a column or 
calculate a break-through curve in the way des- 
cribed in (a) provided that in the effluent 
O< ¢tp<e 


K. Kramer and D. 
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The forms assumed by equations (6), (8) and (10) 
for film and particle diffusion are shown in Table 1. 

The above formulae apply to chromatography 
in general. For ion exchange analogous equations 
may be used. A summary of all the formulae is 
given in the Table. 


COMPARISON OF MEASURED AND 


CALCULATED Break-Turovucu CuRvVEsS 


An idea of the accuracy of column design may 


be gained by means of a comparison of cal- 
culated and measured break-through curves. In 


this section some practical cases will be discussed. 


Dusarit VKG 


Figure 3 shows the calculated and measured 


(a) 


break-through curves for the polyfunctional 
cation exchanger Dusarit VKG in the H* form, 
loaded with sodium chloride solution. 


—— apenmento 


wuloted 


| 


Fic. 3. Calculated and measured stationary break-through 
curves for Dusarit VKG. 


Data : amount of ion exchanger: 200 ml, concentration 
of sodium chloride: OO1N (ay 1), space 
velocity: 251 per hr, cross-sectional area of 
column : 4-5 em®, specific surface area : 110 em™', 
equilibrium curve ;: ef. Part Il [4], mechanism : 
film diffusion, so that 


oa x 
2 ky, a(x — a) (11) 
MJ, 
ky, = 8-0 10° em/sec 
(calulated by equation (19)) 


a, = 0-65 (calulated by means of the equation 
ad« a) dx 
(12) 


. eq ey 


wh 
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56 h (calculated by the equation 


jal y + ISf' (a,)/av}) (13) 


The agreement between the calculated and the 


measured break-through curve (as shown in 


Fig. 3) may be called very good. 


(hb) Dowex 3 

Figure 4 shows the calculated and the measured 
break-through curve for the regenerated weak 
basic exchanger Dowex 3 loaded with a solution 
of hydrochloric acid. 

As shown previously (Part V [2]), the mechan- 
ism of the process in this case depends mainly on 

— 


Expenmentol 
Calouloted 


6 


40 


Fic. 4. Calculated and measured stationary break-through 


curves for Dowex 2. 


particle diffusion. In order not to make our 
calculations too complicated, we shall use the 
simple equation for particle diffusion, in spite of 
the fact that at low concentrations the influence 


of film diffusion ts clearly perceptible. 


pa’ 
f(eo) Ne) — @) 


88:5 ml con 
OOLN, space 


Data: amount of ion exchanger used 
centration of hydrochloric acid 
0-575 1. per hr, cross-sectional area of 


velocity 
38-5 em™!, 


2-1 em?®, specific surface area : 
ef. Part [4) 


cm) sec 


column : 
equilibrium curve : 
kp: 69-107 
¢, : 63m. eq./l 

(calculated by means of equation (4)) 
t, : 82-8hr (ealeulated by equation (2)) 


The agreement between the calculated and the 
measured break-through curve (Fig. 4) is not as 
close as in the previous case. This is due to the 
neglect of the film diffusion and more especially 


exchange 


VI 


to the use of a very simple formula for the 


particle diffusion. 


3. Tue Scatinc-Up PrRoBLem 


Finally we shall consider the scaling up of 
laboratory columns to large-scale plant apparatus. 

Especially the break-through capacity will be 
considered, The kinetic break-through capacity 
W,.(%p) is the number of milli-equivalents of a 
viven ion taken up per litre of exchanger until 
the moment when a certain (low) concentration 


xp a of the ion is observed in the effluent ; it is 


ly 


From a technical point of view only the break- 


detined by 


through capacity on loading is of interest. Tech- 
nical loading is carried out nearly exclusively in 
the case of isotherms having a negative curvature. 
If the break-through curve is stationary, yw, is to 


be written — after substitution of equation (15) in 


the pertinent formulae of the Table — as : 


eq i 


and 
da ls 


(17) 
f av 


Wi 


Le 


In the case of particle diffusion as a rate deter- 
mining factor a technical-scale filter will give the 


same result as a laboratory filter, pro ided_ the 


group 
a’ lS ai a'lat (18) 


is the same for the two columns. It is essential that 


in the laboratory experiment the break-through 


curve is stationary. 

Bearing in mind that the integral in equation (17) has 
a negative value, we see from this equation that the 
break-through capacity decreases in proportion to the 
increase of the concentration in solution and of the space 
velocity (number of unit volumes of liquid transferred 
per unit time and per unit volume of exchanger, v/LS). 

McGarvey and THompson [5] found this relation to 


apply to the break-through capacity of Amberlite IRC 50 


in the case of calcium bicarbonate solutions. We may, 


j 
| | 
: 
4: - - 
—- +> 
23 


K. 


therefore, conclude that in this process particle diffusion 


is the rate-determining step 

VW ith film diffusion as a rate det rmining factor 
the matter is different. We that for 
Reynolds numbers below 10: 


may 


os “4 


k Vo 
p) 


and he nee 


ky (20) 


If the Revnolds number is greater than 50 we 


vet | 6} 


(3) 


(22) 


Similarity mav be expected in the case of film 
diffusion for Re 


constant 


and for Re 


50 uf 


a yrs constant 


only from the 


Not ol 


(break-through capacity) but also with r gard to 


view of loading 


regencration (percentage of exchanger revenerated 


and regeneration efliciencyv) the s ing up 


proble mis interesting. 
Naturally 
the isotherm 
Mathematically, the 
difficult) problem. 


we are here concerned wit! 1 CAN 


in 
which has a positive 

kinetics here 
It 
that in various instances (ef. Part IV [1 
the effluent 


instantaneous establishment 


iture. 
present a very 
however 


the « 


recalled 
if 


of equilibrium had 


curve was found to look 


ais 
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taken the it 


generally permissible in dealing with these cases 


place in column. Therefore. 


Is 


to use the equilibrium methods described earlier 


(Part III [7)). 


NOTATION 


total ion concentration in solution 
particle surface area per unit of bed 
volume 

concentration in the liquid 
diffusion coeflicient 

as detined by equation (1) 
concentration in the adsorbent at 
equilibrium 
concentration in the ion exchanger ; 
equilibrium 

coeflicient of mass transfer in the 
particl 

coeflicrent of mass transfer in the liquuid 
film 

length of column 
concentration in the adsorbent 
cross-sectional area of column 
tirne 

linear velocity of liquid 

th 


distance to the pout where the liquuid 


rate ofl ~ 


of liquid 
enters the column 
equivalent fraction of a given ion in 
solution 
free space per unit length in the 
column 
Viscosity 
density 
» (tp) kinetic break-through capacity 
Indices 
break-through 
discontinuous break-through 
equilibrium 
tungent pont 


condition before the liquid enters the column 
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Perforated-plate performance 


R. A. McAuustrer*, P. H. McGinnis, Jr. and C. A. PLANK? 


North Carolina State College, Raleigh, N.C. 


(Received 2 August 1957) 


Abstract—The effect of the ratio of plate thickness, T, to hole diameter, D,, on dry-plate pressure 
drop has been correlated over a wide range of column size, hole area, and Reynolds number 
through the holes. New data are presented showing the effect of T/ D, on the wet-plate pressure 


drop characteristics (for the air water system) of perforated plates in a 15-in. diameter column 


having a tray spacing of 18 in. An oscillating region, where the aerated liquid mass on the tray 


moves back and forth perpendicularly to the direction of liquid flow, is reported. The mass 


flow rate of gas at which oscillation begins (the oscillation point) has been correlated with the 


height of clear liquid on the operating tray. In every case the oscillation point occurred well 


before flooding. The data presented here agree with the weepage correlation of HuGuMark and 


C’CONNELL [1). 


Résumé—L effet du rapport épaisseur du plateau, T, au diamétre du trou, D,, a été relié a la 


chute de pression a travers chaque plateau sec dans un champ étendu de dimension de colonne, de 


surface de trou, et ce nombre cle Revnolds a travers les trous. Pour une colonne auuN plateaux 


mouilles (systéme air-eau) les nouveaux résultats presentés montrent leffet du rapport 7 dD, 


sur la chute de pression a travers chaque plateau. La colonne utilisée avait 15 pouces de diamétre 


et les plateaux perfores ectaient separes par une distance de 18 pouces Ila été montre que dans 


une certaine région le liquide ventilé sur les plateaux vibrait perpendiculairement a la direction 


de flux. Le flux massique de gas pour lequel les oscillations commencent (le point d’oscillation) 


a été relié a la hauteur de clair liquide sur le plateau opérant. Dans tous les cas le point d’oscillation 


ol sont en accord 


apparait bien avant que le plateau soit submergé. Les résultats présentés 


avee corrélation de HuGumMark et O'CONNELL [1] qui permet de déterminer le flux de gas pour 


lequel le liquide est retenu sur les plateaux. 


Zusammenfassung Die Wirkung des Quotienten aus Plattendicke T zu Lochdurchmesser D, 
wurde in Beziehung zum Druckabfall der trockenen Platte innerhalb eines weiten Bereiches 


der Kolonnengrésse, der Lochtfliche und der Reynoldszahl in den Lochern. Fiir eine Kolonne 


von 38 em Durchmesser mit perforierten Platten bei 46m Bodenabstand werden neue Werte 


mitgeteilt, die den Einfluss von 7/D, auf den Druckabfall der benetzten Platte beim System 


Luft-Wasser zeigen ks wurde gefunden, dass in einem bestimmten Bereich die beliiftete 


Fliissigkeitsmasse auf dem Boden senkrecht zur StrOmunysrichtung ins Schwingen geriit Der 


Massenstrom des Gases, bei dem die Schwingung beginnt (der Schwingungspunkt), wurde in 


Beziehuny gesetzt zu der Fliissigkeitshéhe fiber dem arbeitenden Boden. In allen Fallen wurde 


der Schwingungspunkt vor dem Fluten erreicht. Die hier mitgeteilten Daten stimmen mit der 


Beziehung von Hughmark und O'Connell [1] tiberein, mit der sich jener Gasstrom berechnen 


lisst, bei dem die Flissigkeit auf dem Boden zuriickgehalten wird (weep point). 


Tue number of technical papers appearing in the have many distinct advantages over conventional 


last decade or so [1-20] concerned with perforated- — bubble-cap plates. They are easy to fabricate and 


plate column performance is indicative of a lively are more economical than bubble- ap plates. 


interest in this type of counter-current liquid Mayrie.ip, Cuourcnu, Green, Ler and Rasmussen 


vapour contactor. The investigations have shown [16] and Jones and Py e [8] found that perfor 


that in addition to stable operation, these plates ated-plate efficiencies were greater than those for 


*Present address: Lamar State College of Technology, Beaumont, Texas. 


tPresent address: University of Louisville, Louisville, Ky. 
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bubble-cap plates under similar conditions. 
Keiiy [13] and 
forated-plate columns had greater capacities than 
bubble 


gradient 16}, 


ZENZ [20] indicated that per- 


those containing caps. In addition, 


hydraulic entrainment [8, 12], 


and pressure-drop [8] characteristics for sieve 
plates are, in general, superior to those for bubble- 
cap plates. 

The purpose of the work reported here was to 
obtain new operating data with special attention 
effect of the plat 


hole diameter ratio (7' D,), and on the oscillating 


focused on the thickness to 
region reported by McALuster and PLANK [17]. 

The 
having various values of 7 D,, was first studied 
extensive ly by Kamer et al. [10, 11 
and those of other investigators, indicate that for a 


dry plate pressure drop using 


plates 
Their results, 
and mass flow 


given plate thickness, fre space, 


rate of air, as the hole diameter is decreased the 
value of 
of the 


D,) 


pressure drop decreases until a 
dD, ~ 2:3 is reached. 
hole 


causes the pressure drop to increase. 


kurt he r ck creams 


diameter (increase in the value of 
It is shown in this paper that the increased 
2°3 
calculated 


dry plate pressure drop for values of 7 D, 
results from hole friction and can be 

Fanning friction 
Below 7 dD, 2-0 
negligible and the 


primarily of the 


accurately using the usual 
factor for smooth pipes [22). 
bn comes 


the friction term 


pressure drop is a function 


contraction and expansion losses when the gas 

Hunt, 
Hanson and Wiike [5]. All available dry-plate 
data (7'/D, ranging from 0-08 to 7-9) 
the additional 
ranging from 0-116 to 


passes through the holes as indicated by 


including 
iT D, 


well. 


results presented her 
5-333) correlate 

Studies of pressure drop as a function of \ apour 
different 
types of curves 
Fig. 1. 
1, is that for the dry-plate 


rate indicate several regions displaying 
The 


work are 


types of operation. three 


encountered in this shown in 
The curve of Type I, Fig. 
runs. No liquid was flowing in the column, and 


the pressure drop versus flow rate curve was 
linear on log-log graph paper, and had a slope of 
approximately 2-0. 

By far the most common type of curve with 
liquid running in the column is labelled Type I 


in Fig. 1. The column behaviour with this type 


and C. A, PLANK 


in. LIQUID 


dP, 


log 


TYPE | ORY-PLATE 
OPERATION 


TYPE WET-PLATE 
OPERATION 


log Gy -ft* 


Characteristic operating curves for perforated 


plates. 


of curve is discussed in detail by Arnoup, 
PLANK and SCHOEN BORN [2], and may be described 


Below 
through the 


briefly as follows. point A the liquid 


drains freely plates —the liquid 


‘rains.”’ Above point 4, bubbling of gas through 
the liquid occurs at some of the 
At other perforations the liquid “ weeps *’ 
through the Point B is called the weep 
At point B, all the holes are bubbling, 


and above this point has been considered to be 


perforations. 


he les. 


potnt, 


the normal operating region of sie e-tray columns. 
The curve shown as Type II generally occurs 
only at low liquid rates (less than 350 lb. liquid ‘ 


hr ft* of column) and then only with certain plates, 


In the work reported here only Insert IV, operat- 
ing at liquid mass flow rate of 350 1b hr ft®, 
Type Il. Arnowp et al. [2] 
different plate (with 
0-376-in. holes, 11-5 per cent free space, 0-029-in, 


showed a curve of 


reported data from a 


plate thickness, 270 lb/hr ft® liquid flow rate) 
Below point D, 


liquid is draining through the 


which operated in this manner. 
Curve II, Fig. 
perforations in the plate and no bubbling occurs. 
At point D, which might be called a weep point, 
the gas flow begins to support liquid on the tray 
and bubbling begins in some of the perforations, 


SSIBLE 
| FLOCONG | / 
| (SEE TEXT) 
ad 
| Cc 
\8 
om | | VOL 
9 
is 


The 


sometimes bubble. At still 


other perforations sometimes weep and 
higher gas rates, 
weeping appears visually to stop. However, 
there is no break in the performance curve. 
One other observation over the range of gas 
flow rates noted here, and by AkNoup et al. [2], was 
that oscillation never occurred with a system 
having a Type II performance curve. 

The oscillation behaviour in perforated-plate 
[21], 
and more recently by McAuustrer and PLANK 
[17]. The latter showed that the effect of the 


oscillation increased the plate elliciency, since 


columns was described first by PLANnk 


it decreased back-mixing on the tray. However, 
indications are that a detrimental effect may 
result due to possible flooding in the oscillating 
region. This latter effect is discussed more fully 
below. When oscillation occurs the whok liqquid 
and frothing mass on the tray moves back and 
forth perpendicularly to the direction of flow of 
liquid. That is, if the net liquid flow is from left 
to right, the surging is established from front to 
back. 

For plates and liquid flow rates having per- 
formance curves of Type III, oscillation either 
began at point B, or at a point such as C. In 
this work in general for the thick plates, over yin. 
thick, oscillation began at point B, while for the 
thin plates, 0-029in. thick, 
point C, For the latter plates it is significant to 


oscillation began at 


emphasize that no additional break in the per- 
formance curve occurred at the oscillation point. 
In all cases where oscillation occurred there was 
periodic “ dumping” of the liquid through the 
plates on the same side of the plate as the crest 
of the wave. The holes on this side of the plate 
bubble as the 


reached its peak, and then liquid discharged 


momentarily ceased to wave 
through the holes as the wave began to fall. 

“ Flooding ” in a perforated plate column can 
result from several causes and hence is shown as 
a zone of broken lines in Fig. 1. Flooding can 
result from (1) backing up of liquid in the down- 
comer, (2) expansion of the froth into the tray 
above, (3) splashing of the liquid into the tray 
above due to oscillations. (4) complete lifting of 
the liquid from the tray floor, (5) carry-over of 


all the liquid on the tray as entrainment and, 
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(6) perhaps other possibilities, Types 1 and 2 
are quite common, and result in a vertical AP 
line at the flooding velocity. Depending on the 


tray spacing and physical properties of the 
system, especially the frothing characteristics, 
flooding of this type might occur at any point 
between A and E. 


Once the oscillation point is 


Type 3 is usually only a 
partial flooding. 
reached, further increase in the vas rate at a 
given liquid rate only increases the violence and 
amplitude of the oscillations. In this study the 
splashing was quite often vigorous enough to 
carry liquid up into the tray above, even though 

to of 
PLANK [21] 
noted action of Type 4, where all the liquid was 


an “average” froth height of only 


the tray spacing was observed. 
lifted off the tray and hung as suspended drops 
hetween the plates accompanied by heavy entrain- 
ment. Whether or not this tray action can be 
considered to be a type of flooding is largely a 
matter of definition. This type of flooding is 
quite unstable and very easily changes to flooding 
of Type 5. Still further increase in the gas 
velocity would inevitably result in complet 
carry-over of the liquid in the form of entrain 
ment (Type 5) possibly giving a curve such as 
F in Fig. 1. 


than one type of flooding simultaneously, or in 


There may be combinations of mors 


all likelihood other mechanisms which result in 
flooding. In the work reported here flooding of 


Type 3 only was observed. 


EQUIPMENT 


The rforated-plate column d here WAS the 
Many of the 


details of construction, including some photo 


same one used by ARNOLD ef al. [2]. 


graphs, are given in the above reference. A 
summary of some of the pertinent dimensions of 


the column is given in Table 1, 


Table 1. Perforated-plate column dimensions 
Number of plates 3 
Inside column diameter 15 in 
Tray spacing 18 in, 
Inlet weir height 4 in. 
Outlet weir height 2 in. 
Height of splash baffle (outlet side) 8 in, 
Clearance between outlet weir and splash baffle 2 in. 
Clearance under baffle, from floor of tray } in, 
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Table 2. Plate insert dimensions 


Plate Hole diam Pitch Free space,* 
Plate thickness Number (FA) 
insert T. (in.) of holes per cent 


1-13/64 77 

19/32 

27/64 

14/64 5-33: ziso 

Viv moo OST 160 


* Based on column area 


tinserts 1-4 and 1-5, respectively, used by Annoup. Prank and ScHoennonn 
| 


The perforated plates in this column were plates, TD, ~ 2-3, the magnitude of the friction 
removable inserts for e: nging plate through the holes themselves becomes important, 
thickness (7), hole diameter (D,). and distance requiring a term for the pressure loss due to 
between holes (p). The dimensions of the inserts friction of the Fanning type. 
are given in Tabl 


All the holes were drilled or punched on a tT 
9-0 


‘OF 


triangular pattern. It is important to note that 
there were no holes closer than about Ig in, 
downstream of the inlet weir, or closer than 
Ig in. upstream of the outlet weir. It was noted 
in this investigation, and in others [2, 14-16], 
that if the perforations are too close to cither the 
inlet or the outlet weir these holes will weep 
continually. Also, if the holes are too close to the 
outlet weir, excessive froth passes under the splash 
baflle and over the outlet weir. Inserts I-V 
were made of Plexiglas, whereas Inserts VI and 
VII were made of stainless steel. For the Plexivlas 
plates, holes were drilled and the burrs wer 
carefully removed. The stainless steel plates were 


punched. 


RESULTS 


Dry plate pressure drop. Th dry pl ite pressure ib/hr- tt" 
drop as a function of the mass flow rate of air is Fic. 2. Dry-plate AP as a function of G. 
shown in Fig. 2 for Inserts The pressur 
loss as gas passes through a dry pertorat dl plate 
results from (1) a contraction loss, (2) a friction Following the approach of Hun et al. sudden 
loss in the hole, and (3) an expansion loss. contraction losses in pipes [22] can be expressed 
Hunt, Hanson and Witke [5] at a 7/D, ratio as: 
of 1-0 correlated their results. considk ring only 


the contraction and « xpansion losses. For thicker 


1-25 (1) 


: 
52 
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For expansion, Figure 4 shows a comparison of pressure drops 


4, \* 
(2) 
A,! 2g and experimentally determined pressure drops. 


calculated by the use of equation 4 and Figure 3, 


The average deviation of the points from the 


The corresponding equation for friction losses is:- 


line of equality in Fig. 4 is small. 


2g d. 


Combining these three losses for a perforated 


plate, and substituting plate thickness for 


in equation (3) results in the following relation 


for a single plate : 


(1-25 f(T D,) | 


| 


The factor k has been added in the formulation of 


equation (4) since equations ( 1), (2) and (3) were S ¥ 

derived for a singk hole (a pipe) and not an : 4 yY. 

array of parallel paths. Presumably coetlicient k 406 DO 

can take care of interference between adjacent A ATED OR ATE PF RE OF n. H2C 

jets, and the Couette correction and velocity- Fic. 4. Comparison of calculated and reported dry-plate 
gradient irregularities within the holes. For each pressure drop. 


insert used here, h was plotted against the quantity 
[O-4 (1-25 1, A, } t/ ( T D,) (1 A, A. *) Wet plate pressure drop. Figs. 5-11 show the 


1 2 g. In each case a straight line resulted performance curves ol all the inserts used here. 


having a slope of k. Likewise data from all other The wet-plate pressure drop is given for the 


re adily available sources were treated in the same indicated liquid rates as a function of the mass 
manner. Fig. 3 shows the calculated & values as —— — a = 

a function of 7, D,. The agreement of the points T= V6 in | 
with the line is remarkable in view of the range 0, = 3Ne ir 

of 7 D,. column diameter, free space, and hole fie JF | 


velocity represented in Fig. 3. 


Fic. 3. Dry-plate coefficient k as a function of T/D,. Fic. 5. AP as a function of G — Plate insert 1. 
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ei = 1600 bb/tr 


= 3000 


ib/hr -ft* 


Fic. 6. AP as a function of G — Plate insert I. 


flow-rate of air, G, in lb air/hrft® of column 
cross-sectional area. All of the curves are of the 
Type LLL shown in Fig. 1 with the « xception that 
in Fig. 7 for a liquid mass flow-rate of 350 Ib 
water/hr-ft". This curve is of Type Il. The 
oscillation points are indicated in Fig. 4-10 as 


the small arrows. 


G, 


AP as a function of G 


Plate insert IIT. 


2 


= 350 —ft 
© =1600 - ft? 


| = 2000 Ib/hr -ft? 


; 


= 
oc 1000 
G, 


. 8 AP asa function of G — Plate insert IV. 


In Fig. 12 and 18 total pressure drop minus 
the clear liquid head, as a function of the mass 
flow-rate for Inserts V and VI respectively, 
has been plotted. The solid line shown in these 
Figures is the dry-plate pressure drop. The 
vertical distance from the line to-the wet-plate 


data points is known as the “ residual ” pressure 


[5]. Fig. 14 shows the residual pressure drop 


as a function of 7'/D, at a single mass flow-rate 


= 350 ib/hr -n? 
= 600 ib/hr- ft? 
= 3000 ib/hr-ft? 


1000 
G, 


AP as a function of G — Plate insert V. 


T=V2i 
| | | | 
| 
4 || | | 
| | | | 
| WA 
0-8} = 350 lb/hr 0-8 vA = 
50K 2000 
500 100 1000 20K 
| 
V4 in | T= V2 in 
3°0 } | Op 3/32 in | 
6 | 
oy” 
e £=600 lb/hr | 
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30 


L£=350 
Ib/hr-te? 
4 =1600 |_| 
£=3000 Ib/hr-ft? 


+ + 


6, 


Fic. 10. AP as a function of G — Plate insert VI. 


of air, At other values of G, the residual pressure 
drop increases with increasing G. The data in 
Figs. 12, 13 and 14 show the residual pressure 
drop to be a function of 7'/D,, L and G. An increase 
of L, for thick plates, results in an increase of 
R. Likewise R increases as TD, increases for 
constant values of L. No general correlation 
was made, but these data were given to indicate 
the general trends. 

The results of this investigation and 
obtained by Arnon et al. [2] show that for thin 


those 


10-0; 
9-0 
8-0) 
7-0 
6-0) 
5-0} 


e L=720 ib/nr-ft* 
=1660 Ib/hr-ft? 
= 3200 Ib/hr-ft? 


"306 B00 1000 
G, 


2-01 


AP as a function of G — Plate insert VII. 


Fig, 11. 


Perforated-plate performance 


© = 3000 Ib/hr -ft | 


0-6} in 


0» = 3/32 in 


400 600 800 1000 


G, ib/hr-ft 


Total pressure drop minus the clear liquid head 


as a function of G — Plate insert V. 


plates the residual pressure drop is affected little 
by changes in liquid flow-rate. The ratio of the 
residual pressure drop to the dry-plate pressure 
drop remains about constant for the thin plates. 
Thus the data lie above and in a line parallel 
to the dry-plate curve in Fig. 13. Fig. 24 of 


ref, [2] shows the same trend for the thin plates 


9-OF 7=0-029 in|} 
7-O}-4=0°250 ini 


| 

4+—_+—+—4 + + 


5-0) 


0-2 
200 


300 400 800 1000 2000 


° 
6, Ib/hr —ft* 


Fig. 13. Total pressure drop minus clear liquid head as 


a function of G — Plate insert VI. 
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Table 3. Weep points 


| 


Plate | L G 
insert | (Ib H,O/ her ft*)| (in. H,O)) (Ib air/hr ft*) 
No 


B50 610 
1600 555 


3000 470 


350 2 770 

1600 710 

2-37 620 
Fic. 14. Residual pressure drop as a function of T D, 


and L, 350 7 


used by ARNOLD et all. [2]. For the thicker 910 


plates, however, the residual pressure drop 

350 
increases with an increase of liquid flow-rate. 


The ratio of residual pressure drop to dry-plate eeee 


vas 


pressure drop decreases with increasing 
tlow-rate. 350 
The complete set of data for Inserts I-VI ron) 


giving total pressure drop, clear liquid height, 


froth height, gas and liquuid flow-rate, ete. has ano 
been tabulated by McGinnis |. 1600 


Weepage correlation, Fig. 15 shows the total “HMM 


pressure drop plotted as a function of the vapour 


F-factor through the holes at the weep pot. 


The Vapour F’-factor is the square root of kineti 

1660 
2020 
240 


energy of the gas per unit volume. Table 3 
presents the weep-point data of this investigation. 
The solid line in Fig. 15 is the weepage correlation 
proposed by Hucumarx and [1]. 
The data of the present investigation agree 
reasonably well with the correlation. However, 
some of the data points shown in Fig. 15 are 
not in the range of hole diameters included in 


HuGumark and correlation: thus 


the weep-points for the hin. holes are out of the 


range of the correlation and do not agree with it. 
Also the data for the 3,32 in. holes (T' dD, 5°33) 


WET-PLATE PRESSURE DROP, in. H2C do not agree with the correlation. It is felt that 


15. Weep-point correlation there is a surface tension effect which should 


| 
| 
| 
G6 = 120% b/hr -ft* 
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Table 4. Oscillation points 
Plate | Z.. 
insert L | G AP (in. 


No, |(Ib H,O/hr ft") (ib air /hr ft*) (in. HO) 


350 1220 2-22 0-68 
1600 | 1160 | 242 | O91 
3000 1150 260 | 1-00 


1330 
1260 
1100 


350 
1600 
S000 


350 1470 
1600 1280 2-43 
S000 1220 


1375 


12H) 


1600 
3000 


100 
1600 1050 2-76 


1030 


VI | B50 1238 6-81 
j 


1600 
S000 


Vil 12233 709 0-57 
120 OSS 
1211 7-87 ool 
1380 7-83 O-ST 
1660 111 7:13 O71 
2020 1182 1-14 
2450 1170 
2610 1166 8-03 
3200 1148 751 130 


be included in the correlation, and that this 
dD, ratio, 


presents 


effect may be magnified by the high 7 
Table 


the oscillation point data of this investigation, and 


Oscillation- point correlation. 


Fig. 16 presents the air mass flow-rate as a function 
of the clear liquid height, both at the oscillation 
point. The agreement is quite satisfactory 
except for Insert V (j in. thick, 3 32 in. holes). 
It should be noticed that the weep-point correla- 
tion for this particular insert was also unsatis- 
some surface tension effect 
both 


improve the agreement of the data for Insert V 


factory. Possibly 


incorporated into correlations would 


It is interesting 


with that of the other inserts. 
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to note that the oscillation-point correlation was 
based on the column mass flow-rate whereas the 
weepage correlation required the hole velocity 
to bring the data into the proper order. 


TT TT) 


pi 


Fic. 16, Oscillation-point correlation. 


The clear liquid height, Z,, at the centre of 


the operating tray was used in the oscillation- 


correlation, and was measured with a 
One leg of the 


inserted through the plate from beneath and was 


point 


manometer. manometer was 


made flush with the top surface of the plate. 
The other leg of the manometer was mounted in 
the vapour space above the plate. The mano- 
meter reading itself was considered to be the 
clear liquid height, measured in inches or centi- 
metres of liquid. It represents the hydraulic 
pressure of the liquid and froth on the plate in 
terms of clear liquid head. The clear liquid height, 
along with the froth height has been shown by 
quantity in 


PLaNk [24] to be an important 


characterizing aeration on the tray. Aeration 
in turn is one of the factors which affects mass- 
transfer elliciency in liquid phase controlling 
systems [24]. 

The clear liquid height is a difficult quantity 
to measure accurately. It fluctuates constantly, 


especially in the oscillation region described 
above, and is often less than 1 in. in magnitude, 
even with a 2-in. weir. The clear liquid height 


and the froth height undoubtedly hold many of 


| 
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Il 233 0-60 
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the keys to a complete understanding of tray 
hydraulics, as well as mass-transfer efliciency on 
perforated-plate trays, bubble-cap trays, and 
the like. 


CONCLUSIONS 


(1) The dry-plate pressure drop for perforated 
plates is described well by equation (4) and Fig. 3. 

(2) At constant hole diameter and free area, 
dry-plate pressure drop decreases with increasing 
plate thickness (i.e. 7'/D, increases) until hole 
friction becomes appreciable (at 7D, ~ 2-3) ; 
then a reversal of this trend occurs. 

(3) Four distinct types of plate action were 
encountered in this study -—raining, weeping, 
stable operation, and oscillation. 

(4) Periodic liquid and froth surges or oscilla- 


tions can occur at high gas rates and is accom- 


panied by periodic liquid dumping through the 


perforations near the column walls. 

(5) For thin plates, residual pressure drop is 
shown not to be a function of the liquid tlow- 
rate, it increases with increasing gas flow-rate, 
but the ratio of residual pressure drop to dry plate 
For thick 
plates the residual pressure drop is a function 
of D,, G and L. 

(6) The 
agreement 
Hucumark and O'Connet 

(7) The 
column area, at which oscillation begins is shown 


pressure drop iS ESS¢ ntially constant. 


results of the investigation are in 


correlation of 


with the weepage 


mass flow-rate, G, Ibair/hrft® of 


to be a linear function of the clear liquid height 
at the oscillation point. 


NOTATION 
cross-sectional area of the larger pipe, ft® 
cross-sectional area of the smaller pipe, ft*. 
column cross-sectional area, ft*. 
total perforation area, ft*. 
diameter of the perforations, in. 
diameter of the smaller pipe, ft. 
pressure drop across a plate, in. of fluid. 
Fanning friction factor for smooth pipes, ref. [20] 
total perforation area-to-column area ratio, per 
cent. 
F-factor through the 
gas density, Ib/ft*. 


holes, 1), V pe: where 


p 
m of gravity, ft/sec. 

mass flow-rate of gas, lb gas/hrft® of total 
column cross-sectional area. 

head loss across a dry perforated plate, ft of 
gas flowing. 

head or pressure loss due to contraction, ft of 
fluid. 

head or pressure loss due to expansion, ft of 
fluid. 

head or pressure loss due to friction, ft of fluid 
a dimensionless coeflicient. 

length of pipe, ft. 

mass flow rate of liquid, Ib liquid /hr ft? of total 
column cross-sectional area. 

pitch, the centre-to-centre distance between the 
perforations, in. 

residual pressure loss (total pressure loss minus 
dry plate loss minus Z,). 

plate thickness, in. 

average velocity in the smaller pipe, ft /sec 
average velocity in a perforation, ft/sec. 

clear liquid height, in. 
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Abstract 
function of the drying tin 


In batchwise drying the 
is well 
less parameters, and based on some 
this method of drying which pern 
material and of the drying potential of the 


moisture content of the material that 


us of the location of the material in the drier 


ar 


is to be dried is a 


Using dimension- 


simplifving assumptions, a method of calculation is given for 


its prediction of the distribution of the moisture content of the 


This is possible regardless of the shape of the 


curve which relates rate of drying to moisture content, as obtained from laboratory drying tests 


under constant drying conditions 


Résumé 
la durée de séchage de méme que ci 
puranie tres sans dimensions, et puart 
calculation est présentée que permit de 
sechage de | 


au teneur « 


dle la force cde air 


lie la vitesse ce sechae 


des conditions de séchage invariables 


im 


Funktion der 


Zusammenfassung 


yvutes eine Trocknunyszeit 


dung dimensionsloser Parameter und mit 


Trocknungsart eine Berechnunysmet hoc 


und die trocknende Wirkung der Laft \ 


ohne Riicksicht auf die Kurvenfor 


dem Feuchtivkeitsychalt, wie sie aus 


bedingungen erhalten wird. 


IN convection batch dri rs drying is achieved by 


forced convection of air 
ol 


tray and compartment dryers as 


or th nich 


stationary arrangement wet solid material. 
Examples are 
used in the heavy brick and clay industry, 
truck 


driers for drying thick beds of various 


tray 


and driers for drying wood, sieve tray 


vranular 
materials. 
with this method of drving two 


both relating to the 


In connexion 


main problems arise, uni- 
formity of the moisture content of the products, 
One of them is, how to obtain a uniform distribu 
This 
though 


sometimes dillicult, a satisfac tory solution can be 


tion of the air across or through the stock. 
problem is an aerodynamical one and 
found in most cases, ¢.g., by flow experiments on 
small-scale models to determine the appropriate 


dimensions and arrangements. 


la position du matériel dans le 


clire 


und des Ortes innerhall des 
mitveteilt 


Laboratoriumsversuchen 


Dans le séchage par charge le tencur en eau du matériel A sécher est une fonction de 


s«hoir, En se servant de 


int de quelques suppositions simplifiantes, une méthode de 


la distribution du teneuren eau du matériel et 


Cette méthode est prosst lle suns compter la forme de la curve qui 


eau, résultant des expériments de séchage laboratoires dans 


nweisen Trocknen ist der itsye halt des Trocke n 


lrockners, 


\nnahmen 


nter Verwen 
vereinfachenden 

die Verteilung des Feuchtigkeitsgehalts 
yestattet. Dies 


wird fiir diese 


Voraussage ist mévlich 


ler Beziehung zwischen der Trocknungsgeschwindikgeit und 


unter konstanten Trocknunys 


However even if a uniform distribution of the 
be realized there the 
of the content of the 


that is inherent in this method of drying. It is 


formuity mosture materials 
caused by the increasing humidity of the air on 
its way through the drier and by the fact that the 
rate of drying ce yn nds on the moisture content 
the 


material is a function of drying time and of loca 


of the material. Since morsture content of the 


tion with reference to the entrance of the air. the 
situation is more complicated than in the case of 
continuous driers; so far we know of no general 
calculation methods, as are available for contin 
uous driers. This is perhaps one of the reasons 
why the use of batch driers is gen rally avoided 
if possible, 

which 


However for materials 


require a long 


drying time, such as green brick and wood, batch 
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driers must be used. The problem arises of pre- 
dicting the drying conditions and calculating the 
drying tin required to kee p the moisture content 
of the product between given and generally 
rather close limits. Also certain combinations of 
moisture content, humidity and temperature 
conditions of the air which would cause deteriora- 
tion in the material must be avoided. 

This paper gives a method for calculating the 
moisture content of the material and the humidity 
of the air as functions of drying time and location, 
some idealizing assumptions being made. 

Rate of drving v. moisture content curves 
determined under constant conditions by labora- 
tory drying experiments are assumed to b avail- 
able. It is assumed that the adiabatic saturation 
temperature of the air remains constant. In case 
of recirculation, the time required to return air 
from the exit to the entrance is neglected. 
Furthermore the assumption is made that a one 
dimensional treatment possible. many 
driers of the batch-type the actual configuration 
can be arranged so that the same conditions for the 
material and for the air can be expected in cross 
sections normal to the direction of the air flow, 

Another assumption is discussed in the deriva 
tion of the differential water balance, i.e. the rate 
of change of the water vapour content of the 
air in the drier is neglected with respect to the 
net rate of water-vapour transport by forced 
convection. 

Two cases of initial moisture distribution that 
are believed to cover most practical cases are 
treated in this paper. Mathematical details not 
directly needed for understanding the calculations 
are omitted to reduce the size of the paper. We 
intend to check the validity of the calculations 
by experiments, but present suflicient data are 


not yet available for insertion in this paper. 


Tue Eevarion vor tHe Rare or 

The rate of drying of a wet solid in flowing air 
depends on drying conditions in a complicated 
way. For most engineering purposes, however, 
it can be described with reasonable accuracy by 


the equation : 


ol 


k.V.g(F) (1) 


Adiabatic convection batch drying with recirculation of air 


As a first approximation the various factors 
in this equation are separated into three in- 
dependent groups : the mass transfer coefficient 
k which depends on the air flow conditions, the 
driving force or drying potential V which is a 
function of temperature and humidity of the air 
and finally the dimensionless function g(F) of 
the mean water content which takes into account 
the physical properties of the material to be dried. 
In most cases the functions k and g(F) must be 
obtained from laboratory drying experiments 
under conditions approaching as closely as 
possible the actual conditions in a drier. 

If. as is often the case, the critical moisture 
content F. and the equilibrium moisture content 
| can be assumed to be constants it 1s more con- 
venient to write g(F) f(F F./F, 
this function f having the value 1 for F KF 
and 0 for F 

kor low temperature and humucity V equals 


H). The formula 


= (2) 


P 


| 0-622 In 


has a more extended range of applicability. 
By converting partial pressures into specific 


humidities : 


0-662 
0-662 In (3) 
fi 


If a correction factor 6 is introduced : 


If adiabatic saturation of the air is assumed H, 
is a constant. For the temperature ranges 
occurring in practice 6 will only vary slightly. It 
appears therefor that in cases of high temperature 
and humidity drying can be approximately 
described by apply ing a constant correction factor 
8 to the drying potential (H, H). Hence 


equation (1) can be written as: 


k.B.(H, ms 4 (5) 


Tue DirFERENTIAL WaTeER BALANCE 


In this paper tray and compartment driers 


are treated which, for purposes of calculation, can 
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the 
Of course 


be represented, at least sectionally, by 
mathematical model shown in Fig. 1. 
the possibility of doing this requires an equal 
distribution of the air through the actual drier. 
For the model of Fig. 1 a one-dimensional treat- 


ment is assumed to be applicable. 


Materia! to 
be dred 

Fic. 1. Mathematical model of drier 


Considering an element of volume between 


cross-sections at z and 2 Az during a time inter- 
val Af, a water balance can be written by « qualizing 
the amount of water vaporized from the wet solid 
to the sum of the net amount of water vapour 
transported by convection and the increase of the 
water vapour content of the air in this element. 


If a factor b. Az. At is dropped, the 


equation is obtained : 
h P Hi 
The last term of this equation can b approxi 
(hy,/Q) 5H /dt if y, is the 
the 
temperature. Now in the equation : 


following 
(6) 


mated by specific 


weight of dry air at adiabatic saturation 

hy, 1 
= A 

Q od oa 


the ratio of the coeflicients of the second and third 


(fia) 


term is approximately equal to the linear air 


velocity on a dry air basis. Under conditions 
generally occurring in practice the rate of change 
of the water-vapour content of a given element 
of the drier will be small in comparison with the 
net rate of vapour transport from the element by 
forced convection. For most practical cases it can 
therefore be expected that a good approximation 
is still obtained if the last term of equation (Ga) 
is neglected and the calculations are based on the 


approximate water-balance : 


16F 
Q dt dbz 
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Dimensionless equations 


Equations (5) and (7) can be made dimension- 


less by introducing the parameters : 


het 


8 kz 
AQ? 


Equation (5) transforms into 
(dh) 
and equation (7) into 


(10) 


The parameter 7 can be eliminated from (0) 
and (10) giving : 
"4d df 


» > 
‘ } ‘ ‘ 


(11) 


With an arbitrary function P (0) inter- 


mediate integral of (11) is : 


Another equation that is needed gives the 
drying potential of the air entering the drier as a 
function of the drying potential of the fresh air 
and the recirculation ratio. Again it is assumed 
that the adiabatic saturation temperature is a 
constant. Furthermore the time needed to trans- 
port the air leaving the drier back to the entrance, 
where it is mixed with fresh air. is neglected : 
(1 (13) 


+ aly 


Here Vi, and respectively represent the 
values of V at the entrance of the drier. of the 
fresh air and of the air at the exit of the drier. 
CALCULATION OF THE DrytIne Process 
Initial moisture content above the critical moisture 


content 


Now the case is considered when drying starts 
at all 
points in the drier, Three periods can be dis- 


with the moisture content greater than F. 


tinguished : 


om 4 
d 
(9) 
v7 vi VOL 
9 
19528/ 
: 
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First period : Everywhere F > F. Drying rate 
independent of F. A moisture gradient develops 
through the drier. This period terminates when 
F = F, at the entrance. 

Second period : Moisture contents both higher 
and lower than F, are present. The critical 
moisture content moves through the drier. 


Third period: F < F, everywhere in the drier. 


First period 

In general the initial moisture distribution will 
be uniform: F = F,. Sometimes, however, 
after the period of heating up a non-uniform 
distribution is caused by condensation of water- 
vapour on the wet material at the moment at 
which drying really starts. In order to include 
such a case a more general initial distribution ts 


assumed : 
(Cf, 0) dy, d, (14) 


in which 4, can be positive, negative or zero as 
long as 0) for A. During this 
first period we can take f (¢) 1. Then equation 
(12) takes the form : 


From equations (9), (10) and (13) : 


7) (16) 
l 


and with equation (14) : 


7) dy, 
l ac 


Of course 7 appears to be independent of + 
and for a given ©, @ decreases linearly with +r. 
It is further noted that a time 7 can always be 
found for which 46 — ¢d,. With (16) and (17) the 
first period is fully described. 

The second period starts at + r, for which 


(0, 1. It follows that : 


T (dh 1) (18) 


and the distribution of the moisture content at 


this moment is given by : 


T,) dy (dy, (19) 


Adiabatic convection batch drying with recirculation of air 


This distribution is independent of ¢,, a, 7% 


and A. 


Second period 
From the initial distribution for this period and 


from equation (12) it follows that : 


1 4 (24 
- - @ 20 
Fd vt + 0 ) 
in which again f(¢) 1 for ¢ = 1, whereas for 


0 A 1, the function {(¢) has an arbitrary 
shape which is known from drying experiments. 
By integration of (20) : 
ds 7, 
M (r) (21 
(dy 8) (8) ) 
If the critical moisture content ¢ 1 is 
situated at (., it is easily derived from (21): 
M (+r) 3 


ds 
(dy s) f(s) 


(22) 


It is noted that 4 depends on & 


distribution for € > ¢ is given by : 


Cord = do — (% — I) 22a) 


The value of 7 is calculated from equations (9) 
and (22) to be : 
dl 


= 23 
(23) 


(dy 
in which dC, dr Is the dimensionless velocity of 
the critical moisture content. This velocity ts 
calculated from the recirculation equation (13) 
by substitution of 7 from equation (23) for ¢ = 0 


and ¢ as 
(24) 
dr dy (0, a Dy, (A, T) 


The values of 46 (0, and 4(A, r) are found as 
functions of [ from equation (22). 
Equation (24) can be substituted in (23) to 


give : 
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, (28a) 


+) OtA. tr) 


The value of ( as a function of 7 is easily eal- 


culated afterwards by graphical integration of 


dr, with respect to 


Third pe riod 
If for + Te 
third period begins. It is clear, however, that the 


has reached the value A. the 


calculations for this period follow exactly the 
If a 
physical meaning is no longer attached to ¢ 
there 


separate one. 


same lines as those for the second x riod. 


is no necessity to treat this period as a 


Some remarks can be made that consice rably 
reduce the amount of work for calculation of the 
follows that 
and the 
parallel shift of a 


distributions. From equation (22) it 
the distributions ¢ during the second 
third period follow from the 
single curve lving between 4 0 and 4 


and approaching these limits asymptotically 


for © and 0. 
For the second yo riod the distributions for 
od 1 are found by parallel shift of the distribu 
tion viven by equation (19). The easiest way to 
complete these distributions for the second 
period is by calculation of the distribution for 
A from equation (22) and shifting this until 
it fits into the given value of ©. 
The distribution of 7 is found from equation 
23a) or, according to equation (24), by deriving 
dr from the 


distributions and multiplication with 4, 5 as 


the value of dé graph of moisture 


obtained from the same graph. 


If for values of 4 to zero the graphical 


integration of equation (22) becomes diflicult on 


account if the great values of the inteerand. the 
function f(d) can be approximated in this range 
by a linear function permitting analytical integra 
tion. 

Linear falling rate section : 

In the falling rate section we can put for many 


materials as a good approximation / (¢ »b. In 


this case the graphical integration can be re placed 


by an analytical one and some formulae can be 


derived. 


VAN 


Of course for the first period the formuale given 
before remain unchanged. 


For the second ype riod integration gives : 


To determine (Ur) from equation (23) 


dt. must be calculated from equation (13): 


dr 


The following observation can be useful. For 
some materials it is nece ssary, in order to avoid 
deleterious effects, to keep the driving force V 
value as long as the moisture 


below a certain 


content is above FF. If this maximum allowabk 
value is applied at the beginning of the second 
period it can be shown that, for high values of the 
that 


some time the drying pot ntial at © rises above 


recirculation ratio, there is danger after 
this value. If no recirculation is applied, the value 
of at is given by 

(28) 


risk 


function of 


monotonic 
the 


and there is no since this is a 


decreasing which tends to 


limit 


In fact this allows a further raising of 7, after 


third period formulae (26) and (27 


some time with nearly a factor 4, /¢ 
For the 


become 


o 


(2%) 


(la), 
ae“ (26) 
Integration gives as a function of : 
> (dy, Ll) 
(1 (7 T,) In VOI 
7 
~ 27 
a(dy — (27) 1952/ 
Vo 
and 
3 
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calculated before as a function of 4,, equation (35) 


(l — Tt.) = In 


can be integrated to give o, aS a function of +r: 


(dy 
aln (30) db, (s) 


f(s) 


A simple formula for the total drying time +r, 


needed to obtain a value Py ° I at > \ for The calculation ol 7T proceeds AS follows : 


drying under constant conditions can be cd rived From equations (9) and (34) : 


from equations (18), (27) and (30) ; 


hence with (35) : 


(1 a) (38) 


alnd In 


which formula is the Same as (23a): 7 can now he 


If d,,. this can be simplitied to: 


calculated from the set of distributions 6. 


(32 ) 
I Linear falling rate curve 


In case of f(s s, the formulas eiven in this 


Homogeneous mnitial distribution below the critical section can he integrated directly Formula (33) 


moisture content IVES 


d.. d 
If the initial distribution is denoted by De, | 1) A PO Pi € 39) 
equation (20) can be used again. A treatment Po Py Dy 
slightly diffe ring from the fore poihng one is given Krom this equation 7 is easily found by means ot 
below. Kquation (20) 1s now integrated to : equation (38), 
kor (s) formula 36) inte crates to: 


with ca +r). a Dy € 

(¢ 1) In 

Taking an arbitrary value of 4, ( d,) the a — dy ¢ 
distribution is easily completed by graphical d a d 

= in — 4 (e*e* 1) In (40) 
integration of equation (33). In this wav a com . vs 
a Pi 

plete set of distribution curves can be obtained. 
lhe drying time can be calculated as follows. CHANGE Dry inc CONDITIONS 


From equations (9) and (33) : 
According to equation (22) the drving potential 


l dd l dd, 34 Vand the recirculation ratio a do not affect the 
(34) 
(dy (d,, ?,) dr shape ol the moisture distribution >. as long as the 
, air velocity Q remains unchanged. Formula (24) 
From (34) and (13), if 4 d(A, r): 
shows that drying time is inversely proportional 
l dy 35) to as 1s to lhe exper ted from equation so 
f(d,) dr dy, d, the effect of changing this parameter can easily 
be taken into account. \s mentioned earlier 
(tr r dy d and —" ii. (0, and r) can be regarded as functions 
of «& In case of a recirculation ratio a b ing a 


Since ¢, is known trom the distribution é function of time. integration of equation (24) 


(do — 1) + 1 
a ds (1 a)m,, (36) 
do, d, ] Id, 
4 — a) — 1) +O — 
dy, 
VOL. 
9 
58/59 d,, A 
4 
a 
(do — 8) f (s) arate” 
t1 
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remains possible in principle, resulting in the 
drying time + corresponding to a given value of ©. 
However, in case of stepwise variation of a, the 
new drying time can more easily be determined 
from repeated graphical integration of equation 
(24). 

In case of varying air velocity Q the situation 
becomes somewhat more diflicult since in general 
the mass transfer coefficient k will vary too. As 
@ and k enter into the equations (8) for + and {, 
both the drying time and the distributions of the 


water content will be changed. However, for 


oh air velocities k 


Q”* 


will vary approximately as 
and therefore the effect of varying Q onk /Q in 


will be rather small. As a first approximation the 


distributions of the water content can be assumed 


to be independent of Q under these conditions. 


EXAMPLI 
that the 


simple 


Suppose actual drier can be 


to the 


reduced 
shay of Fig. 1 for purposes of 
The following data are given: L 

Olm; Fy =—725kg/m*; F, = 3-62 


ko 0-72 ko 10,000 kg 


calculation. 


h 


m*h, 


corresponding with a linear air velocity of about 


200 ke m*h. V 


\ir conditions: fresh air supply at 105°C; 


ad.sat.temp. 80°C ; 


0-063 
0-622 
0-622 


“S511 
0-53 ; 


The dimensionless parameters are : 


F ‘72 0-538 200 


2-9 


10.000 
0-106 2: 
0-53 200 


A 0-106 10 1-06. 


Drving is assumed to start from a homogeneous 
initial distribution F 


7°25 
= 


VAN MEEL 


Recirculation ratio a 5/6. 
The case of drying under constant inlet conditions, 


The 


curve {(¢) is supposed to have the shape given in 


as given above will now be treated. rate 


Fig. 2 as determined by laboratory experiments. 


Fic. 2. Dimensionless rate curve. 


First period 


At the end of this period the distribution of the 
water content follows from equation (19) : 


7) 2°25 


125 e 
Drying time according to equation (18) : 
1-46 


or h. 


(y, 


{Curvel 


a | | 

*2 + + > + + + 
04 O8 16 : 
g 
1958/ 

6 

< 0-063 
6 

|| 

30 ae T 3-0 

| | 

0-5! 2-0 

1-8 

1°6 “ ; 

ae + + + +> + + 6 
0-2 

O2 04 06 O08 TO 
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¢ 


Fic. 4. Dimensionless distribution of water content. 


0-026; From equation (16): O-OLAT Ge or 


0-024) | | | ] IV O-OOTR ¢ If no recirculation is applied thre 


0-020) drying time will be 7, IDStort, — O-54h. 


0-016) 
0-012 
0-008) 
0:004} 


| 
| 


Second and third period 
From f (¢) as given by Fig. 2 the function 
l 
(do — f (¢) 


is plotted as curve 1 in Fig. 3. By integration the 
function 


43 


2°46 
2-0 
VOL, 
9 
52/59 
7 x 
"6 24 3+2 5.05 T=257 
Fic, 5. 
T=225 
7=210 
02 O68 O8 TOA 
Curve I 
Fic. 6. 
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ds 
s) f(s) 
is derived from this curve as shown by curve 2 
This 
distributions 4, if 4 
At the end of the second pe riod 


equation (22) it follows that curve 2 (Fig. 


in Fig. 3. curve forms the basis of the 

From 
3) gives 
the desired value of 4 at a given value of © In 
Fig. 4 the distribution for ¢ 
are plotted first. By parallel shift of the first one 


to the right and of the second one to the left, the 


Oand for ¢ 1-06 


intermediate distributions for ¢ Oz: Ob: 


found. From curve 2 (Fig. 3) distributions 
for the third 


1-2; 1°35; 


O-sS are 


found by inserting 


riod are 
These 


partially be found by parallel shift 


distributions can 
Ones Fig, 


” and 3. 


) 


being completed the values of ¢ 


of 4 (A, +) can be obtained from this graph. 


This enables us to calculate dC. dras a function of 


C. as shown in Fig. 5 According to equation (23) 


the distribution of the drying potential ts found 


by multipheation of 4 é@ as taken from Fig. 4. 


with the corresponding value of di it read 


from Fig. 5. The result is given in Fig. 7. \s 


before, the recirculation causes a 


The 


mentioned 


steady increase of a at ¢. after some time, 


calculations are completed by integration of 
dr dt (curve 1, Fig. 6), 


viving the time r r, belonging to a given value 


of 


Finan REMARK 
Starting with the equations (9) to (12) 
treated, For a 


linear 


general cases could be rate of 


drying curve possessing a falling rate 


VAN MeEEL 


section, calculations have been made by the 
author for the general case of an arbitrary initial 
distribution. In general. however, there will be 
no need for this, and suflicient information will 
be given by a study of the two cases described 
above. 

It will be that, by 


obvious modifications, cases such as the drying 


observed means of some 


of thick beds of granular material by means of 


through circulation can be treated along the 


same lines, 


No. ATION 


part of the mass rate of 


that is 


recirculation ratio 
flow through the 


dimensionless 


drier recirculated; 
correction factor, dimensionless 

moisture content, kg /m® of exposed surface 
initial value of F 

critical moisture content 

equilibrium moisture content 

height of drier, m 
specific humidity, kg H,0 kg of dry air, 
specific humidity at saturation 

mass transfer coeflicient, kg /m*h, | 
length of drier, m 

total pressure, kg /m* 

partial H,O-pressure in bulk of air, kg/m? 
partial H,O-pressure at wet surface, kg /m* 
drvy-air mass-velocity, kg 

absolute temperature, A 

drying potential, dimensionless 

drying potential of fresh air 

drying potential of recirculated air 
co-ordinate, m 

correction factor, dimensionless 


specific weight of dry air, kg/m ,, 


dimensionless parameters 


defined by equation (8) 


value of Cat 2 


: 
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Diffusion in binary mixtures 


C. O. BENNETT? 


Purdue University, West Lafayette, Indiana, U.S.A. 


(Received 6 May 1957; in revised form 1 November 1957) 


Abstract—The effect on the ordinary coefficient of diffusion of the mobilities of two interdiffusing, 


non-ionic components is evaluated by a straightforward application of the equations of irreversible 


thermodynamics along the lines suggested by pe Groor [1 The form of the resultant equations 


for the diffusion coefficient indicates that this approach does not agree with experiment so well 


as that of considering the effect of the different mobilities according to the equations proposed 


by DarkKEN [2] and others, The analysis given in the text indicates that a binary diffusion 


coetlicient is best calculated from the measured rates of diffusion of radio-active tracers by the 


equations of DARKEN [2 


Résumé effet de la mobilité de deux composés non-ioniques, diffusant iproquement, sur 


le coefficient ordinaire de diffusion, est évaluc par une application directe des équations de ther 


modynamique irréversible, en accord avec les suggestions de Dre Groor [7 La forme des équa- 


tions, du coefficient de diffusion, qui en découlent, montre que cette approche du probleme ne 


concorde pas AUSSI bien xperience, que cell qui consiste a considérer action des différentes 


mobilités daprés les équations proposées par Darken [4] ete... L’analyse reproduite dans 


ce texte indique qu-il est préférable de calculer le coeflicient binaire de diffusion a partir des 


mesures de vitesses de diffusion de traceurs radioactifs, d’aprés les équations de DARKEN 4]. 


Zusammenfassung —Die Wirkunyg der Beweglichkeiten von zwei ineinander diffundierenden, 


nicht ionisierten Komponenten auf den gewOhnlichen Diffusionskoeflizienten wird durch direkte 


Anwendung der Gleichungen der irreversiblen Thermodynamik gemiiss den Uberlegungen 


von Dr Groor [7] berechnet. Die Form der Endgleichungen fiir den Diffusionskoeflizienten 


zeigt, dass dieses Ergebnis nicht so gut mit dem Experiment tibereinstimmt als die Betrachtung 


der Wirkung der verschiedenen Beweglichkeiten gemiiss den Vorschligen von Darken [4] 


und anderen. Die im Text mitgeteilte Untersuchung zeigt, dass ein bindrer Diffusionskoeflizient 


am besten aus den vgemessenen Diffusionsgeschwindigkeiten radioaktiver Indikatoren durch 


die Gleichungen von Darken [4] berechnet werden kann 


IN common with many other natural processes, also been used to some extent for measuring 


diffusion can be studied by postulating a diffusion in binary mixtures. It is with particular 


mechanism and attempting to relate the observed regard to this latter situation that the equations 


rate of diffusion and the coeflicient of diffusivity to follow have been developed. Although 
derived therefrom to the basic, known para developed in connexion with the measurements 


meters for the molecules or atoms involved. On of gaseous diffusion at high pressure, because of 


the other hand, the diffusion can be described — their general thermodynamic nature, the relations 


from the macroscopic point of view by the use hold also for liquids and in certain cases for solids. 


of the thermodynamics of irreversible processes. The coetticient of diffusion of ** 1 °° in a mixture 


Whereas the former procedure is probably more of “1” and “2,” D,,, will be detined for a non- 


fruitful for the relatively simple cases for which ionic system by the equation : 


the theory is satisfactory, at present it is necessary 


J D,, (de, dx) (1) 


to be satisfied with the phenomenological descrip 


tion in many non-ideal situations. for flow in one dimension at constant temperature 


Radio-active tracers have proved very useful and density with no driving forces present except 


in measuring self-diffusion coeflicients and have the concentration gradient. J, is the flux in 
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O. 


the molar 


velocity. For the case where p is not independent 


¢ moles ‘em? sec relative to average 
of x (.e. of y,) D,. may conveniently be defined 
by 


J p Dy, (dy, ‘da (2) 


For flow relative to a stationary 


define a tlux N, by 


Je +(N, 


In addition. 


(t}) 


1-0, D., and ire identical: 


diffusion coetlicient 


na Sith’ 


there is only one | the type 


detined by equation (1) needed ior a binary 


The 
but they are presented here again in order to mak« 
The 


factors 


mixture. above results are well known, 


the basis of the following discussion clear. 


need for a careful definition of the basic 
in a diffusion problem has been emphasized at 


[3). 
Note that if the density were a function of 4. 


some length by 
equ (1) would vive 


y (7) 


dp \ dy, 


dy, 


ind many equations in the subsequent develop 
ment would be complicated by the 
of the Therefore it is 


convenient to define by equation (2) for the 


arance 


term involving dp, dy,. 


unusual case of a variable ce nsity. 


The imtroduction of equations (3) and (4) is 


consistent with (6), which is true in 


veneral, By the 


measure self-diffusion in either a pure component 


equation 


use of radio-active tracers to 


or a uniform mixture, experimental conditions 


which N We, 


The measured tlux of the 


are obtained for and so 
J 
thus he simply related to the determination of 


is introduced into the same system as previously 


tracer can 


Note that if there is no tracer, and a gradient 
chemical 


used with the tracer so that ordinary ~* 


diffusion results, N, does not nece ssarily any 


longer equal N,. The actual relation between 


N, and N, 


the diffusion-measuring system and the way in 


will depend upon the geometry of 


which p is related to the composition. Problems 
involving the evaluation of N, and N, have been 
Biro [3]. Jost ([4], page 441) 
applies equations (3) and (4) to solvents evaporat- 


discussed by 


ing through a capillary or to the slow dissolution 
of a salt in a similar geometrical arrangement. 
The between N, and N, 


mined by the change in volume of the salt as it 


relation will be deter- 
dissolves and diffuses out, for example. 

After this preliminary discussion, let us turn 
to the relation between D,, in a binary mixture 
and the measured rates of diffusion of radio-active 
the 
mixture ; that is, a mixture with no concentration 
that the 
shall not here be concerned with the 


tracers through otherwise homogeneous 


gradient due to radio-active 
tracer. We 


mathematics 


pt 
involved in obtaining a diffusion 
coctiicrent from measured concentration gradients 
in an unsteady-state experiment. This topic is 
well covered for many situations in the book by 
Crank [5], for example. 


the 


thermo- 


what understanding of 
the 


For a mixture 


Let us first see 


problem can be obtained from 
dynamics of trreversibk 


({1], pages 100-104) : 


processes. 


of 1 and 2, 


d dy. 
(8) 
ar 


dy, djty 


L + Lies 
dz 


(9) 


21 


Since it argued that there is no 
coupling effect in ordinary binary diffusion [6), 
might set L,, 


For example, doing so and considering an ideal 


might be 


on in equation (8) equal to zero, 


gas leads to the equation ; 


Ly, D,, 


10 
RT 


For an ideal gas D,, is independent of concen- 
tration to a first approximation [7], and in many 
greatly with 


other cases D,. does not vary 


concentration. Therefore, it seems reasonable to 


detine a new quantity, G,, by : 


G, = (11) 


st of axes, 
J (N, \ (3) 

— Also, for the other component, 
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and also, 


(12) 


Gy = 


Although it will be seen below that equations 
(13) and (16) to be obtained from equations (8) 
and (9) with L,, 
definitions of G, and G, given by equations (11) 
and (12) are generally useful and will be retained 


0 are not correct, the simple 


For an ideal gas, with L,, zero, 


(13) 


from here on. 


G, =G D,, RT 


For a real solution, not necessarily gaseous, 


dy, RT d\n a, (14) 


and a, is conveniently defined as : 


ay Yih (15) 
so that, 


G, Ly, dIn 


on 


G, RT 


for 0. 

Before considering the case where L,, 4 0, 
let us consider another approach ({8], [9], [10}). If 
du, /dx is viewed as a “ virtual force,” a mobility 
molecular 


G,, can be defined as the average 


velocity in the x-direction under unit force. 


molecules passing a unit 


Then the number of 

area per second is : 
G, (du, 

and so a flux in terms of moles is 


G, dy, 
17 

N dx (17) 
The coeflicient 
The 


mobility may be regarded as defining a flux of 


where N is Avogadro’s number. 
G, already introduced is merely G/N. 


molecules of one component relative to a set 
of axes fixed in the bulk of the fluid. Because 
of the obvious differences in mobilities of the 
molecules of the components of a gaseous or 
liquid mixture, and the existence of the 
KIRKENDALL effect ({11], [12]) equation (16) can- 
not be expected to hold. Since using a tracer is 
really a direct way (see below) of measuring 
the mobility, and since the mobilities found by 
then active “2” 


using first active “1° and 


as tracers are different, equation (16) would have 


Diffusion in binary mixtures 


the impossible result of giving two values of 
D,, for the same conditions. 

Let us now reconsider equations (8) and (9) 
with L,, #4 0 and see whether the term involving 
L,. can account for the observed differences in 


mobility. Using equation (11), the Gibbs-Duhem 


equation, L,, = L.,, and the fact that J, must 
equal — J,, there follows : 
[1 I 
G, = 1S) 
Yeo 


= Bi (20) 
4) 
where, 
d In 
| 
Om 
and so: 
dy, RT , dy, 
Wr (21) 
da da 


According to pE Groo' {1} and also, for 
Jost ({4], page 498) or ordinary binary 


so that: 


example, 


diffusion L,, Lie = Les 
Li. 


(22) 


and equation (20) reduces to, 


RT - RT G,/y, (23) 


Equation (22) follows from a mathematical 
treatment of the equations of irreversible ther 
modynamics ({1], pages 101-104). It is shown 
that if du, /dx is set equal to du,/dr, L,, equals 

L,., and since the L’s are not functions of the 
driving forces (du, /dx and dy, /dx), L,, should equal 
L,, for other values of the driving forces. Two facts 
lead one to question the applicability of this 
reasoning to ordinary binary diffusion. First, 
no physical situation exists for which dy, da: 
equals du, /dx except the trivial case where both 
equal zero. Second, in order to change the relative 
make 


them unequal) the concentration must be varied. 
it would 


magnitudes of dy, dx and di, dx (i.e. 


Thus, in order to retain Z,, 


By combining (8) and (18), 
J, = — (19) 
Yeo a 
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the L’s to be 
concentration in addition to driving force. 
For the 


a tracer “2” in a 


be necessary for 


measurements using 
the 


and Wy and 


sper ial of 


bulk of normal l ot 


same chemical species, = 


both equal unity, so according to equation (20) 


D,. RT G, RT G, 
coethcrent 
ill it 


However. 


Is rigorously ZeTo, D, is the 
let tis 


and | 


of self-diffusion in this case 


This is a quite reasonable result 


given the reasonable 


that 


physically assumption 
this case, 


further 


essentially equals lor 
satistied 


that 


(22) is rot 


equation 


doubt thre statement 


rite dow! 


J Pp j 


a2 


cise tissiot 


is of interest in the 
if one of the 


This case present 


components ( 2 Ss a tracer of 


the same chemical species as the normal ** ] 


3 represents thre 


If the 


homoge neous 


normal second chemical 
»an other 
muxture of and “3 
J, and J (), Also 


0. Therefore the equations (25) give 


species, tracer 1s introduced inte 


there 


results J, 


once more equation 
that 


19). with the added require 


ment 


“13 


Since OL 
Indeed, since L,; 


IS ZeTO ; L, is not necessarily 


zero. is actually the same as 


the cross coeflicient for the two normal! 
it should 


(1%) for 


spec ICs, 


not be zero, In addition quation 


~ 0 reduces to 


J, (dy, da) 


PY» (du, da 


independent of 


Equation (27) shows that the rate of diffusion 
of a tracer in an otherwise homogeneous binary 
G., 
a tracer of the other component 


The 


mixture gives a direct measurement of 


and by using 


G, can be obtained. ordinary diffusivity 


is then calculated from 


Ys 


The equations developed require a some what 


surprising vartation in G, and G, with Y). In 


particular at y, G.: the behaviour 


is not im accord with measurements on the 


and active gold in a 50 


diffusion of active silver 


pel cent Ag-Au 


given low che this case adequately 


illoy Ihe equation to 


and permit G, to differ from G,, at y, 0-5, as 1s 


observed xperime ntally. 
of the 


\ sible variation 


various coellicients as a function of Y, 


required by equation (28) is given in Fig. 1, 


ae : Schematic curves of the Variation of various 


coetlicients ording to equ ition (28) 


that 


here do 


It seems correct to state irreversible ther- 


modynamics as applied not 


vive a 


reasonable interpretation — of the situation, 


although there are at present few data on which 


to base conclusions. 


4 
; and, : 
“11 "22 “12 
9 
In a ternary mixture we 1958/ 
ad da la 
) 
J. | Yo L > (Po) 
ae 
‘ 
| 
2 
J, ~ RT (de, de) (27) 
: 2 
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Alternative ly ({1), [2]. [5], (13]) the whole prob- 


Now consider the case of the chemical diffusion 
of “1° through “3° where there is a known lem can be treated by leaving out any cross coeffi 


portion of active “1” present, which we shall cl nts such as L,.. ete. A flux p,. re lative to an axis 


eall * 2. In this case all three of the equations moving at the molar average velocity can be 
(25) are needed, and djg dx + 0. detined as : 
However, y, y. is a constant, so that : 
owe ‘ fy fo d, dD, p ly, dv) (36) 
dy, dif, 


and, noting that d. and J. are referred to the 


i 


y,da ya 


Same AXIS, 


Also, 


da J, b.) Ys (38) 


in a binary svstem., where Is positive in the 


where J, 4 1-0 positive r-direction and J, * Notice that 


the term d.) y, arises from the difference 


Thus : 


of the molecular mobilities ; it has nothing to do 


HO) with whether \, \ or not, J, may be 
(la Yo da ada 


viewed as the sum of ¢,. the basic diffusion flux 


and from equation (25) : arising from the gradient of chemical potential, 


y\ RT dy and p, $2) Yi» a flux caused by displace 
J, (Gry, p | ‘) (31) ment to satisfy J, J.. In addition, 
Ys da 
RT du \, (d, Ds) (N, N,) (39) 
J, Gos p Les ‘| (32) 
Ys) Wy di The term N, V, results from the bulk, hydro- 
se dynamuc flow of the system in which diffusion is 
J, (Ly; Los Gs p Y;) ahd yp, dy, (33) taking place ; it is influenced by the geome try of 
" da the apparatus and the way In which the density 


of the mixture varies with composition. The 
If, as postulated in irreversible thermodyna- : 

term (d, represents a flow on a molecular 

mics, the coetli ents L depend on the composition 

scale relative to an axis fixed at the “* centre of 


and not on the value of the driving forces. for 


moles of the fluid, i.e., an axis moving with a 
1 mixture having the same composition a la velocity related to the flux N, N,. rhe 


considered in the previous section it ts pe rmussible 


magnitude of (d, ¢,) is tixed by the mobili- 


to use the values of L alre addy determined. From 


ties of the diffusing molecules — i.e. by the state 


the previous results, L,, = 0 and L,, is given 


of the system, and it is not influenced by the 


by equation (29). Thus : 


external arrangements which cause hydrodyna- 


dy, mic flow. Possible relative magnitudes of these 
J, R1 Wy ( y, Ys rd P da (34) quantities are illustrated in Fig. 2. The arrows 
and, are vectors for the one-dimensional case being 
D,, RT Gy) Lis (35) considered. For the diffusion in a solid, 
Y Wa py the movement of inert markers is related to the 

flux (N, N, .). 
The value of D,, obtained from equation (28) By combining equations (1), (36) and (37) one 


and from equation (35) would be the same only if obtains : 


L,./py, is negligible, but since “1” and “2 
are chemically identical and have been assumed 
to diffuse at the same relative rates, it may be D dy, D dy, : p, a 


that the term actually is negligible. 


> 
5 | 
dia 
a and, 
| 
RT 
ne 
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pee case for N, different from No 


Stationary axis 


+a.) 


Fig. 2. 


the various fluxes at a point at which 4) O2 and 


by| >| dy). 


and, 


dD yD, 


If the flux ¢, is given by 


¢, = Gc, du,/da (42) 


as suggested by the discussion leading to ¢ quation 
(17), the quantities G, and G, can be determined 
by measurements with tracers of “1° and “2” 
in an otherwise homogeneous mixture. From 


(36) and (42), 


dD, RT yb, ( 43) 
and, 


D, RT G, (44) 


Thus : 
d In 


RT 


G,) (45) 
In contrast to equation (20), equations (41) and (45) 


certainly seem to require a variation of the various 


Vectors representing the relative magnitudes of 


which 
intuitively satisfactory, and moreover which is in 
accord with available data ({1 +}, [15]). A sketch 
of these curves is given in Fig. ¢ 


coefficients with concentration appears 


| 
| 
| 


seif diffusion 


Mole fraction 


Fic. 3. 


coeflicients with concentration according to equation (37). 


Schematic curves of the variation of the various 


CONCLUSION 


The equations proposed for steady-state, one- 


dimensional diffusion in a binary mixture are : 
(N, 

Dis p (dy, 

RT (y,G, + y,G,) 


where the mobilities G, and G, can be obtained 


N, =J, N,) 


(45) 


from measurements on the rate of diffusion of 


tracers of 1 and 2 in otherwise homogeneous 
mixtures. These equations are essentially those 
of Darken [2]. 

The basic diffusion equations involved in the 
from 


relations are obtained 


(40) and (42), giving : 


above equations 


J, GC, Yo (dp, (46) 


On the other hand, equations (8) and (9) from 
irreversible thermodynamics lead to : 


(47) 


ids 
| 2 ‘ | 
| 0,, self diffusion 
| RT! 
| l-y G. RT, 
(hod 
| | : 
° 
> N N, +N 5 
| = VOL 
; ly N.+N | 105R2/ 
| 
2 
| 
(41) 
3) 
dy, 
50 


[2 
[3] 
[4] 
[5] 
[6] 
7) 
[9] 
[10] 
(11) 
{12} 
[13] 
[14] 
{15} 


If pe 


or, using equation (18), 


R. B. 


Jost W. Diffusion in Solids, Liquids, Gases. 
Crank J. The Mathematics of Diffusion. 


K, 


DarkEN L. S. Trans. 


Ordnance 


Ge, dy, — Gyesy,| (48) 
ar Ye ar 
Groort’s assertion that is 
accepted, equation (48) becomes : 
Diy 
l dy,’ 
J, = (49) D 
dx da 


Although equation (46) seems to be in reason- 
able accord with experimental results ((14) and G, 
(15), for example), neither equation (48) nor (49) 
appear to be in accord with the results obtained 


in binary systems involving ordinary diffusion. Jj 
The conclusions which may be drawn from the 
above remarks may be summarized as follows : I 
ik 
1. Although equations (28), (48) and (49) 
are not in agreement with experimental N; 
data, only limited conclusions should be R 
drawn from this point. The Onsager re- . 
lations and equations (8) and (9) have a; 
proven their validity in many other cases, Cj 
such as the treatment of thermal diffu- z 
sion, to name only one. Mi 
2. Equation (45) and the related equations Me 
certainly can be viewed as owing much in ¢; 
their derivation to one of the basic ideas of 
irreversible thermodynamics, namely, that 
in diffusion there is a linear relation be- 
tween the flux and the gradient of chemical 
potential. This idea is at least as attrac- vj 
tive as that of a “ virtual force.” p 
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NOTATION 
coefficient of diffusion, em*/sec, defined by the 


equation J; Diy. p (dy; /dax) 
coefficient of diffusion, cm*/sec, defined by the 


equation ¢; D; p 
dha 
mobility, g mole em/sec-dyn” defined by the 
1 
equation ¢; Gy; “Ki. 
da 


Je» 
an axis moving at molar 


flux relative to an axis across which J, 
g mole /em? see i.e., 
average velocity 
coefficient in Onsager’s equations, 
(g mole)*/sec-dyn-cm* 

flux relative to a stationary axis, 
g. mole /cm? sec 

gas constant, dyn-cm/g mole-°K 
absolute temperature, “K 
activity 

concentration, g mole em 
distance, cm 

mole fraction 

activity coefficient 

chemical potential, dyn-cm/g mole 

flux relative to at the 
average velocity and related to the molecular 


Gc; (du, /dx), g mole /cm? sec. 


an axis moving molar 
mobility by >; 
It accounts for only part of J;, to which it is 
related by equation (40). 
din Yi 

diny; 


density, mole /em® 
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Heat transfer between a fluidized bed and a horizontal tube 
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Abstract— Measurements of coeflicients of heat transfer from fluidized beds to horizontal water 
cooled tubes are discussed The bed diameter was 0-565 m The following quantitities were 
varied (1) bed temperature velocity of the fluidizing air (3) mean particle diameter 

(4) particle shape; (5) particle density: (6) water-tube diameter. The second third and fifth 
quantities were varied down t ich lower values than before | the sixth quantity had not 


previously been varied 


Phe results of old and new experiments are correlated by two different formulae, one for the 
range where the viscous for nm the particles are predominant and an analogy between the 
fluidized bed and a flowing gas or li juid is assumed, and another for the range of coarser and 
heavier particles where inertia effects prevail. The distinction between the two ranges is 


characterized by the value of Reynolds group. The limitations of the applicability of the 
correlations are discussed 


Résumé— Lo auteur discute des mesures des coefficients du transfert de chaleur a partir de lits 
Hluidisés vers des tubes refroidis 4 eau. Le diamétre du lit était de 0-565 m. On t fait varier les 
facteurs suivants : (1) température « 2) Vitesse massique de lair fluidisant; (3) diameétre 
moven des particules; (4) forme des particules: (5) densité des particules; (6) diamétre des tubes 

“uu Les facteurs sous (2), (3) et (5) ont été variés jusquia des valeurs beaucoup plus basses 


qu auparavant [1 Le facteur sous (6) navait pas été varié avant 


Les résultats des expériences nterieures et presentes sont mis en corrélation a aide ce deux 
formules differentes, une pour ntervalle of dominent les forces Visquetises agissant sur les 
particules et ou Pon admet qu thexiste une analogie entre le lit fluidisé et un courant de gaz ou de 
liquide, et une autre pour Pinter lle des particules plus grossiéres et plus lourdes of ayissent des 


effets d inert ic Lua distinetion entre les deux est curacterisce par la valeur d'un he vnolds 


On discute les limitations de lapplicabilité des corrélations 


Zusammenfassung Messunyven der Koeflizienten der Warmeiibe rtragung von Fliessbetten an 
waagerechte wassergekihite Rohre werden be sprochen. Der Bett-Durchmesser betrug 0.565 m 
Folgende Folgende Grossen wurden veriindert : (1) Bett Temperatur, (2) Luftgeschwindigkeit 


(3) Mittlere ‘Teilchendurchmesses Teilchenform, (5) Teilehendichte, (6) Durch sser des 
W asserrohres 


Die zweite, dritte und fiinfte Grosse wurden bis zu viel kleineren Grossen vartiert als bisher 


dlie sechste Grosse wurde bisher nicht yeindert 


Die Ergebnisse der alten und neuen Versuche wurden durch zwe i verschiedene Gleichungen 
korreliert, cine fiir jenen Bereich, in dem dic Zihigkeitskrifte an den Teilehen vorherrschen 
und eine Analogie zwischen dem Fliessbett und einem strémenden Gas oder ciner Fliissigkeit 
ingenotmmnen werden kann, und eine andere fiir den Bereich gréberer und schwererer Te ilchen, in 
dem die Trigheitswirkungen vorherrschen. Die Unterselh idung zwischen den beiden Bereichen 
ist durch den Wert einer Reynolds-Zahl gvekennzeichnet Die Grenze der Anwendbarkeit det 
Beziehungen werden besprochen 
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Heat transfer between a fluidized bed and a horizontal tube 


l. INTRODUCTION investigated. The diameter of the tube was not 


A FEW years ago, the results of some measurements varied. Many data were obtained on the influence 


of the coefficient of heat transfer from a large- of bed temperature and mass velocity of the 


scale fluidized bed to a horizontal water-cooled fluidizing air. It was found that the heat-transfer 


tube were presented [1]. Beds consisting of — coeflicient was independent of the air pressure when 


various kinds of coarse sand and iron ore were the linear velocity of the air in the bed remained 


Table 1. Experiments with a horizontal tube 0.85 m above the sieve plate. The height of the bed was 1-2 m 


in the non-fluidized state the bed diameter was 0-565 m 


Vinimum valu 


Vean of product of mass 
vel Weight Particle relocity and kine Partich Outer Vass Curve in 
Material surface mean density matic viscosity at shape diameter velocity Figs. 1, 2. 3. 
party li diameter \3 chich fluidization of tebe of air 5. and 6 
diameter (3 was observed 


ni=se 


ores 


ors 


O152 


O42 


O-O153 i 


Reyenerated 
eracking 63-5 Os 1.60 O-065 Sharp 


catalyst 


140 


oO510 0-006 J J 


O457 


O-O169 OOTTI 


Sand used 
for sand- 187 0-42 Sharp O-0459 T 7 


blasting 


Dune sand 213 20 2 O65 Round O-0336 00-0562 \ \ 


aoe 
161 
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constant. The influences of bed-height and location 
of the than 20 


per cent. 


diametrical tube were smaller 


Since then some complementary measurements 
with much finer and lighter particles have been 
made, and stainless-steel tubes with various 
diameters have been applied. The bed te mperature 
and the mass velocity of the air were varied in the 
same way the bed 


height, and the position of the tube were kept 


as before. The air pressure, 


constant. 


In what follows old and new results will be 


correlated by two different formulae, one for the 


J/m* sec °C 
h 


300 °C 


b 


Fic. 1. Coefficients of heat transfer A from beds of catalyst 
at various mass velocities to horizontal tubes with different 
diameters, as functions of the bed temperature 4, (see 


Table 1). 


range of fine and light particles, and another for 
the range of coarse and heavy particles. The 
transition between the ranges of validity of the 
two expressions will be discussed. 

In a subsequent publication similar correlations 
for heat transfer between a fluidized bed and a 
tubx 
already published by us [2] and those of other 


vertical will be derived from the data 


investigators. 


2. Hear-TraNsrer MEASUREMENTS 


The experimental setup, the test procedure, 


and the method of dealing with the data have 


J/m* sec °C 


200 300 °C 


Fic. 2. Coeflicients of heat transfer A from beds of sand 
at various mass velocities to horizontal tubes with different 
diameters, as functions of the bed temperature t, (see 


Table 1). 
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Table 1 
marizes the new series of experiments, the results 


been described previously [1]. sum- 


being shown in Figs. 1 and 2. The increase in the 


heat-transfer coeflicient h with increasing bed 
temperature or mass velocity has already been 
discussed {1}. A similar increase is found at 
constant mass velocity, when the tube diameter is 
decreased ; this effect is analogous to the well- 
known increase in the coeflicient of heat transfer 
between a tube and a fluid flowing normally to 


that tube with decreasing tube diameter. 


3. PRESSURE MEASUREMENTS 


Pressure measurements were made in the same 
The 


obtained with the four different horizontal tubes 


way as previously described [2]. results 
did not appear to differ. 

With all kinds of sand the pressure drop per 
unit length proved to be independ nt of the loca- 
tion in the bed, and the porosity « could be derived 
from it in the same way as before [2]. This pro- 
cedure was justified because the total pressure 
drop was equal to the total weight of the bed per 


unit of area. Some results are given in Fig. 3. 


€ 050 T T 


. 0 100 200 300°C 
Fic. 3. The porosities « of beds of sand at various mass 


velocities, as functions of the bed temperature 4, (see 
Table 1). 


The absence of any measurable difference be- 
tween the two curves observed with dune sand 
may be owing to experimental error. 

With catalyst the pressure drop per unit length 
was by no means independent of the location in 


Heat transfer between a fluidized bed and a horizontal tube 


the bed. The results obtained in the upper part 
of the bed agreed in some measure with those 
obtained with sand, but in the lower part some 
remarkable features were observed. The apparent 
values €, of the porosity, derived in the usual way 
from the pressure measurements in the lower part 
of the bed, decreased with increasing mass velocity 
G and, at low values of G, also with rising bed 
temperature. It was concluded that €; did not 
represent the true porosity ¢«, the pressure drop 
being unequal to the weight of the bed per unit of 
area. All these anomalies could be ascribed to a 


combination of two phenomena : 


Segregation of the catalyst particles, so that the 
majority of the coarser particles will move to 
the lower part of the bed, and the majority of 
of the finer particles to the upper part. This 
effect will cause the porosities « and ¢«, to be 
higher in the upper part of the bed than in the 
lower part ; the difference will decrease with 
increasing mass velocity G owing to better 
mixing of all particles throughout the bed. 
It was found, however, still to exist even at 
the highest values of G applied in the present 


investigation. 


Channelling is an effect which has been described 
by various authors [4, 5]. Leva [4] states that it 
is promoted by low mass velocities and low 
mean particle diameters, by large bed dia- 
meters, and by the presence of moisture. All 
but the last of these factors are favourable to 
the occurrence of channelling in our experiments 
with catalyst. It causes a decrease in the 
pressure drop, which becomes lower than the 
weight of the bed per unit of area [4, 5] ; differ- 
ences up to 28 per cent were observed. Accord- 
ingly, «, will become higher than «. The main 
reason is that a great part of the particles 
are supported by the sieve plate instead of 
being borne by the flow of fluidizing gas. It 
follows that the effect is most important in the 
lower part of the bed, and that it decreases with 
increasing mass velocity G. At the highest 
value of G applied it was absent even in the 
lower part of the bed, since the total pressure 
drop was found to be equal to the total weight 


of the bed per unit of area. 
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The 


horizontal tube 


distance between the sieve plate incl the 


was so large that channelling did 


not occur at the location of the Cherefore. 
the local value of « could be derived in the usual 
wav from pre ssure measurements in thy r 


part of the bed. Results are given in Fig. 4. No 


€ 060-—— T T T 
x 
| 
x 
| ail 
2 0S > kg * se 
Fic. 4. The porosity « of the upper part of a bed of 


catalyst as a function of the mass velocity G 


svstematn influence of the hed temper iture On 


was obse rved possibly brew tine the normal 


with mnecreasing bed temperature 


compensated for a decrease in’ the effeet of 


segregation. The large influence of t] 


to t hee 


Tass 


Ve low itv is duc range ol of 


covered in the « xperime nts with catalyst 


CORRELATION 


Since the experiune ital data obtained with fine 


parti les, such as catalyst. and with tubes of 


Various diameters, cannot be re presented by an 
EXPression as simple as the correlation of previous 
results [1]. a discussion of new correlation 
formulae covering old and new data will follow. 

It has proved to be worth while to describe the 
new results with fine and light parti les. namely 
catalyst, sand used for sand blasting, and dune 
those with coarser and 


sand, s« parate ly from 


heavier materials. This distinction between two 


groups ol particles may be connected with a 


difference in the nature of their motion in the gas 


currents, and the assumption will be made that 
the tine and light particles will almost exactly 


follow the paths of the gas, whereas the coarser 
and heavier particles, owing to their inertia, will 
have a motion that is less dependent on these 


paths. Of course, such an assumption may over- 
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simplify the picture, but a reasonable de scription 
of the experimental results obtained with fine and 
light particles may be derived from it. 

If the paths of the fine and light parti les do 
not differ from those of the gas. the mixture of 
fluid 
fhuidizing gas. 


hn al, but the 


gas and particles can be considered as one 
moving in the same wav as the 
This cannot be true throughout th 
greater part of the horizontal water tube will be 
in the both the 


gas are moving upward [2], so 


mner zone of the bed where 
particles and the 
that there it may be correct to consider the 
mixture as one fluid 

Coetlicients of heat transfer between a tube and 
a gas or licqguid flowing normally to that tube can 
Nusselt group 
divided by Pr,, as a function of the Re vnolds 
|. Therefore the analogy 


between a fluidized bed of fine and light particles 


hie represented by plotting the 


group [6 assumed 
and a fluid could be checked by attempting to 
represent the experimental results obtained with 


these particles by a re lationship hetween : 


AD, GD, p, (1 €) 
, and 


(C je kr) pile 


The Reynolds group GD, p €) con 


tains the bed density which, since the contribu 
tion pe of the gas is negligible as compared with 
the contribution p. (1 ©) of the particles, ts 
equal to the latter contribution, and the inter 


stitial gas velocity G pe which has to be used 


since, in general, the parti le velocity ina fluidized 
bed is unknown. Similar arguments apply to the 
other quantitne s occurring in the Reynolds group, 
the Nusselt AD k, and the Prandtl 


yroup Cy 


Fig. 5 proves that it is indeed possible to give, 
by means of a relation between the two groups 
considered, a fair re presentation of the results of 
all the measurements with fine and light particles 
specified in Table 1. The experimental curves do 
not deviate by more than 42 per cent from the 


straight line 


Al 0-66 (GD, (1 
(¢ | pe } 
\ good correlation of all the data for coarse and 


heavy particles, namely the new type of Maas sand 
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Correlation of results of experiments with fine and light particles (see Table 1) 
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Fic. 6. Correlation of results of experiments with coarse and heavy particles (see Table 1 and Ref. [1], or Table 2) 
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and the materials previously specified [1], was 
obtained by adding the dimensionless group : 


D» 


to the two groups mentioned above. This group 
will be discussed further in the next section. It 
argued that a porosity like 
(1 e)/e should be but it found 


that the introduction of a porosity term into the 


might be term 


added, was 
correlation did not improve the description of the 
experimental data for coarse and heavy particles. 
It was therefore omitted from the third group, 
and also from the Reynolds group. The re pre- 
sentation of the results of all the measurements 
with coarse and heavy particles by means of a 
simple relationship between the three groups 
The 


not deviate by more than 29 per cent from the 


thus obtained, is given in Fig. 6. curves do 


straight line : 


AD, k 
(¢ 


10 GDip, 


If the mean particle diameter D,, is increased 
and it is required to maintain the same bed tem- 
perature, the same mass velocity G, the same 
tube diameter D,, the same particle density p,, 
the same specific heat of the particles c, a similar 
particle shape, and a similar particle size dis- 
tribution, there will be a decrease in the por sity €: 
hence the correlation for fine and light particles 
predicts an increase in the heat-transfer coefticient 
h. According to the correlation for coarse and 
heavy particles, however, the heat-transfer co- 
efficient A is proportional to — so that it 
should decrease. 

When D, is increased and the degree of fluidiza- 
tion* is to be kept constant while the same bed 
temperature, the same tube diameter, the same 
particle density, the same specific heat, a similar 
particle shape and a similar particle size dis- 
tribution are maintained, it will be necessary to 
increase the mass velocity G. At the same time, 
the increase in D, will cause a smal! decrease in 

* The degree of fluidization can be defined as the ratio 
of the product of mass velocity and kinematic viscosity 
of the fluidizing gas to the minimum value of this product 
at which fluidization can occur [1). 


the porosity «. For these two reasons the correla- 
tion for fine and light particles again predicts an 


increase in A. However, since the effect of the 


Increase in D» will exceed the effect of the 


increase in G, the correlation for coarse and heavy 
particles again predicts a decrease in h, 

When the particle density p, is increased and 
the other properties of the particles, the bed 
temperature, the mass velocity, and the tube 
ke pt 


decrease in the porosity because of the decreasing 


diameter are constant, there will be a 


degree of fluidization. For fine and light particles, 
both the Increase In p, and the decrease in « will 
cause an increase in hk, For coarse and heavy 
particles, however, the heat-transfer coeflicient is 
proportional to p,°*, so that it will decrease. 


When p, is 


properties, 


increased and the other particle 
the bed 


temperature, and the tube diameter are to be 


the degree of fluidization, 
kept constant, it will be necessary to increase G 
8}. 


this means that 


proportionally top, Since p is 
negligible as compared with p,, 
light 


coetlicient is proportional to p, 


for fine and particles the heat-transfer 


will 


* so that it 
increase. For coarse and heavy particles, however, 
no change in h is to be expected. 

It should be borne in mind that the effects of an 
increase in D), or p, on h were derived on the 
assumption that the increase would not involve 
the light 
particles to that of coarse and heavy materials, 
This the 


section. It follows from the above that 


a transition from range of fine and 


limitation will be discussed in next 
under 
definite conditions maximum heat-transfer co- 
eflicients can be obtained by choosing values of 
D),, and p, that belong to the transition between 
the two ranges. 

Similar statements about the effects of other 
quantities contained in the correlation formulae 
can easily be derived. 

VALIDITY OF 


RANGE OF THE 


CORRELATIONS 
After a collision between two particles, the 
resulting velocities may be in any direction. 
When external forces are absent, they will retain 
their directions up to the next collision. When, 
however, the ambient gas moves relatively to 
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the particle, a force will be exerted on the particle values the paths of the particles will be bent 
owing to the viscosity of the gas and it will bend towards the direction of the gas flow, and the 
the path of the particle towards the direction assumption under which the correlation for fine 


of the gas flow. This force acts on the particle and light particles has been derived may be 


along a distance which, on an average, is equal expected to be valid. This conception is com- 


to the mean path between two successive colli- patible with the mechanism of heat transfer 


sions. The corresponding work is proportional proposed by Micktey and Fairsanks [9] and 
to the product of this path and the force and, if it Kramers [10]. For high values of the Reynolds 
is assumed that the particle velocity is propor- group the correlation for coarse and heavy 
tional to the interstitial gas velocity G pe,theratio particles may apply; it would appear that in this 
between the kinetic energy of the particle and case the group ,* Dy? p, 8 affects the stirring 


this amount of work can be shown to be propor- factor introduced by MickLey and FainBanks [9]. 


tional to GD, p, pu. The porosity « does not The group po g is the ratio between the 


occur in this ratio. Froude group G*/D,, p* g [11] and the square of 


Since the Reynolds group GD,p, py is a the Reynolds group GD,, p, py. 
measure of the ratio between the actions of The range of the values of the Reynolds group 


inertia and the viscous forces, the latter will be pertaining to all the relevant experiments is 


predominant for low values of the group, and given in Table 2. A survey of these values leads 


the inertia effects for high values. Thus for low to the conclusion that the correlation in Fig. 5 


Table 2. Applicability of the two correlations 


Mean Range 


Range of 


volume- Range of bed Results Curves in Correla 
Material surface Particle of mass tem- in Fig. relevant tion in 
particle density velocity G perature group a D pps or in figures Fig. no 
diameter t, py Ref. no. 


D, Ps 


(m 10° (kg/m*® 10%) |) (kg /m*sec) 


Regenerated 
catalyst 1-60 0-010-0:149 40-260 740 Fig. 1 |G,H, I,J, K 5 


Sand used for 


sand-blasting 136 2-66 0-046. 0-077 40 200 670— 1290 Fig. 2 5 
Dune sand 213 2-66 0-056-0-096 40-240 1270- 2550 Fig. 2 A,B 5 


New type of Maas 


sand 316 2-66 0-099-0-163 20-260 3360- 6610 Fig. 2 C.D 6 
Sharp silver sand 215 2-66 0-089-0-:151 40-260 2050- 4060 Ref. [1]| E, F,V,W 6 
Old type of Maas 

sand 353 2-66 0148-0242 40-220 5620-10900 Ref. [1] oO, P 6 


Mixture of silver 


and old type of 
Maas sand 260 2-6 0-118-0-198 40-340 3270— 6690 Ref. [1] Q, RLS 6 
Round silver sand 259 2-66 0-091-0-:150 40-280  2520— 5500) Ref. [1] M,N 6 


Ground iron ore, 


wide sieve frac- 
tion 130 515 0-234 40-220 6160— 7250 Ref. [1] XxX 6 


Ground iron ore, 


more homogen- 


' 
~ 


eous fraction 247 515 239 40-240 | 12900-15100 Ref. [1] Y 6 
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is applicable when the value is lower than 2050, 
and that the correlation in Fig. 6 applies when it 
2550. Values of the 


2050 and 


is higher than Revnolds 


group between 2550 are included in 


both correlations. For this transitional range it 
is recommended that the average of the amounts 
the 


good prediction of the heat-transfer 


consid red a 


coctlicient A. 


given by two correlations b 

Since D, and p, occur in the numerator of the 
Reynolds group, it is obvious that the criterion 
given can be considered as a distinction between 
fine and light parti les on the one hand. and coarse 
and heavy particles on the other. However, this 
conception Is too crude because, for instance, 
even with catalyst the value of the group may 
exceed 2050 if extremely high values of G or low 
In these cases the 


procedure should bn followed as wit 


values of p are applied. san 


coarsctT and 
heavier partic les, since it is not permissible to 
neglect inertia effects. 


Since the correlation in Fig. 5 has been derived 


from an assumed analogy betwee he thuidized 


bed and a flowing gas or liquid, it is recommended 
that it be applied in all cases where a coctlicient 
of heat bed 


diametrically inserted tube has to be calculated, 


transfer between a tlhuidized and a 
so long as the value of the Revnolds 


than When the bed 


than 0-2 m the value of « 


yroup Is le SS 


2050. diameter is greater 


occurring in the cor 
relation can be 


Figs. 3 and 4. 


from the data 
Leva [7, 12, 13, 14] 


derived eiven mn 


presents data 


The 


influence of the bed diameter on the porosity has 


for bed diameters smaller than 0-2 m. 


already been discussed [2, 15]. 


The dimensionless yroup un 


g has been 


introduced into the correlation in Fig. 6 only 


because a good ck scription of the ‘ xperimental data 
The use of this correla 


could thus bn obtained. 


tion can therefore safe ly be recommended only 


when the value of the Reynolds group ts higher 
than 2550 and when, moreover, all the variables 
are within or nearly within the range covered by 


the « x perime nts. 


NOTATION 


specific heat of particles at bed temperature (J kg 'C) 
mean volume-surface particle diameter [3 

outer diameter of tube 

mass velocity of thuidizing was 

acceleration of gravity 

coeflicient of heat transfer between bed 
und tubs (J 
thermal conductivity of gas at bed tem 
perature (J see mC) 
prandtl group 

mean bed temperature 

absence of 


fraction of voids (in the 


channelling « ) 


4 
apparent fraction of voids, derived from 
pressure drop 

dynamic viscosity of gas at bed temper- 
ature 

density of gas at bed t« mperature 


density of particles 
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The problem of smoothing out concentration fluctuations by passage through cascades 


of stirred mixers is formulated from the point of view of probability. The formulation is in terms 


of the theory 


side by side with spectral ten hniques The two methods are illustrated by 


particular problems, some of which have 


and others. 


Résumé— Les probléme de la repartition 


passage a travers une série de mélangeurs agités, est formuk 


of random processes, and techniques based on 


previously 


uniforme des 


covarhince functions are introduced 


ipplication to various 


been treated by Danckwerrs ([4 5}) 


fluctuations de concentration, par 


du point de vue de la probabilite 


La formule est établie a partir de la théorie des re partitions desordonnées. Des techniques ayant 


wur base les fonctions covariantes sont présentées en relation aves 


Les deux méthodes sont illustrées par leur application a divers problémes ; 


ayant été discutés par Danckwenrrs [4 


Problem 


Rithrekesselkaskaden auszugleichen 


Zusammenfassung — [Das 


Konzentrationsschwankungen 


des te« hniques sper trales. 


certains dentre eux 


et dautres auteurs. 


beim Durchgang durch 


wird vom Standpunkt der Wahrscheinlichkeit formuliert. 


wird clit The or cler Zufallsproze ase he ranyvezoven una he he Tithve n, die auf kovarianten 


Funktionen beruhen, werden 


Methoden werden an verschiedenen 


durch Daneckwerts ({4, und andere 


I INtRopUCTION 


Tuts paper studies the performance of cascades 
of stirred tanks designed to smooth out random 
concentration fluctuations in a continuous process 
stream. The point of view adopted is Statist- 
ical in the following sense: the concentration 
history of some material in a stream is considered 
not as a single function of the time ¢, but rather 
as an ensemble of such functions indexed by a 
parameter «. Such an ensemble of time histories 
IS called a random process. The ensemble index 
x may be thought of as ranging over duplicate 
pieces of equipment operating under the same 
statistical conditions. 

The statistics arise in establishing probability 
Thus, 


distributions over the ensemble index «. 


for an ensemble of concentration histories 


e (t, a), one has the univariate probability density 


c 


zusammen mit) Spektralmethoden 


behandelt worden sind. 


angewandt Diese beiden 


\nwendunygsfillen demonstriert, von denen einige bereits 


t) da 


the probability that a c(t, x) 


da 


the fraction of all possible x for 


which < c(t, a) <a da 
the bivariate probability density Q., 


ys 7; t) da dy 


c(t, a) 


the probability that 2 
dx and y c(t + +. a) 7 


dy simultaneously 


the fraction of all possible « for which a 


c(t. x) r-+dzx and y < c(t rT, a) y 


dy simultaneously 


and so on. 

We shall be concerned here not with the full 
scope of the probability distributions P, @, etc., 
but only with their first and second moments. 
We consider first the mean value of the random 


process c (t, «), 
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m, (t) E [e (t, x)] 


t) da 


E, 


value, denotes an averaging over the ensemble 


where the operation taking the expected 


index «. A measure of the variation of ¢ around 
its mean value is furnished by its variance : 


m, (t) 


m,(t) |? P.(a; t) da 


Primary information about the statistical rela- 
tionship between values of ¢ at different times ts 


furnished by the auto-covariance function : 


m, (t)} 


mu 


lly 


It may be noted that: 


» (t) 


In mixing studies one has to consider not one 
random process ¢ (f, ~) In isolation, but rather two 
such processes ¢ (t, a), c' (t, a) jointly. The pro- 
cesses ¢ and c' are, in general, the concentration 
ol some key component entering and leaving a 
mixing system. The joint statistics of ¢ and c' 
generalizations of the 


are defined by suitable 


P_.(a;a';t)da dx 


distributions / Thus, one considers : 


the probability that 2 


(t, a) a dz and x' c' (t.a) 2, 


simultaneously 


ete. information about the statistical 


re lationship between values of ¢ and ¢! at different 


Primary 


times is furnished by their cross-covariance 


function : 


m, (t); 


(tf 


r)}] 


It will be important in what follows to single 


>» &) 


out from the general class of random processes 


those which are called stationary. The stationary 


62 


those for which the probability 
P, Q, ete. are independent of the 


' that is, those processes ¢ for which: 


processes are 


distributions 
origin of time 


— P,(a), and 


@. @. 7) 


etc., independent of «. For a stationary process 


c, the mean and variance are constants : 


of ¢, and the auto-covariance 


function depends only on the lag +: 


independent 


independe nt of ¢. 

If two random processes, ¢ (t. x) and c' (f,.a) are 
separately stationary, then their separate prob- 
the 


.on the other, are se parately 


ability distributions on one 
hand, and P., 
ndent of t, with the conse quences noted in 
the last paragraph. If, further, their joint pro- 
bability distributions ete. are independent 
of t, the two processes may be said to be jointly 


stationary. As a consequence of this, their cross- 


covariance function depends only on the la 


0) (r) 


independent of t. 


In the body of this paper, the criterion of 
performance of a given mixing system will be 
taken to be the variance of a key concentration 
the outlet 
not because it is thought to be the most realistic, 
but 


mathematically. 


in stream. This criterion ts adopted 


because it seems to be the most accessible 
derived 
of the outlet 
stream concentration in terms of the covariance 


See- 


Expressions will be 


in typical eases for the variance 
structure of the feed stream concentration. 
tion IT will review the transient performance of 
Section 


IIl will discuss a typical problem with non- 


cascades of continuous stirred mixers. 


vive 


stationary feed stream concentration, and g 


the design relationships between the inlet and 
outlet 
Section IV will discuss a similar problem with 


concentration covariance functions. 


stationary feed stream concentration, and show 


ere 
2 
= pe (7) 
Tf 
158/ 
pet r)}] 
i 
Pee* 
: 


how the design relationships can be simplified 
by introducing the spectra of the inlet and outlet 
concentrations. Section V will consider a problem 
in discrete, rather than continuous time, which 
can be handled within the same statistical 


framework. 


Il. Cascades oF STIRRED MIXERS 


Consider a cascade of continuous well-stirred 
tanks — the first feeding the second, the second 
feeding the third, and so on, all at the rate 
Wt® hr. Let the volume of the nth tank be 


Vy, ft®, and its nominal residence time therefore, 


0 V/W hr. 


n n 


A straightforward material balance then shows 
that the concentration ¢, (t) moles ft® of some 
material in (and leaving) the nth tank satisfies, 
as a function of time, 4, hr, the differential 
equation 


dt 0 


In (1), the function ¢, ¢, (t) is the concentration 
of the material in question in the feed to the 
first tank. 

We want from (1) information as to how varia- 
tions in ¢, (¢) are propagated through the cascade 
to produce variations in ¢, (t). But this informa- 
tion should not include extraneous transients 
such as arise, for example, in start up. To 
achieve this, imagine the cascade to have been 
started up mn the remote past, and solve the 
differential equations (1) under the initial 


conditions. 
¢, (t) bounded as t + Opn 


Under suitable regularity conditions on ¢, (t) the 
mathematical system (1), (2) has the unique 


solution : 


Cy (t) = | H,, (x) eg (t — 2) de (3) 
where H, (x) is the convolution of elementary 
exponentials (cf., e.g., [1]) 


1 
H,, (2) e e e */"n (4) 
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that is: 
1 (x) 0, ¢ (5) 
with 


1 
H.(a H, (x)* 
» (2) 1 r) 0, 


| H, — u) du 


0 


Os 
(a u) u/% 
é au 


0 


and so on. It might be pointed out that in the 


case of equal residence times in all stages : 


6, = 0, -=@ 7] 


n 


the kernel i, (2) has the explicit form : 


r z/@ 
H, (2) = (6) 
(n 1)!@ 


In Sections III and IV we shall treat the inlet 
and outlet concentrations as random processes 
in the sense of Section I, Accordingly we introduce 
their dependence on the ensemble index « 
directly into (3) to give : 

(t, a) = | H, (a) (t — a, «) dx (7) 


0 


This relation is of the general form studied 
extensively in certain control systems applications 
(cf., e.g., [2]), where the inlet concentration So 


is thought of as a generalized driving force : 
Cy (t, a) = (t, x) 


and the outlet concentration ¢, aS a generalized 
response, 


c,, (t, a) = (t, a) 


For reasons of technical convenience, we write 
the connexion (7) in the form : 
c' (t, a) K (x) c(t —a, «) dr (8) 


where : 
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In (9) the function H, (2) is given in general by 


(9) 


(4), with the important spe ial cases (5), corres- 


ponding to only on mixing stage, and (6), 


COTTE sponding to n identical stages. 
The kernel AK (7) in (8) may be characterized 


by its Fourier transform (cf., e.g., [3}) 


r 
| A (rye 


By (9). 


L 


with the particular forms (cf. [1]) : 


. for the i, 


" 


. for the i, 


_» for the i, 


The function L(qw) is the transfer function 
of the system which transforms ¢ into ec’ (ef. [2}). 
Section III 


non-stationary inlet concentrations c, and Section 


will discuss the relation (8) for 


LV, for Stationary ¢. 


We now consider the feed to a mixing system 
a key 


where the 


Non-STATIONARY INLET CONCENTRATIONS 


to have component with concentration 


c(t, a), random process ft. a) has 


mean and auto-covariance function 


31), as defined in Section L. It may then be 
seen from (8) that the random process c(t, a) 
describing the outlet concentration of the key 


component has mean : 


m,* (t) K (x) m, (t 


vr) da (12) 


cross-covariance function with the inlet concentra- 
tion : 


| (x) p. (+ 


r:t)da 
and auto-covariance function : 


t) | A (x) 


r) da (14) 


To illustrate the use of these formulas we con- 


sider the feed stream to be made up in batches 


in such a way that the concentration c of the key 
component suffers a discrete change every h hr. 


imagine the batch concentrations 


Specitically, we 


to be drawn independently from a population 


with mean m, and variance o,*, both independent 
oft. If we Ss) nchronize the change-over times with 
then 


the auto-covariance function of the feed stream 


the integral multiples of h on the t-seale, 


concentration ¢ Is: 


iftand t r lie between 


the 


tiples of h 


successive mul- 


otherwise 


Finally we imagine the feed stream passed through 
a single stage of mixing with residence time @ hr 


so that, from (9) and (5), 


and ask how large we should take @ to keep the 
of the 


assigned level. 


Variance outlet concentration c' below an 


Wi that, 


kernel (16) is 1, the constant mean concentration 


note tirst since 


the integral of the 
mm. passes unchanged (as it ought) through the 
mixer ; that is, by (12), m,, m.. Next, inserting 
(15) into (13) gives the cross-covariance function 


in the form : 


a, h) (17) 


where : 


|_| | i, da 
VOL. 
- 
l+iwd, 1+iwd, | 
| | | | 
| 1+ | 
l 
a 
‘ 0 ; u 0 


In (17), 7 is the fractional part of t/h: that is, 


(18) 


where [=] is the largest inte ver not exceeding 


Substituting (17) into (14) gives the 


variance of 


the outlet concentration as: 


(t) 


2nh/e\ 


tanh é (19) 


Ze 


In (19), » is given by (18) and varies between 


Oand 1. The 


at the endpoints of the 7 


maximum value of (19) is assumed 


interval, giving the 


(best) inequality 


Thus, to ensure that for all ¢ 


(7) 


between O and choose 


for some assigned é 


to make 


Oh 
are tanh 


Similar calculations can be made for cascades 


of mixers, but the comple xity of the calculation 


mounts rapidly with the size of the cascade. 


IV. Srarionary CONCENTRATIONS 


If the random process ¢ (¢, «) describing a key 
feed to 
stationary in the sense of Section I, then general 
(cf, that the 


a mixing system ts 


concentration in the 


considerations show corres 


ponding outlet concentration a) given by 


(8) is also stationary, and that the inlet and outlet 
One 


are jointly stationary. 
that the 


proce SSCS, ¢ and c', 


can of course § se¢ directly constant 


stationary ¢ unchanged 


mean mM, ol a passes 


(as it ought) through the mixing system ; from 


(12), 


| A da 


‘ 


since, by (10) and (11), 


| K (x) da L (0) l 
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Also, the general covariance relationships (13) 
and (14) show directly that the functions Pees 
and p,, are independent of their second arguments 


if p, is thus : 


The form of the equations (20) can be con- 
siderably simplitied, and the tedious calculations 
in solving problems such as those in Section HI 
in some measure avoided, by introducing Fourier 
transforms, when the y exist, of the covariance 


functions p, 


The functions S, and S. are the power spectral 
densities of ¢ and c', respectively; the function 
is the Cross-power spt etral ce nsity of ¢ to 
The 


to express the covariances in terms of 


transforms (21) can be inverted 
(et. 


these Spe ctral densities : 


T) | S_(w)¢ des 
c 


(a) ¢ dw (22) 


| dw | 


The convolution form of the equations (20) 
permits writing them (cf. [3]) in terms of the 


spectra S as follows : 


(w) L (w) S, (w) 
(23) 
So (w) L (w) (w) 


where L(w) is the transfer function defined in 
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: 
| 
h 
| 
p rT) | A (2) p r) da 
(20) 
t 
Pp, ) | A (xr) p r) da 
VOL. 
P 
9 
Te S, (w) 
‘ 
(aw) |p, (r)e "adr > (21) 
l 
S io (r)é¢ ar 
0, 
| 


(aw) (w)/* Ss, (aw) 


and from (22): 
Yr 


| IL (w)|? S, (w) ¢ 


Per (7) 


so that taking + 0 gives : 


On | iL 
In (24), |L (w)|* 
real 


(w) dw (24) 


is the sum of the squares of the 


and imaginary parts of L(w), as will be 
illustrated in the sample problem to follow. 


To illustrate the consider, as 


Ill, a feed 


stream in which the fluctuations in concentration 


use of (24), we 
in the sample problem of section 
cofa key component are to be smoothed out by 
If this mixer 
then the 
function L (w) is given by (11) as: 


passaye through one mixing stage. 


has residence time @ hours. transfer 


(w) (25) 
so that: 
LL (26) 
l+ & 


The transfer function (25) corresponds to the 
kernel (16) om the sample problem of Section IIT. 
We take the inlet concentration ¢ to have variance 
a.” and characteristic time lag A hr for the ce Cay 
of the dependence of neighbouring values of ¢ 


on each other: specifically, we take 


where here 


denotes mere ly the numerical 


value of the should 
take @ to keep the variance of the outlet con- 


centration c' 


lag T. We ask how large we 


below an assigned level. 
Substituting (27) into (21) gives the spectrum 
of the inlet concentration as : 
S. (w) a (28) 
Substituting (26) and (28) into (24) gives the 


variance of the outlet concentration as : 


S. Katz 


(10). In (23) the bar denotes the operation of 
taking the complex conjugate value. From (28), 


(29) 


1+0h 


« for some assigned 


to 


Thus, to ensure that o,1 


between 0 and choose make 


h 


here for 
the 
analytical complications of Section ILL, Suppose 


Similar caleulations can be made 


cascades of mixers without encountering 


for example that in the foregoing problem the 
taken to be 


residence time 


mixing mixers in 
0 2hr, so that 


the transfer function of the system Is, by (11): 


system is two 


cascade, cach of 


L (w) 
(l 


(30) 


rather than (25). Substituting (28) and (30) into 


(24) gives the variance of the outlet stream 
concentration as: 
0 
4 
h 


which of course is always larger than the corres- 


ponding variance (29), 


Vv. A 


To exhibit the scope of the statistical methods 


In Discrete Tims 


conclude with an illustrative 
We imagine that 
v ft* 


ake ribed he we 
stage batch mixing probl m. 
batches of volume 


material, of uniform 


are prepared at times indexed by 
n 2, & 


that the nth batch 


¢(n) moles ft® in some critical component, As 


and has the concentration 
each batch is prepared it is fed into a blending 
tank containing V ft" of past material; the tank 
is thoroughly mixed, and then the same batch 
Let c' (n) moles ft® be 


the key concentration in the batch withdrawn just 


volume v is withdrawn. 


after the nth batch has been fed into the mixer 
and blended. A straightforward material balance 
then shows that c'(n) satisfies the difference 


equation 


c' (n) c' (n 1) t Vfe(n) c' (n)} 


(31) 


: 
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As in the cascade equations of Section II we 
want the solution of this difference equation 
which shows the dependence of c'(n) on c(n) 
in a form uncontaminated by extraneous tran- 
sients due to start up. This is achieved as before 
by imagining the mixing system to have been 
started in the remote past, and solving (31) 


under the initial condition : 
c'(n) bounded as n -» xX (32) 


Under suitable regularity conditions on c (n), 


the mathematical system (31), (32) has the 
unique solution : 


H (i)e(n 


c' (mn) 


i) (33) 


where 


H (i) (34) 


As in the continuous problems, the statistics 


arise in thinking of the concentrations as random 
processes ¢(n,«) for the inlet concentration, 
and ec'(n, «) for the outlet concentration, where 
x is the ensemble index which we may think of 
as ranging over duplicate physical systems 
operating under the same statistical conditions. 


We rewrite (33), (34) for the random processes 


Ge 
c' (n. x) K (i) e(n i, (35) 
— 
with 
| 0 2, 3. 
(i) vo Vy’. (36) 
0.1, 2, 


in analogy with (8) and (9), 


as in Section I the moments of these random 


pr 


m, [e (n, x) | 
— E [e' (n, «)] 
p, (v) [fe (n, x) m.| 
pe (n m, p, (0) (37) 
pp (v) E (n, a) m, 
(n v, Me» |; per (0) 
Poor E [fe(n, x) ; 
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Further we define 


all independent of the time index n. The opera- 
tion E denotes, as before, an averaging over the 
ensemble index «. (These definitions are made 
only for what corresponds to the stationary case 
of Section I). 

Applying (37) to (35) shows first that : 


m m 
c 
as it ought, since by (36) : 
K (i) l 


and second that 


Combining the two equa- 


in analogy with (20). 


tions (38) gives : 

‘ A (i) A (j)p.(v +2 j) 
— 


which, on taking v 0, gives for the variance of 


the outlet concentration c' 


AK (i) AK (j)p.(@—j) (89) 
— 


We can, by analogy with Section IV, introduce 


forms of the spectral densities S appropriate 
to the discrete time processes (ny, c' «). 
These are defined to be the Fourier sums, when 


they exist : 


(A) Pp, (v) e' 
1 = 
yy 


3 
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corresponding to (21), with the inverses : corresponding to (24). We shall not, however, 


use this formulation in the example which 


( (A) "ada follows. 

, To illustrate the application of (39), we con- 
sider the simple case where the ke Vv concentrations 
c(n, a) in the feeds are independent from batch 
to batch, that is, where : 


corresponding to (22). If we introduce further 
a discrete time analogue to the system transfer Equations (39) and (36) then give 


function (10): 


Thus, to ensure that CO. * « for some assigned 
then equations (38) give bn tween O and 0,. choose to make 


L (A) 8. (A) 


L(A) S. (A) 


corresponding to (23), and equation (39) gives : = 
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Abstract —Isobaric vapour-liquid equilibrium data of the ternary system benzene-cyclohexane 


methyl isobutyl ketone were determined at 760 mm Hy total pressure using a modified Colburn 


still, The activity coeflicient data for all the experimental points were correlated satisfactorily by 


Wohl’s three suflix Margules ternary equations with a ternary constant. It was concluded that 


methyl isobutyl ketone might be used as an entrainer in the separation of benzene—cyclohexane 


mixtures into pure components. No ternary azeotrope was observed. 


des isobares vapeur-liquide du systéme ternaire benzene 


Résumé —Les données déquilibre 


cyclohexane-methy1] isobutyl cétone ont été déterminées a une pression totale de 760 mm Hg en 


utilisant Pappa eil de Colburn moditix Les cocllicients dactivite pour tous les points experi 


mentaux ont été reliés de fagon satisfaisante par les equations 3 aires de MARGULEs a 3 indices 


de Wout, avee une constante ternaire. Loauteur conclu que la méthyl isobutyl cétone peut 


étre utilisée comme entraineur dans la separation des melanges benzéne-cyclohexane en leurs 


composants purs Il na pis observe d azcotrop. termuire 


Zusammenfassung Dic Werte fiir das isobare Dampf-Flissigkeits-Gleichgewicht des ternaren 


Svstems Benzol-Cvclohexan-Methyvlisobutviketon wurden bei 760 Torr Gesamtdruck in einer 


abgeiinderten Colburn-Apparatur bestimmt Die Werte des Aktivititskoeflizienten fiir alle 


Messpunkte wurden befriedigend korreliert durch ternire Margulesgleichungen vom Drei-Suflix- 


Typ nach Wohl mit einer terniren Konstanten. Ks wurde geschlossen, dass Methylisobutylketon 


als Schleppmittel zur Zerlegung von Benzol-Cyclohexan-Mischungen in die reinen Komponenten 


geeignet ist. Kein terniires \zeotrop wurde beobachtet 


In the purification of benzene from small amounts isobutyl ketone as a solvent. As the vapour 


of higher non-aromatic hydrocarbons a lower liquid equilibria of all the constituent binaries 


non-aromatic hvdrocarbon such aus lol benzene lohexane [1], benzene methyl] iso- 


[2], and  cyclohexane—methyl 
(2) are available in the literature, 


which forms an azeotrope with benzene is added butyl ketone 


and this azeotropic mixture is removed as an isobutyl ketone 
overhead product, while the higher boiling the ternary vapour-liquid equilibrium data of 


impurities are left as residue. The hydrocarbon — the system benzene—cyclohexane—methyl tsobuty! 


products obtained at the top of the column can ketone were obtained at 760 mm total pressure. 


then be fractionated into pure compo nts adopt- 
tech- Materials used 


ing azeotropic or extractive distillation 


niques, using a polar solvent such as alcohol or Benzene: Baker analysed reagent grade 
Letone. benzene was fractionated in a glass column 


and the fraction distilling over at 80-10°C 


Design of the plant for separation of this type 


requires prior knowledge of vapour liquid was collected and used. 

equilibria of the mixtures encountered. This Cyclohexane : Cyclohexane supplied by Light 
investigation was undertaken with a view to & Co. was subjected fractional distillation 
studying the separation of benzene cyclohexane and the fraction distilling over at 80-6°C 


mixtures into pure components using methyl was collected and used. 
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Table 1. 
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Physical properties 


Density at: 
Chemical 


Lit, Obs, 


Benzene [8 40 


Cyclohexane 
Methyl isobutyl ketone 


IS 


Refractive index at 30°C Boiling point at 760 mm 


Lit. Obs. Lit. Obs. 


1 


14050 
14210 


13012 


1 


SO-7OS8 1S) 80-6 


[21 


Methyl isobutyl ketone 
fractionated and 
at 115-9 


Laboratory reagent 


grade methyl tsobutyl ketone was 


the 


was collected and used. 


The physical prope rties of the ch micals used 


in this investigation are compared with the 


literature values and are given in Table 1. 


ANALYTICAL MetTuODS D 


Both refractive index and specific gravity were 
used to fix the composition of a te rnary mixture. 
In brief the method consists of mixing known 
quantities of the three components in which the 
mole ratio of two components was held constant 
\ 


physical 


and determining the two physical properties. 
50 ol t hue 


s. the composition of the component, 


large scale plot, 50 em, 


‘ 


property 


whose mole number is at variance with the mok 
ratio of the other two compone nts as parameter, 
resulting in a family of curves, was prepared, 

The physical properties were next smoothed at 
even values of the composition and cross plotted 
on a 50 em side tsosceles triangl this chart was 
used to analyse the unknown ternary mixtures, 
To test the of method 


several known mixtures wer analysed and the 


validity this analytical 


results showed an average deviation of 1 per cent, 
The at 
30 0-1 C for sodium light using Abbe’s precision 
refractometer, of 0-0005, 
Water was circulated through the prism jackets 


refractive indices were determined 


capable reading to 
maintamed at 30 C by Fisher's Unitized constant 
temperature bath. 


The equilibrium still used and th technique 


of its operation adopted are the same as re port d 


reviously The vapour liquid equilibrium 


fraction distilling over 


data and the values of experimental and cal 
culated activity coeflicients are given in Table 2. 


The smoothed analytical data are shown in Fig. 1. 


Lines of constont ref. index 
Lines of constant spe. gravity 


benzene 


mole 


Cyclonexane methyl isobutyl ketone 
mole / MLB. K. 
Fic. 1. 


yravity 


Smoothed data of refractive index and specific 
at the 


methyl isobutyl ketone. 


v. composition in svstem benzene 


cyclohexane 


COEFFICIENTS 


be 
critically evaluated readily and extended by the 
When the 


phase is considered ideal the activity coeflicient 


LATION OF ACTIVITY 


Vapour licpuid equilibrium data 


use of activity coelflicients, vapour 


of a component in a mixture of » components is 
given by 


For cases when the deviations of the o: 


is phase 


could not be ignored, Benepicr et al. [4] expressed 
the activity coellicient as : 


: 
Benzene vol 
100+ 9 
~ 
1958/ 
.0°85 
60 
144 
0-81. 
0-78 
~ 
70 


Ternary vapour-liquid equilibria in the system : 


Table 2. 
total p 


Vapour-liquid equilibria in the system benzene cyclohexane methyl isobutyl ketone at 760 


benzene—« yclohexane—methyl isobutyl ketone 


mm Hg. 


ressure 


Temp. 
S. No. 
95-2 
O5 
SS 


‘1 "1 
Obes. ‘ale. Obs. 
O-OV515 

Oo70 

0200 

O10 

O491 

9729 

Oslo 

O550 

0410 

0510 

O130 

9620 

00350 

9790 

9407 

9917 

1560 

0510 

280 

OS00 

7550 

1270 

6180 

3480 

1580 
2510 
1800 
1940 
2630 
1740 

2590 

9725 

1050 

0670 

$220) 

3410 

0760 


O67 
O28 
003 
07 
003 
O42 
Ol4 
026 
OST 
131 
103 
120 
O72 
O28 
OOS 
ool 
996 


O41 
O25 
O66 
O19 
055 
OOS 
O28 
O45 
O72 


O58 


99903 
0020 
0020 
O830 
370 
1200 
OTM) 
O500 
9900 
9915 
99002 
0030 
OOS80 
OO80 
(280 
0510 
9938 
0100 
0516 
0960 
S770 
6270 
+460 
£350 
4170 
2540 
2800 
1950 
0560 
0510 
(960 
1200 
‘1700 
1460 
1060 
$510 
1050 
1090 
‘1900 


Osu 
O54 
OGS 
ool 
21 
O37 
O55 
OOS 
ool 
167 


147 
O75 
O36 
Oso 
oO 
mol 
O56 
O57 


O75 
o70 
O45 
O76 
ool 
O46 
O61 
O51 
00S 


1-6 
1-71 
1-33 
1-2 
l 
l 
1 
l 
l 


{ P, | 
2-303 RT 


The terms p,, v;, 8; are all functions of tempera- 
ture and values of these were found at various 


temperatures between 80 and 116 C, Consequently 


it was possible to express Z directly in terms of 


temperature. The second virial coefficient B was 


calculated using the equation given by Wout [5] : 


The critical constants for benzene and cyclo- 


RT 


0-146 


0-012 T, 
[327 


(0-197 


hexane are obtained from Perry [6] and those of 
methyl isobutyl ketone are calculated by the 


Y3 
72 41 Cale. 
2-000 150 67-25 4°25 l 
38-75 1-75 62-00 00-9970 
es 37-25 1-50 61-00 6-50 0-9972 
5-75 $25 67-75 25 0 
44-25 6-00 64-75 12-75 0 
: $1-25 12-50 54-75 25-25 
& a7 38-50 174M) 10-50 35°25 
11 2425 3-75 20-00 50-00) 
iz St 20-50 4-75 24-75 61-75 
13 85-5 19-25 M75 24-75 57-75 
S64 37-75 2400) 57-75 
VOL, 16 100-0 10-25 9-25 19-75 32°25 0 
17 93-2 31-400 38-50 20-50 ( 
52/59 Is 90-0 5-00 65-0 
20 57-50 575 74-75 10-50 
21 86-6 62-75 3-75 77-25 70 
22 66-00 3-75 7°75 
23 85-5 56-00 O75 10-50 
ee 24 84-9 $2.25 23-00 50-00 36°25 
25 91-2 17-25 250) 2525 ( 
26 00-4 37-50 6-00 67-50 r 
27 86-4 7:75 19-06 10-50 71-75 
28 86-5 3-25 53-25 5) 76-25 
32 80-8 23-00 61-75 27-00 66-25 
35 34-50 25 54°25 | 
36 S140) $2-75 1-75 
37 1-75 30-75 meow 
3S 32-75 525 
30 88-00) 5-00) 7-50 
80 82-25 10-00 85-50 11-75 
41 71-25 16-25 Ta 2540 
65-00 2450 67-25 20-00 
61-50 22-50 63-75 
5 7% 47-000 15-00 10-25 
48 87-50 27-00 58-75 35-75 0 
P 4, 
Yi 
where 
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method given by Houcex, Watson and Racarz 
7 The data aol 
ponents 


butvl ketone 


| Vapour pressure three com 


benzene, cyclohexane il methyl iso 


ure taken from the literature. 


woo 7 


was drawn for 


For each component log P v, was plotted 


and a mean straight line each 


case ; these lines ar represented by the following 
equations, 


1685-7266 
Benzene [8] log P 
(/ 273) 


1203-526 


Cyclohexan [9] log P OSELOS 

Methyl isobutyl ketone 
[10] log P 
273 
CORRELATION AND Discus or THI 
DATA 

Among the various equations [I], to 15 


developed for correlating ternary 


Won 


pour liquid 


equilibrium data those of IZ) are com 


prehensive and find wide application. The 
appheability of the t rnary Wont [12], Li and 
15 or Wuer {11 equations is) ore 


stricted by the condition : 


which is not satisfied in the present case, Hene 
an atte mip was made to correlate t} lata, with 
Wout’s 12) three suflix Margules equations, 


which utilize six binary constants and on 


rnary constant. 


The end values of the system benzene—cyelo 
hexane were obtained by correlatin data of 
Ricuarps and Harcerraves [1 Wout's 
two sullix equations; the other tw have been 
reported earlier [2] and are given below 


158 


ition C was a immed to be 
ilated : 


obtained between the 


As a first approxin 
zero and the activity coeflicients wer 
good agreement was not 
riment der to 


abxnut 


calculated and exyx 
obtanu more satisia 


13 values of ¢ 


were computed from experimental 


M. 


JAVYARAMA Rao, 


hia. 2 


free The 


V. R. Acharya and C, Venkata Rao 

ternary measurements using appropriate equations 
and compositions as suggested by Wom [12]; 
0-26, 
into the 
sullix Margules equations of Wou [12] 


the following equations are obtained : 


thev were averaged to vive a value of 


By inserting these constants 


three 


seven 


log Y1 O.15S ty | 0-25 7, 0-627 
0-26 (1 27,) 
OL re (0-1 O-G2 7,) r, (0-158 0-568 7.) 
0-267, 7, (1 22.) 
log y, (0-41 O-G2 Vy) (0-252 — 0-62 


7, (1 22.) 


In Table thy 


activity coellicient data for each experimental 


experimental and calculated 


run are compared ; satisfactory agreement was 


obtained in most case s. The ave rave error between 
caleulated and « xperime ntal activ ity coetlicients 


and Yo are respectively 2-4 per cent, 
t per cent and 4:5 per cent and these errors ar 
within reasonable limits of « xperimental ace uracy, 


Lene the thre« sullix M 


represent the data fairly well with the 


rovules equations 
ternary 
constant ¢ 


No ternary 


pre sent svstem. 


azcotrope was observed in the 


a 
x 
< 
4 


qid(M.1.8.K 


Mook 


benzene in liquid (methyl isobutyl ketone 


elect of methyl isobutyl ketone on 


benzene-cyclohexane vapour-liquid equilibria 


VOI 
he 9 
: 
| 
| 
¢ 
e 
+ 
d 
Pond 
4c 60 80 100 
roraneter; mole benzene m 
free basis) : 
al 
— 
7. 


Ternary vapour-liquid equilibria in the system: benz ne—cy lohexane 


The experimental — ternary vapour-liquid constituent binaries as an 


equilibrium data were smoothed and plotte dona reliable establishment « 


methyl isobutyl ketone 


approximation, for the 


vapour—liquid equili 


to 


a few ternary 


account for 


bye used 


ol benzene 


components 


Mmiponent 


fer 


ais 


th 


by earl r works rs 


No ternary az otrope was observed. 


an 


cyclo 


by 


of the activity coefficient 


pure components at 


ne svstem 


of eon ponent 


i f 


mmponent in liquid 


mponent in vapour 


solvent free basis. that is, on a basis free from brium relationships at least 
methyl isobutyl ketone. Fig. 2 shows that the measurements are necessary 
benzene cyclohexane azcotrope was completely ternary effects. as observed 
eliminated when the solvent concentration was 2. 16, to 20 
about OO px reent. Thus the solvent depre sses thy Methy!] isobutvl ketone might 
relative volatility of benzene and increases that entrainer ins the separation 
of cyclohexane, thus causing a s« paration to take hexam mixtures into pure 
place. However further increase of the solvent extractive distillation. 
concentration did not the relative 
NOTATION 
volatility of the binary mixture to any consider 
able extent. a the end values 
curves 
CONCLUSIONS ( ternary « stant 
Ternary Vapour liquid equilibrium data of the iral coe! it of 
svstem benzene cyclohexane methyl isobutyl) total pressure of tl 
ketone were determined at 760 mm He. total 
pressure using a modified Colburn still The 
activity coctlicient data of all the experimental / ihsolute temperature 
runs were correlated satisfactorily by Wounur’s wtivity coeflicient 
three sullix Maroules equations employing six the system 
binary constants and one ternary constant. It mane Ernction of 
" mole fract 
was found that there is a ternary effect and off 
prediction of ternary data from binaries alone factor introduced t 
would lead to ¢ rror, phase deviations 


Although ternary data can be predicted from Subscripts, 1, 2 
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A note on mean residence-times in steady flows of arbitrary 
complexity 


D. B. SPALDING 


limperial College 


of Science and Technology, London 


(Received 24 January 1958) 


Abstract— The 
arbitrary complexity is shown t 


planes divided by the volumetric flow rate, 


upstream of the injection plane 


position but not with concentration or time. 


mean residence-time of a tracer injected into a steady incompressible flow of 
be equal to the vessel volume enclosed between inlet and outlet 
provided that the tracer is prevented from diffusing 
Velocities and effective diffusion coeflicients may vary with 


If the tracer is not prevented from diffusing upstream of the injector plane, the mean residence- 


time exceeds the above by an al 


the inlet 
auteur 


montre que le 


re che complenxite 


rapport du volume du récipient 


volumique 


arbitraire 


uint equal to the effective longitudinal diffusion coeflicient in 
divided by the square of the velocity in that pipe. 


temps de contact moyen d'un traceur introduit dans un 


lors de écoulement permanent, est égal au 


compris entre les plans d’entrée et de sortie et de la vitesse 


a condition dempécher la diffusion en amont du plan d injection, 


Les vitesses et coeflicients effectifs de diffusion peuvent varier avec la position, mais ne varient 


pas en fonction de la concentration ou du temps. 


| ke 


decrite précédemment dune quantit 


la tubulure d’entrée divisé par la carré cde 


Zusammenfassung Die mittlere 
belibiger An 
Kintritts — und 


pressible Stromung 
tervolumen zwischen 
nicht von det 


Liffuss 


oder der Zeit variieren., 


vesetzt lass der Indikator 
windigkeiten und effektive 


Konzentration 


Wenn der 


Verweilzeit 


Indikator von der 


vrosser als dix 


Verweilzeit eines Indikators 
rdnung injiziert wird, ist, wit 
\ust rit tse bene 


Injektionse bene 


Injektionsstelle 


oben angegebene, 


traceur diffuse en amont du plan d injection, la durée de contact est supérieure a celle 
ite égale au coeflicient effectif de diffusion longitudinale dans 


la vitesse dans ce tube. 


der in cine stationire, inkom- 
yleich dem Behal- 
geteilt durch den Volumenstrom, voraus- 


yezeigt wird 


stromaufwirts diffundiert. Gesch- 


nskoeflizienten kOnnen von Ort zu Ort, nicht aber mit der 


stromaufwirts diffundieren kann, wird die 
und zwar um einen Betrag gleich dem 


effektiven Langsdiffusionskoeflizienten im Ejintrittsrohr, geteilt durch das Quadrat der Gesch- 


windigkeit in diesem Rohr. 


Nort! 


PURPOSE oO} 


If a pulse of tracer material is injected uniformly 


across the inlet duct to a vessel of volume V, 


through which incompressible fluid tlows steadily, 
and the measured 


concentration of tracer 1s 


vs. time at the outlet section, a simple formula 
[1] for the 


instant of injection to instant of appearance at 


mean residence-time @,, (time from 


exit of an clement of tracer, averaged over all 


tracer clements) 


0., = V/v (1) 

where v is the volumetric rate of fluid flow. 
This the 
which this formula is valid. It 


conditions under 
will 


that it holds if the tracer material is prevented 


note examines 


be shown 


from diffusing upstream at the injector plane, 
provided that the effective diffusion coefficients 
and velocities are independent of concentration, 
even for flows of 


arbitrary complexity. If 


upstream diffusion at the injector plane can occur, 


we: 
. 
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the formula underestimates the mean residence- 3. ANALYSIS 


time. Equation of conservation 


It will be supposed that the tracer trave - 


2. STATEMENT OF ror INgecror PLANE 
convection and diffusion. Then its concer 


IMPERMEABLE TO TRACER 


¢ at any point obeys the equation 


Consider first the vessel shown in Fig. 1. At 
time @ = 0 a small tracer quantity S is suddenly d¢/)0 = V(DVe) — Vige) 3 
injected uniformly into the fluid crossing the inlet 


where q \s the local velocity vector. and D 


plane 1. A 


sieve’ prevents the tracer from 


local effective (laminar and turbulent) diffus 


diffusing upstream. An analysis device measures 


Both gq and D vary with position, but not wv 


the concentration of tracer ¢c, at the outlet plane 2. 


concentration or time. 


The volume of vessel and duct between the two 


planes is V. The flow velocities can be regarded 


as uniform across each of the planes 1 and 2. Theorem | 


Integrating the equation from @ 0 


-+4-2 0 ©, (3B) becomes 


Now at all points, « 0: for at long er 


times no tracer material remains. Also. for po 


Fic. 1. The flow system considered, 


downstream of — th injector plane 0 


for the tracer takes a finite time to reac! ul 


points in finite concentrations. The L.H.S 


The analysis device gives a result such as is 


of { +) Is therefore ZCTO, 


illustrated in Fig. 2. We have to find the value 


For convenience, put 


Then (4) becomes 


V(DVT) Viql) 0 6 


be 


This 


diffusion equation in the absence of sources 


will recognised as the steady-state 


Without formal analysis (which could of cours 


be given), it is easy to see that, since gradient 


Fic. 2. Typical concentration time diagram at outlet. 


of c, and therefore of 7, are zero at the vess« 


walls, the only possible solution is 


of the mean residence-time which, in accordance oe a 
with the above definition, is given by eile steele 
a 
| Cy 0 dg Equation (7) represents the first theorem, 
0 o (2) which is obvious enough. It means among other 
' . “’ things that, if a steady stream (train of equal 
| c, dé pulses) of tracer is injected, the concentration is 


uniform throughout the vessel. 


\ 
ation 
1 
| = | —Vi(q| edo) (4 
o 
| 
nt 
I cd0 
4 
C, | 
| | | 
| 
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Theorem 2 


To approach @,,, we first multiply 


then 


equation 


(3) by @, and integrate from @ 0 to 


©, obtaining 


The may written 


svmbo the 


N.B. 


square 


In orc | the 


bracket sha represent zere t the 


upper 
finite 

Putting. for 


convenicnce, 


equation is now becomes 


V(DVI)— Viqd (i 


Equation (11) will be 


recoonimscd is 


a disti 


he Ste 


state wit! ted source 


adilfusion equatior 


term /, which Theorem 1 has she 


constant magnituce solution te equation 


an again immedia without 


analvsis which of course uld agaim 


Thinking of J as 


enthalpy and J as a heat source 


formal 


easily be given). epresenting 

that 
wall gradients of c. and J, are is 
evident that the problem is akin to that of thud 
through a heater. Since the total “ heat 
is JV, it is clear that th 


and J at 


flow 


input rate difference 


between J at outlet inlet is given by 


(J J, I \ (12) 


Now 


upstream . 


take 


vidently 


for clarity a plane just 


zero. ‘Therefore 


| (13) 


Reverting to the earlier notation, this can be 


written as 


SPALDING 


c, dd (14) 


Introducing @,, from (2), we have therefore 
proved that 


0 (15) 


m 


which is our * simple formula ” (1). 

Theorem 2. \ different 
prool has previously been given by DANCKWERTS 
jl}. It that the holds 
irrespective of the pattern in the vessel. 
There 


the vessel 


This is somewhat 


will be noted theorem 
flow 
regions of recirculation in 


may even bn 


where the streamlines are closed. 


Changes in flow pattern will introduce changes 


in the shape of the curve of Fig. 2, but we 


have 
not had to make any assumption about this shape, 
with the 


being concerned only position of the 


centroid of the area beneath it. 

Tur ¢ 1 PERM! 

We now carrv out 

for the « in 

inlet 


ABLI 


the corr sponding analysis 


which the ‘sieve’ is absent 


from the duct, so that tracer material is 


free to diffuse upstream, 


Distribution al di 


No change in the 
Fheorem 1 for the region down 
There I Ss still holds. We 


distribution of J in 


analysis is needed in the 
derivation of 
stream of plane 1, 
now have to consider the 
the inlet 

Let the of the 
flow velocity isv A 
duct 


equation (6) for this region runs 


duct. 
duct be A,, so that the 
Let the diffusion coetlicient 


area 
the uniform value dD. 


in the have 


D v dl 0 


16 
da A, da 


where a distance along the inlet duct in the 
stream direction, 

The solution of this equation, with the boundary 
conditions S/v ata 
0 at a 


, (inlet plane), and 


(far upstream), 1s 


v) exp (2 — D, (17) 


. . 
. 

: 
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This means that J falls off exponentially in an 
upstream direction. 


Determination of 0,,: Theorem 3 


The previous analyses leading to (11) is still 


valid; the solution for J, will be seen to differ 
however, For ease of physical understanding, we 
will continue to regard J as a heat input term ; 
no rigour is lost thereby. 


The total “ heat 


Input rate now contains 
two terms, one for the vessel as before. namely 
SV v, and another for the inlet duct. The latter 
must be written as i / A, de, Introducing 
equation (17), valid in this region, we have 

1, da SA, {exp (a ',) D,A,)da 

SD, AS (18) 

The solution of (11) is therefore 
J, =SV/v + SD,A2 (19) 
Introducing the original notation as before. 


a new expression for the mean residence-time is 


obtained, namely 


/ ) 


m 


(V0) 


This is Theorem 3. The second term may be 


interpreted as a correction to the simple 


formula,” necessary because some of the tracer 
material lengthens its residence-time by diffusing 
hye fore 


DAS 


upstream being downstream. 


Alternative ly 


swept 
could he treated as a 
added to the 


correction which must be vessel 


volume V to account for the permeability ot the 


inlet plane. 


5. 


SS TON 


The result contained in equation (20) mav be 


regarded as a generalisation to flows of arbitrary 


compl xity ot a result reported recently by 


Levensrie. and Ssirn [2]. for the case of 


constant-area flow with a uniform diffusion 


coetlicient, 


However. the suggestion of these authors 


that, in complex flows, an average value of 


diffusion coetlicient should used in the 


residence-time correction. appears to be ill 
founded ; it is the diffusion coefficient in the 
inlet pipe which matters. This will often be 
much smaller than the average, 

It also appears that the time correction 


increases as the of the This 


effect of an Imperme able injector 


square intake area. 


means that the 
plane can always by 


obtained by reducing the 


cross-section of the inlet duet In proportion to 
the other dimensions of the vessel, In many 
practh al cases therefore the correction will be 


entire ly negligible. 
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Local residence-times in continuous-flow systems 


P 
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V. DANCKWERTS 


(Received 24 January 1958) 


Summary — If a tracer is injected at the inlet of a continuous-flow system in the form of a step 


r pulse signal, the course of the subsequent change in concentration at a point within the system, 


be used to determine quantitatively the average “ age’ as an indication of the degree 


ynation of the material at this point. 


Résumé d'un traceur, Valimentation dun systéme en écoulement continu, 


us la forme dun signal discontinu ou dun signal dimpulsion permet de déterminer quantita- 


tivement la vitesse de remplacement ou le degré de stagnation de la matiére, en un point du 


par suite des fluctuations de concentration en ce point. 


Zusammenfassung Wenn ein Indikator stossweise in den Eintritt eines kontinuierlichen 


Str ungssystems injiziert wird, kann man den Verlauf der nachfolgenden Konzentration- 
lerung an einem Punkt innerhalb des Systems zur quantitativen Bestimmung der Verdriin- 


eschwindigkeit oder der Aufenthaltsdauer des Stoffes in diesem Punkt verwenden. 


EisewHere in this issue (p. 74), Spacpinc shows material at the point of measurement could be 


that if pulse consisting of a quantity of tracer determined by evaluating graphically the integral 


ected at the inlet of a steady-state flow 


0, then the quantity 


In practice this may be hard to evaluate because 
of uncertainties about the tail of the curve. For 


ere ¢ is the local concentration of tracer at many purposes the median age, 7, will be just 


time ¢ after injection) will have the same value as useful as an index of the local rate of replace- 


it every point in the system, This suggests a ment, Half the material at the point has an age 


ethod of determining the rate of replacement greater than 7; hence 


or degree of stagnation of the material in different 


parts ol the svstem 
c -dt 
SPALDING's proposition is discussed further in 2v 
the Appendix, where it ts shown that the distribu- ° 
tion of aves’ of the material in a volume- The value of 7+ is easily found by plotting ¢ 


element around the point at which ¢ is measured against ¢ until the area under the curve is equal 


can be deduced from the variation of ¢ with t. to @/2v. 


The “age” of a particle of material is the The same information can be obtained by 


time elapsed since it entered the system). The raising the concentration of tracer in the inlet 


fraction of the material having ages between — stream at f 0 from zero to a steady value C. 


, (f+ dt) is (ve Q). dt, where v is the rate of The value of ¢/C at the point of measurement at 


flow of material through the system, and @ the — time ¢ is then equal to the fraction of the material 


quantity of tracer injected. having an age less than ¢, and the median age 
It follows that the average age, 0, of the is equal to the value of ¢ at which ¢/C reaches }. 


: 
= 
9 
Se, 
1958/ 
stem at faa 
ct, at, 
0 
| dt 
= 
~ 
2 


Local residence-times in continuous-flow systems 


These techniques offer a quantitative method 
of investigating rates of replacement of material 
in various of a continuous-flow 


parts plant 


(or of air in a room or water in a pond). For 
instance, a contour map of median ages through- 
out the plant might be made, and would show 
up clearly dead spots and by-passing effects. In 
general one would expect the median or average 
ages near the inlet to approach zero, and their 
values near the outlet to approach those for the 
material in the whole stream leaving the system. 
the 
be considerably greater than 


age of 


However, in stagnant pockets the 
material might 
that of the leaving stream, and channelling or 
by-passing effects would be shown up as regions 


where the age was relatively low. 


APPENDIX 


A quantity @ of tracer is injected into the inflowing 


stream. If Q is finite, the injection cannot be 


instantaneous : it will be 


truly 
assumed that it occurs at a 
time-interval Let the 


distribution function for the ages (@) of the particles of 


constant rate over a short 


material in the volume element, which contains a quantity 
5 m of material be f(@). At any time ¢ after injection the 
age of the tracer will be between ¢, (1 57), and hence 
the fraction of the material in the volume-element having 


the same age as the tracer is f(t). d¢ 


Since a quantity 
v.4¢ entered the system with the pulse of tracer, and 
therefore has the same age, the ratio of tracer to other 
Provided 


therefore that tracer and other material migrate tovether 


material entering during the pulse was Q/v . 51 


the quantity of tracer.in the specified volume-element 


at time ¢ is 


( A ( t) 
f(t). dt. bm 


D.¢ 


i 


and hence the ratio c, of tracer to the whole of the 


material in the volume element is given by 


Q f(t) 
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Lateral diffusion with liquid flow through a packed bed of ion-exchange particles 


(Received 6 March 1958) 


INTRODUCTION >. and of the radial distance from the 


uXts, rT, Is approxi 
with the ctio mately given by 


Some difficulties encountered 


the distribution of feed streams into 


(1) 
ion-exchange resin particles led us t it mm 
ipparent radial svitics mm 


In this In the of this equation it has been assumed that 
of a solution conta iNtal transport | nevhuwible and that the 
fixed bed of known iintainmed radial diffusion infinitely thin line 


In 


with the same com «ol 1 int ier lilly in wtual 
ditronal ount of ton ; i through experiment the mi ei stream has a finite cross-sectional 
the aXts ed and 


so that int er eivhbourheod of the injection 

t hee The aver uld valid Further, the bed 
is made equal to the supertic treat wlially confined. i limdrical with a radius 
After a sufliciently long wever the reg contaming the bulk of the diffusing 
conditions the ton com iu er peartiod iterial ith respect to the dimension of the in 
in equilibrium with tl teat ll with respect to the bed diameter 
n of the latter is i } equation (1) ts ery good approximation of the actual 


mentration distritution to the 
expertments to te 


' i . leseribed, equation (1) will be used 
f inert basis for further caleulation Finalls 
valid only for the ise where the approach velocities of 

the ‘ ‘ breve } 


the myected hqunl and of the maim stream are equal 
of thy ne 


this 


md under equilibriun 


equation (1) ts 
etri structure “ the 


medition was obeyed in the experiments 
letermined b u with a noent ted solution Tre the present me of the exchange of univalent ions 
img « ind by subsequent titration his with a monofunctional resin the equilibrium relationship 


otal amount ts a measure of the spreading the myected between the ij 


m compositions in the liquid and in the 


liquid and theref an expressed as 


re of the yparent ra liffusivity resin ¢ 


Phe measurements reported here were px rmed with 


narrow size fraction of spherical 
catron-exchanyge particles with an iver liameter 
d O75 mm main fl 


of Natl in water 


w comsisted solution 


(concentration where and 


are the relative concentrations or fractions 
contamed the same ar unt elati of the displacing 


ions in the liquid and resin phase respec 
resin tively (see 


high excess concentration 
originally charged with Na* o e of the With 
total concentration near the injection point the resin 


particles would tend to shrink somewhat n the other 


the displacement of Na* by It the parti les 
to swell The above conditions were chose n order to and equation (2) we may write for 
pra les so that 


no local high porosities and channelling was ¢ 


" 
secure a small net swelling effect of the 
expected 
fie (3) 
he relationship between the total amount of H*-ions 
in the resin phase and the radial diffusivit DD 


given in the next section. 


will 


The total amount of H*-ions present in the resin price 


in the length interval 0 ‘ L is given by 
THEORY 


Under the steady state conditions described above. 


the excess, Ac, in H*-ion concentration in the liquid as 


a function of the axial distance from the injection point, 
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where q¢ is the capacity of the ion-exchanger per unit 


volume of bed 


After substitution of equations (1) and (3) into (4) the 


latter can be integrated to give 


my 


with 


Ke, DL 


\ 
Ae 


When A changes from 0 to &, the expression between 


brackets in equation (5) varies from 0 to 2. If the conditions 


of the « xperiment Land v) and the properties 


of the system (q and A) are known equations (5) and (6) 


can be used for the calculation of D_ from the experimen 


tally determined value of may 


Apart from the experimental requirements connected 


with the validity of equation (1), the method postulates 


that the spreading of the injected excess concentration 


can be described by means of a constant radial diffusion 


coellicient 


The values of the various quantities for these experi 


ments were 


i, 10 (inner diameter of injection 
tubs 2-07 mm) 

~ 0-010 kg-equiv, Na*-ions /m® (except for run no. 10 
where Ce O-10) 

Ac, 1) equiv. H*-ions 

A 1-7 (determined s« parately, see also [2}) 

q 21 kg-equiv./m? (also determined separately) 

‘ 40 v) 


The other experimental conditions (the tube diameter 


2H, L and v,) are given in Table 1, together with the 


Lateral diffusion with liquid flow through a packed bed of ion-exchange particles 


Table 


m,,-values found and the values of A and v /D, calculated 


from the data by means of equations (5) and (6). The 
values of D. obtained vary between 0-4 and 3-3 107 
m* /sec, ie. a few orders of magnitude greater than the 
coeflicients of molecular diffusion ot electrolytes in water. 
Since radial diffusion in packed beds is frequently expressed 
in terms of the dimensionless group vd,,/D, (a Péclet 
number), the value of this number is indicated in Table 1 
together with Re), the Reynolds number based on the 
particle diameter and the superficial flow velocity. 

It is seen from Table 1 that the values of vd,,/D, are 
fairly constant and that no significant dependence on 
Reynolds number exists. Runs no. 3, 8, 10 and 11. with 
practically equal results 


for vd, D,, although 3 and 8 differ in tube diameter. 8 


nearly equal values of 1 


and 10 in the value of €@ and 8 and 11 in the length L. 
This indicates that very probably the spreading of the 
excess concentration may be treated as a diffusional 
process with constant diffusivity and it justifies the 
application of equation (1) to these experiments. The 


latter can also be checked by more detailed calculations 


A comparison of these results with other work is of 
interest and Witne [3] measured lateral 
diffusivities in packed beds in the region 5 Re 2.000 


The value of vd, D proved to be of the order of 10. 
but variations were observed between 5 and 20. Later on. 
LATINEN [4, 5] continued this work his experiments 
showed that vd,/D, had a constant value of about 11 
for It ; 200, but that at lower Re-numbers this value 
increased The latter results have been represented in 
Fig. 1. Latinen’s explanation for this effect, ie. that as 
the flow becomes more laminar, the injected liquid is 
less thoroughly mixed betwee nm the parti les and becomes 
more striated with a preference to its original lateral 
position, is entirely plausible ; this was also borne out by 
visual observation with a model. It is, therefore. sur- 


prising that the average value of 1 d,, D.,. of about 26, 


Temp L t 
No. (10°73 m) (m) 3 m 


(10°° ke equiv. — 


v/D d, D v, d,, 


~ 


= 
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ResuLtTs 
(10M, — — 
sec 
17 140 34-8 340 23 1-05 
2 19 37-5 0-229 1-09 33-0 O44 3-3 25 1-46 
3 17 37-5 2-09 33-0 O45 3-4 26 2-08 
4 14 62-2 0-210 0-62 28-3 0-39 3-4 26 0-40 
5 16 62-2 O-s2 28-9 Ow 3-35 25 
16 62-2 210 1G 25-7 31 
7 15 62-2 0-209 1-83 26-2 29 
8 14 62-2 0-200 3-07 29-2 0-42 3-2 24 1-06 
14 62-2 O210 1-47 28-2 O38 3-5 26 2-86 
17 62-2 B07 6-35 3-70 3-6 27 2-13 
11 17 62-2 O-B17 B07 51-8 3-5 26 2-13 
| 


. E. Harrman, C. J. H. Wevers and H. Kramers 


volue for Rep > 200 
4 10 40 100 
Yo dp/v 


Lat inen [4] | this paper 
| symbol | 6 o ° . ° 

| [ 0-018 0036 0-067 0-10 0020 0012 


Fic. 1. Measured radial diffusivities compared with Larinen’s results. 


average particle diameter (m) 


obtained with the present experiments in a range of still d 


smaller Re-numbers, is much smaller than would be D. apparent lateral (or radial) diffusivity (m*/sec) 
expected on the basis of Latinen’s results and considera- AK equilibrium constant, equation (2) 
tions. Since, however, we cannot find any reason for which L = height of the packed bed, measured from injection 
the method used here and its results would be entirely point (m) 
3 wrong, we think that a further investigation along these m,;, = total amount of H’-ions bound by the resin 
. lines would be interesting. (kg-equiv.) 
M. E. Hartman q = capacity of the ion-exchanger (kg-equiv./m' of 
Laboratorium voor Fysische Technologie C.J. U. Wevers fixed bed) 
Te hnische Hogeschool, Delft, Netherlands H. Kramers r = radial distance from the axis (m) 
v = average liquid velocity in the bed (m/sec) 
: V_ = liquid velocity based on the empty cross-section 
NOTATION (m/sec) 
A; = cross-sectional area of injection tube (m?*) x = relative concentration of H’-ions in the liquid ( ) 
f, = concentration of Na*-ions in main and injected y = relative concentration of H*-ions in the resin ( ) 
stream (kg-equiv./m*) z = distance downstream from the injection point (m) 
Ac = excess concentration of H’-ions in the liquid « = porosity of the bed ( ) 
(kg-equiv. /m*) Re, = 09 d,/v ( —) 
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J. Zernike: Chemical Phase Theory. Klumer's 
Pub. Co. Ltd., Netherlands 1957, 666 diagrams, 493 pp. 

Tuts book gives a very thorough account of the subject 
which in English is usually referred to as “ The phase 
rule 
in systems containing one or more components. Up 


that is to say it deals with the phase relations 


till very recently the most comprehensive treatment of 
the subject available was probably the series of volumes 
edited by H. W. B. Roozenoom, published in German ; 
however, beginning as it did in 1901, the data given in 
this work are far from up to date. Zrrntke’s book 
contains an account of the subject in some 493 pages 
which, according to the author, covers literature references 
up to 1954, and refers to no less than about 680 actual 
systems in illustration of the theory. The subject matter 
is treated as follows: derivation of the phase rule 18 
pp.; One-component systems 17 pp.; binary systems 
213 pp. ; ternary systems 139 pp.; more complex systems 


39 pp.: metastable systems 26 pp. ; and answers to 


problems 20 pp. 

The data to which reference is made above serve as 
illustrations to the general theoretical discussion of the 
subject : it would appear that the theoretical development 
and discussion is treated throughout in a_ scientifically 
satisfactory manner. On the other hand the English 


style is not above criticism and makes the reading of 


the book at times somewhat diflicult. 


As examples of awkward turns of phrase two sentences 
from p. 11 may be quoted: “ through the membrane 
the carbon dioxide diffuses in this form that two ions 
of HCO, go from right to left, and 
have a very pronounced example that... .the vapour 


we 


pressure ....%.... 


The reviewer also has a minor quarrel with the author 
over the use of certain words: for instance (p. 120) 


“monotonous curve” should surely be “ monotonic 


curve * (p. 63) “have a common tangent with the 
liquidus in the extreme “ is surely apt to be misread, as 
the author is referring to an extreme value; nor does 


he really like the word “ multijugary ” (multi-com- 


ponent’). The author introduces “ singularly” for 


‘one-component and “chaoticus” for gas-phase 


boundary curve these may be useful innovations, 
although * chaoticus has got to compete with vaporus” 


introduced by other authors. 


One of the most interesting chapters is the one on 
metastable states, a somewhat unusual aspect of the 
subject. This deals for instance at considerable length 
with the sulphur system; it also discusses a number of 
other less important systems with metastable phases, 
including what the author calls pseudo-mixed crystals, 
but which are normally described as cases of adsorption 


on the solid phase. 


As pointed out by the author, the text book by Ricc1 
The Phase Rule and Heterogeneous Equilibrium appeared 
whilst Zernike’s work was still in manuscript. The 
two books undoubtedly cover much the same field, but 
in such an all-embracing subject as phase theory there 
is certainly room for two text books in the English 


language. 


To sum up Zernike’s book constitutes a compre- 
hensive and sound scientific survey of the phase theory, 
very fully illustrated with examples from actual systems. 
Some criticism can be levelled against the book on the 
grounds of style, and it may be hoped that this will be 


improved in a future edition. 


R. 
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Am 8. Oktober finden Besichtigungen von Industriebetricben Berlins statt. Naihere Auskiinfte tiber die Tagung erteilt 
die Geschiiftsstelle der VDI-Fachgruppe Vertahrenstechnik, Frankfurt/Main, Rheingau Allee 25. 
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